Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,155,459 A→C 30.1% V8V (GTA→GTC RCAP_rcc02000 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,155,4590A→C30.1% 236.7 ‑0.5 ‑0.2249V8V (GTA→GTCRCAP_rcc02000hypothetical protein
Reads supporting (aligned to +/- strand):  new base (36/39):  ref base (91/83):  total (127/122)
Fisher's exact test strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.11e-01

Only 100 of 266 total aligned reads displayed.

CGCCCTCACCTTCTGATGCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGC  >  NC_014034/2155359‑2155559
                                                                                                    |                                                                                                    
cccccCCACCTTCTGAGCCCCCCGGCCGCCGCCACCGTCCCCGCCCGCCGGGCGGAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTa                                                                                                      <  2:3164335/99‑1 (MQ=255)
  cccTCCCCTTCTGATGCCCTCGGGCGCCGTCCCCGTCACCGCCCGCCGGGCGGAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTaga                                                                                                    <  2:2764786/101‑1 (MQ=255)
  cccTCACCTTCTGATGCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTaga                                                                                                    >  2:2333980/1‑101 (MQ=255)
   ccTCACCTTCTGATGCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTagag                                                                                                   >  2:4411295/1‑101 (MQ=255)
     cccccTTCTGATGCCCTCGGGCGCCGTCCCCGTCCCCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCt                                                                                                 <  2:421335/98‑1 (MQ=255)
        ctttCTGATGCCCGGGGGCGCCGCCACCGTCCCCGCCCGCGGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCtcgt                                                                                              <  2:806474/99‑1 (MQ=255)
          ttcggacccccGCGGGCGCCGCCACCGTCCCCCCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCTTAGAGCtcgtcg                                                                                            <  1:2733105/93‑1 (MQ=255)
          ttCTGAGCCCCTGGGGCGCCGCCCCCGCCCCCCCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCtcgtcg                                                                                            <  2:3633903/101‑1 (MQ=255)
           ttggcggCCCGCGGGCCCCGCCCCCGCCACCCCCCCCCGGGCGGAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGa                                                                                           <  2:2860112/95‑1 (MQ=255)
           tCTGATGCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGa                                                                                           >  2:4308350/1‑101 (MQ=255)
           tCTGATGCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGa                                                                                           >  2:3446505/1‑101 (MQ=255)
           tCTGATCCCCCGGGGCGCCGCCCCCGCCCCCCCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGa                                                                                           <  1:10003/101‑1 (MQ=255)
             tGATGCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATg                                                                                         >  2:4556896/1‑101 (MQ=255)
               acccccTCGGGCCCCGCCCCCGCCCCCCCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGcc                                                                                       <  1:1372810/98‑1 (MQ=255)
               aTGCCCTCGGGCGCCGTCACCGTCCCCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGAGATGCAAGGGGCCTCGTCGAGCTCGCCGATGcc                                                                                       >  1:4348606/1‑101 (MQ=255)
                     tCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTgg                                                                                 >  1:1932655/1‑101 (MQ=255)
                       gggCGCCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTgggt                                                                               >  2:3396800/1‑100 (MQ=255)
                         gcgcCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTggggcg                                                                             >  2:4350691/1‑101 (MQ=255)
                              gtcaccgtcaccGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTggggcgggggc                                                                        >  2:3490405/1‑99 (MQ=255)
                              gtcaccgtcaccGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGcc                                                                        >  1:3099941/1‑101 (MQ=255)
                              gtcaccgtcaccGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCACCTggggcggggcg                                                                        >  2:2471868/1‑100 (MQ=255)
                              gcccccgcccccgcccgcCGGGCGGAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTCGGGCGGGGcc                                                                        <  1:4365658/91‑1 (MQ=255)
                                ccccgcccccgccggcGGGGCGGAAAGCCGCGAAAGGGCGGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGt                                                                      <  2:1846845/93‑1 (MQ=255)
                                   cgtcaccgCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTcac                                                                   >  1:373763/1‑101 (MQ=255)
                                    gcccccccccgccGCGTGGAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTcacc                                                                  <  1:1525600/94‑1 (MQ=255)
                                           cccggcGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATgg                                                           <  2:2101711/101‑1 (MQ=255)
                                            ccgcggGGCGGAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATggc                                                          <  1:2334144/101‑1 (MQ=255)
                                               ccgGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATggcggc                                                       >  1:4565002/1‑101 (MQ=255)
                                               ccgGGCGCCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATggcggc                                                       <  2:3248492/101‑1 (MQ=255)
                                                 gggCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATggcggcgg                                                     >  1:2236545/1‑101 (MQ=255)
                                                      gcAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTc                                                >  1:3475983/1‑101 (MQ=255)
                                                      gcAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTc                                                <  2:4516518/101‑1 (MQ=255)
                                                        aaGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGc                                              <  2:3860977/101‑1 (MQ=255)
                                                         aGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGcc                                             <  1:558785/101‑1 (MQ=255)
                                                          gCCGCGAAAGGGCAGGATCATGAGCGATGCGAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCg                                            <  1:3239468/101‑1 (MQ=255)
                                                              cgAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGggcg                                        >  2:3030400/1‑101 (MQ=255)
                                                               gAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGggcgg                                       >  1:2411906/1‑101 (MQ=255)
                                                                    gagaaaaaaaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGAt                                  <  1:3399939/92‑1 (MQ=255)
                                                                     gCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGAtg                                 >  1:1468239/1‑101 (MQ=255)
                                                                     agaaaaaaaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGAtg                                 >  1:1260311/9‑101 (MQ=255)
                                                                      cAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGAtgg                                >  2:715198/1‑100 (MQ=255)
                                                                       aGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGAtgtg                               >  2:1095847/1‑101 (MQ=255)
                                                                        aaaaaaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTgg                              >  2:2773070/6‑101 (MQ=255)
                                                                         gATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGc                             <  1:2454842/101‑1 (MQ=255)
                                                                          aaaaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGAGGTGTCa                            >  2:3145764/4‑101 (MQ=255)
                                                                          aTCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCa                            >  2:4846322/1‑101 (MQ=255)
                                                                           tCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGGGGCGGCGCTCGCCGGGCGGTTGATGTGGCAt                           >  2:4723196/1‑101 (MQ=255)
                                                                           aaaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCAt                           >  2:582911/3‑101 (MQ=255)
                                                                            cATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCAtt                          <  2:1932518/101‑1 (MQ=255)
                                                                            cATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCAtt                          >  2:1390765/1‑101 (MQ=255)
                                                                            aaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGGGGCAtt                          >  2:1635381/2‑101 (MQ=255)
                                                                             aTGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTc                         <  1:1683617/101‑1 (MQ=255)
                                                                             aTGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTc                         >  1:3650219/1‑101 (MQ=255)
                                                                             aTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTc                         <  1:1839449/101‑1 (MQ=255)
                                                                             aTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTc                         <  1:3505755/101‑1 (MQ=255)
                                                                             aTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTc                         <  2:298538/101‑1 (MQ=255)
                                                                             aTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTc                         <  1:4615716/101‑1 (MQ=255)
                                                                             aTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTc                         <  2:1269636/101‑1 (MQ=255)
                                                                             aTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTc                         <  1:670948/101‑1 (MQ=255)
                                                                              tGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCg                        <  1:4846649/101‑1 (MQ=255)
                                                                              tGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCg                        <  1:4293254/101‑1 (MQ=255)
                                                                               gaccGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGa                       <  2:4778579/98‑1 (MQ=255)
                                                                                aGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGAc                      <  2:373710/101‑1 (MQ=255)
                                                                                 gCGATGCAAGGGGCCTCGGAGAGGTCGTCGAGGCCATCGGGGGCGGGGCCGCCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGAcc                     <  1:1345252/101‑1 (MQ=255)
                                                                                 ccGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGAcc                     <  2:2201768/100‑1 (MQ=255)
                                                                                  cGATGCAAGGGGCCTCGTAGAGCGCGGCGATGCCATCTGGGGCGGGGCCGTCACCACGCTGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCg                    <  2:2577201/101‑1 (MQ=255)
                                                                                  cGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCg                    <  2:428305/101‑1 (MQ=255)
                                                                                    ataccagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAg                  >  1:3310773/9‑101 (MQ=255)
                                                                                    ataccagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAg                  >  1:2688014/9‑101 (MQ=255)
                                                                                    aTGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAg                  >  2:3538856/1‑101 (MQ=255)
                                                                                     taccagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAgg                 >  2:2690260/8‑101 (MQ=255)
                                                                                     taccagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAgg                 >  1:3941633/8‑101 (MQ=255)
                                                                                      gCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGt                <  2:1764130/101‑1 (MQ=255)
                                                                                      accagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGt                >  2:872156/7‑101 (MQ=255)
                                                                                       cAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAgggg               >  2:3956348/1‑99 (MQ=255)
                                                                                       cAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTg               <  1:1976870/101‑1 (MQ=255)
                                                                                        aaGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGc              >  2:4660243/1‑101 (MQ=255)
                                                                                         gggggCCTCGTCGGGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCg             <  1:3051605/100‑1 (MQ=255)
                                                                                          ggggCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCgg            >  1:3829645/1‑101 (MQ=255)
                                                                                          ggggCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCgg            >  2:253168/1‑101 (MQ=255)
                                                                                          gaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCgg            >  2:1078773/3‑101 (MQ=255)
                                                                                          gaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCgg            <  1:3530322/99‑1 (MQ=255)
                                                                                           aggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGGGGGGCggg           >  1:4481366/2‑101 (MQ=255)
                                                                                            ggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgc          <  1:4066761/101‑1 (MQ=255)
                                                                                            ggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgc          <  1:251614/101‑1 (MQ=255)
                                                                                            ggCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGgc          >  1:1099830/1‑101 (MQ=255)
                                                                                            ggCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGgc          <  2:519624/101‑1 (MQ=255)
                                                                                             gCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGAGGATGGCGGCGGGGCGGGCCGGGCGCCTGATGCGGCATTCGACCGAGGTGCGGgcg         >  2:4148862/1‑101 (MQ=255)
                                                                                             gCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcg         >  1:4821429/1‑101 (MQ=255)
                                                                                               cTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcgcg       >  1:49693/1‑101 (MQ=255)
                                                                                                tCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGgcgcgc      <  1:2819493/101‑1 (MQ=255)
                                                                                                 cGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCa     <  2:1478791/101‑1 (MQ=255)
                                                                                                 cGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCa     <  1:4185574/101‑1 (MQ=255)
                                                                                                  gTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCaa    <  1:841870/101‑1 (MQ=255)
                                                                                                  gTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCaa    >  2:2349094/1‑101 (MQ=255)
                                                                                                  gTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCaa    <  2:3586950/101‑1 (MQ=255)
                                                                                                   tAGAGCTCGTCGNNNNCATCNGGGGCGGGGCCGTCNNCACGATGGCGGCGGCGCTCGCCGGNNGGCTGATGTGGCATTCGANCGAGGTGCGGGCGCNCANg   <  2:4927982/101‑1 (MQ=255)
                                                                                                   tAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAg   >  1:2074494/1‑101 (MQ=255)
                                                                                                   tAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAg   >  1:4720245/1‑101 (MQ=255)
                                                                                                   tAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAg   >  2:2937976/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CGCCCTCACCTTCTGATGCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGGCAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGC  >  NC_014034/2155359‑2155559

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG