Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 773,138 G→A 5.3% G33S (GGC→AGC)  cydB → cytochrome d terminal oxidase, subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913773,1380GA 5.3% 555.7 / 12.9 168G33S (GGC→AGC) cydBcytochrome d terminal oxidase, subunit II
Reads supporting (aligned to +/- strand):  ref base G (88/71);  new base A (4/5);  total (92/76)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

Only 100 of 173 total aligned reads displayed.

TGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGC  >  NC_000913/772990‑773286
                                                                                                                                                    |                                                                                                                                                    
 gAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTggg                                                                                                                                                     <  4:1356531/149‑1 (MQ=255)
     aGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATg                                                                                                                                                 >  3:222497/1‑149 (MQ=255)
      gTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTgg                                                                                                                                                      <  3:537676/143‑1 (MQ=255)
      gTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTgg                                                                                                                                                      >  4:537676/1‑143 (MQ=255)
      gTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGc                                                                                                                                                <  2:136608/149‑1 (MQ=255)
      gTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGc                                                                                                                                                <  2:537478/149‑1 (MQ=255)
       tCTTCCACGACTACGCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCt                                                                                                                                               >  1:1111878/1‑149 (MQ=255)
         ttCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCa                                                                                                                                             >  2:188152/1‑149 (MQ=255)
         ttCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCa                                                                                                                                             >  4:128958/1‑149 (MQ=255)
             aCGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGTATGCTCACCCg                                                                                                                                         >  3:946988/1‑149 (MQ=255)
             aCGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATCGGGGTGGGCATGCTCACCCg                                                                                                                                         >  3:1198510/1‑149 (MQ=255)
               gACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGtt                                                                                                                                       >  3:318390/1‑149 (MQ=255)
                actactCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGttt                                                                                                                                      <  1:303724/149‑1 (MQ=255)
                actactCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGttt                                                                                                                                      <  4:1223009/149‑1 (MQ=255)
                 ctactCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTc                                                                                                                                     >  1:1306346/1‑149 (MQ=255)
                 ctactCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTc                                                                                                                                     >  4:119405/1‑149 (MQ=255)
                  tactCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTcc                                                                                                                                    >  1:641408/1‑149 (MQ=255)
                  tactCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTcc                                                                                                                                    >  2:627748/1‑149 (MQ=255)
                  tactCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTcc                                                                                                                                    >  2:677620/1‑149 (MQ=255)
                    ctCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTc                                                                                                                                  >  2:228527/1‑149 (MQ=255)
                     tCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCg                                                                                                                                 <  4:615988/149‑1 (MQ=255)
                        gCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTc                                                                                                                              >  2:764646/1‑149 (MQ=255)
                        gCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTc                                                                                                                              >  2:1010239/1‑149 (MQ=255)
                        gCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTc                                                                                                                              >  1:1466565/1‑149 (MQ=255)
                        gCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCAGTTTCCTCGGTc                                                                                                                              <  3:108553/149‑1 (MQ=255)
                         ccGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCg                                                                                                                             <  3:119405/149‑1 (MQ=255)
                            gCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGAGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTaa                                                                                                                          <  3:1258765/149‑1 (MQ=255)
                                gCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGac                                                                                                                      >  4:841253/1‑149 (MQ=255)
                                gCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGAGCATGCCCAcccc                                                                                                                                         <  4:176869/130‑2 (MQ=255)
                                gCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGAGCATGCCCAcccc                                                                                                                                         >  3:176869/1‑129 (MQ=255)
                                  tAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGacac                                                                                                                    >  2:733056/1‑149 (MQ=255)
                                          gAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTc                                                                                                            <  3:140849/149‑1 (MQ=255)
                                          gAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTc                                                                                                            >  2:265193/1‑149 (MQ=255)
                                          gAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTATTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTc                                                                                                            >  3:1233457/1‑149 (MQ=255)
                                             tcgtcAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGaa                                                                                                         <  2:1111878/149‑1 (MQ=255)
                                             tcgtcAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGaa                                                                                                         <  3:41519/149‑1 (MQ=255)
                                              cgtcAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAAt                                                                                                        <  1:964621/149‑1 (MQ=255)
                                               gtcAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAAtt                                                                                                       >  1:1245402/1‑149 (MQ=255)
                                               gtcAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAAtt                                                                                                       >  1:217249/1‑149 (MQ=255)
                                               gtcAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAAtt                                                                                                       >  3:562834/1‑149 (MQ=255)
                                               gtcAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAAtt                                                                                                       >  2:1072614/1‑149 (MQ=255)
                                                 cAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTAt                                                                                                     <  1:1072614/149‑1 (MQ=255)
                                                 cAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTAt                                                                                                     <  2:826413/149‑1 (MQ=255)
                                                  aaaTGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATg                                                                                                    <  4:584976/149‑1 (MQ=255)
                                                   aaTGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTc                                                                                                                              <  1:519194/122‑1 (MQ=255)
                                                   aaTGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGa                                                                                                   >  3:1164760/1‑149 (MQ=255)
                                                   aaTGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGa                                                                                                   >  2:1205623/1‑149 (MQ=255)
                                                    aTGATCGATTATGGAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGAt                                                                                                  <  4:1127344/149‑1 (MQ=255)
                                                       aTCGATTATGAAGTATTGCGTTTTATCTTGTGGCTGCTGGTTGGCGTTCTGCTGATTAGTGTTGCAGTCACGGACGGTTTCGACATGGGGGTGGGCATGCTAACCCGTTTCCTAGGTCGTAACGACACCGAGCGTCGAATTATGATTaa                                                                                               >  2:422426/1‑149 (MQ=255)
                                                       aTCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTaa                                                                                               >  4:264403/1‑149 (MQ=255)
                                                       aTCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTaa                                                                                               >  2:1424896/1‑149 (MQ=255)
                                                        tCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGGATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAAc                                                                                              >  1:256766/1‑149 (MQ=255)
                                                           aTTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTcc                                                                                           <  1:733056/149‑1 (MQ=255)
                                                           aTTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTc                                                                                            >  1:1180780/1‑148 (MQ=255)
                                                           aTTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTc                                                                                            <  2:1180780/148‑1 (MQ=255)
                                                              aTGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCAtt                                                                                        <  3:128958/149‑1 (MQ=255)
                                                               tGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTg                                                                                       >  3:676568/1‑149 (MQ=255)
                                                               tGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTAGAATTATGAtt                                                                                                 <  3:826432/139‑1 (MQ=255)
                                                                gAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGa                                                                                                                 <  4:726608/122‑1 (MQ=255)
                                                                gAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGa                                                                                                                 >  3:726608/1‑122 (MQ=255)
                                                                 aaGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGca                                                                                     <  4:1233457/149‑1 (MQ=255)
                                                                         cGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTg                                                                             <  1:27951/149‑1 (MQ=255)
                                                                         cGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTg                                                                             <  4:827086/149‑1 (MQ=255)
                                                                           ttgtATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCATACTagg                                                                           >  3:970634/4‑149 (MQ=255)
                                                                               aTCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACgg                                                                       >  1:694272/1‑149 (MQ=255)
                                                                                      ggCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCTTTGCACCACACTGGGACGGTAACCAg                                                                <  2:1417624/149‑1 (MQ=255)
                                                                                      ggCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAg                                                                >  1:716839/1‑149 (MQ=255)
                                                                                      ggCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAg                                                                >  2:1092862/1‑149 (MQ=255)
                                                                                         tgctgGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGtt                                                             >  3:1456835/1‑149 (MQ=255)
                                                                                            tgGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATAAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTgg                                                          >  2:843239/1‑149 (MQ=255)
                                                                                             ggTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGc                                                         <  2:1306346/149‑1 (MQ=255)
                                                                                               ttGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTg                                                       >  3:837014/1‑149 (MQ=255)
                                                                                                tGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGa                                                      <  3:235936/149‑1 (MQ=255)
                                                                                                tGGCGTTCTCCTGATTGGTTTTGCAGTCACGGACGGTTTCCACATGGGGGCGGGAATGCTAACCCGTTTCCTCCGTCGTTTCGACACCGATAGTGGAATTATGATTAACTCCATTGCACTACACTGGGATGGTAAGCAGGTTTGGCTGa                                                      <  2:462587/149‑1 (MQ=255)
                                                                                                 ggCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGAt                                                     <  2:889378/149‑1 (MQ=255)
                                                                                                       ctgctgATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACcgc                                               <  3:662032/149‑1 (MQ=255)
                                                                                                          ctgATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCgcggg                                            <  1:1010239/149‑1 (MQ=255)
                                                                                                           tgATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCgcgggc                                           >  1:1029973/1‑149 (MQ=255)
                                                                                                             aTTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCgcgggcgg                                         >  1:601163/1‑149 (MQ=255)
                                                                                                             aTTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCgcgggcgg                                         <  4:389775/149‑1 (MQ=255)
                                                                                                             aTTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCgcgggcgg                                         <  4:222497/149‑1 (MQ=255)
                                                                                                             aTTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCgcgggcgg                                         >  4:1309386/1‑149 (MQ=255)
                                                                                                                 gTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGc                                                         <  4:506249/129‑1 (MQ=255)
                                                                                                                      gCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCtt                                <  2:1466565/149‑1 (MQ=255)
                                                                                                                      gCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCtt                                <  1:764646/149‑1 (MQ=255)
                                                                                                                       cAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCttt                               <  4:1164760/149‑1 (MQ=255)
                                                                                                                        aGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTtg                              <  3:108205/149‑1 (MQ=255)
                                                                                                                        aGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTtg                              >  2:1148504/1‑149 (MQ=255)
                                                                                                                        aGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTtg                              <  3:1124084/149‑1 (MQ=255)
                                                                                                                        aGTAACTGACGCTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCTATTGCACCACACTGCGACGGTAACAAGGGTTAGCTGATCACCGCCTGCGGCGCACTCGTtg                              >  2:661958/1‑149 (MQ=255)
                                                                                                                         gTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTtgc                             >  3:336428/1‑149 (MQ=255)
                                                                                                                              tGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCt                        >  1:455968/1‑149 (MQ=255)
                                                                                                                                   gTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCg                   >  2:384238/1‑149 (MQ=255)
                                                                                                                                        gACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGAt                                                     <  4:1096232/110‑1 (MQ=255)
                                                                                                                                         aCATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTc             <  1:1148504/149‑1 (MQ=255)
                                                                                                                                           aTGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATAAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTa           <  1:843239/149‑1 (MQ=255)
                                                                                                                                            aacgggTGAGCATGCCCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTAt          >  3:884834/4‑149 (MQ=255)
                                                                                                                                            aacgggTGAGCATGCCCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGATGCCTGGCCGATGGTCTAt          >  4:720464/4‑149 (MQ=255)
                                                                                                                                                  tGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCt    >  3:1426717/1‑149 (MQ=255)
                                                                                                                                                  tGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCt    >  1:1389069/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGC  >  NC_000913/772990‑773286

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: