Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 448,966 T→C 6.0% H559R (CAC→CGC) ‡ cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913448,9660TC 6.0% 329.7 / 9.2 117H559R (CAC→CGC) ‡cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base T (46/63);  new base C (4/3);  total (50/67)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

Only 100 of 119 total aligned reads displayed.

GCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAG  >  NC_000913/448817‑449114
                                                                                                                                                     |                                                                                                                                                    
gCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGt                                                                                                                                                      >  2:1465021/1‑150 (MQ=255)
gCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGt                                                                                                                                                      >  2:669256/1‑150 (MQ=255)
 ccAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTg                                                                                                                                                     >  2:425399/1‑150 (MQ=255)
    gATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTATCATTTCCCAGAATGCATCACGTTCGTGAAc                                                                                                                                                  >  1:902097/1‑150 (MQ=255)
      tCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGt                                                                                                                                                <  2:217679/150‑1 (MQ=255)
       cATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTg                                                                                                                                               >  1:662563/1‑150 (MQ=255)
        aTGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGc                                                                                                                                              >  1:148214/1‑150 (MQ=255)
         tGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCg                                                                                                                                             <  2:1567963/150‑1 (MQ=255)
          ggCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCAAGAATGCATCACGTTCGTGAACGTGCgg                                                                                                                                            >  2:1427165/1‑150 (MQ=255)
           gCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGc                                                                                                                                           <  1:1169264/150‑1 (MQ=255)
           gCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGACTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGc                                                                                                                                           <  2:1275061/150‑1 (MQ=255)
           gCGAAACCGAAGATGGGTGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCTTAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTgatgc                                                                                                                                           <  1:654010/150‑5 (MQ=255)
             gAATCCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCAc                                                                                                                                         <  1:1015763/150‑1 (MQ=255)
             gAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTTCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCAc                                                                                                                                         <  2:18623/150‑1 (MQ=255)
             gAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGGTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCAc                                                                                                                                         <  2:1644675/150‑1 (MQ=255)
             gAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCAc                                                                                                                                         >  2:855473/1‑150 (MQ=255)
               aaCCGTAGATGGTGGAGAAAGCTGCAATGTCGATACCTGCTCCGCTGTTTTTCGTCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCTTCACGTTCGTGAACGTGCGGCACTa                                                                                                                                       <  1:43252/150‑1 (MQ=255)
               aaCCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTa                                                                                                                                       >  2:737328/1‑150 (MQ=255)
                 ccGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCTCTACg                                                                                                                                     >  2:474326/1‑150 (MQ=255)
                  cGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACgg                                                                                                                                    <  1:441243/150‑1 (MQ=255)
                  cGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACgg                                                                                                                                    <  1:781809/150‑1 (MQ=255)
                   gAAGATGGTGGAGTAAGCGGCTATGTCGATACCTGCACCGCTGTTGTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGc                                                                                                                                   <  2:350219/150‑1 (MQ=255)
                   gAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGc                                                                                                                                   <  1:729822/150‑1 (MQ=255)
                   gAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGc                                                                                                                                   <  2:745211/150‑1 (MQ=255)
                      gATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGGTTTTCGGCATATGAATTTCTTCATAGTGGTCAGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAg                                                                                                                               <  1:841305/150‑1 (MQ=255)
                       aTGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAg                                                                                                                               <  1:996963/150‑1 (MQ=255)
                       aTGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAg                                                                                                                               >  2:573652/1‑150 (MQ=255)
                            ggAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTGCGCCTTACTGTTTCATTTCCCAGAATGCATCACGTTCGTTACAGGGCGGCCCTACGGCAAAGTtata                                                                                                                          >  1:1073898/1‑150 (MQ=255)
                              agaAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGa                                                                                                                        >  1:1265268/1‑150 (MQ=255)
                                  aGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGACACTACGGCAAAGTTATAGAACgg                                                                                                                    <  2:666323/150‑1 (MQ=255)
                                    cTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAg                                                                                                                  >  2:452695/1‑150 (MQ=255)
                                      gCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGc                                                                                                                >  1:1651287/1‑150 (MQ=255)
                                      gCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGc                                                                                                                >  2:542233/1‑150 (MQ=255)
                                      gCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGc                                                                                                                >  2:526651/1‑150 (MQ=255)
                                         aTGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTgatg                                                                                                                                            >  1:1073991/1‑116 (MQ=255)
                                         aTGACGATACCTGCACCGATTTTTTTCGTCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCGAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCggg                                                                                                             <  2:1546523/150‑1 (MQ=255)
                                               aTACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAagag                                                                                                       <  1:9812/150‑1 (MQ=255)
                                                tACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAgg                                                                                                      <  1:452695/150‑1 (MQ=255)
                                                tACCTGCACCGCTGTTTTTCGGCAAATGAATTTCTTCATAGTGGTCAGGCTGTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCAGCACGTTCGTGAACGCGCGGCACTACGGAAAAGTTATAGAACGGAGGCGGGGAAGCgg                                                                                                      <  2:1464206/150‑1 (MQ=255)
                                                  ccTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGtt                                                                                                    >  1:1394750/1‑150 (MQ=255)
                                                   cTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTg                                                                                                   <  1:600962/150‑1 (MQ=255)
                                                   cTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTg                                                                                                   <  2:848258/150‑1 (MQ=255)
                                                     gCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGcc                                                                                                 <  2:1035572/150‑1 (MQ=255)
                                                     gCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCGTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGcc                                                                                                 <  2:669514/150‑1 (MQ=255)
                                                      cACCGCTGTTTTTCTGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGccc                                                                                                <  1:67646/150‑1 (MQ=255)
                                                      cACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTTTAGAACGGAGGCGGGGAAGAGGTTGccc                                                                                                <  2:1640021/150‑1 (MQ=255)
                                                      cACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGccc                                                                                                >  1:55184/1‑150 (MQ=255)
                                                      cACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGccc                                                                                                <  1:984384/150‑1 (MQ=255)
                                                      cACCGCTGTTTTTCGGCATAGGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTGCGCCTTTCTATTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGccc                                                                                                <  2:310023/150‑1 (MQ=255)
                                                           cTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTcc                                                                                           >  2:1312455/1‑150 (MQ=255)
                                                              tttttCGGCATATGAATATCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGATCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGc                                                                                        <  1:54638/150‑1 (MQ=255)
                                                                      cATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGcc                                                                                >  1:1013387/1‑150 (MQ=255)
                                                                      cATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGcc                                                                                >  1:1200934/1‑150 (MQ=255)
                                                                      cATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATAACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGCCGGGGAAGAGGTTGCCCACTCCAGCGTACGGcc                                                                                >  2:1333695/1‑150 (MQ=255)
                                                                        tatGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCAc                                                                              >  2:1702011/1‑150 (MQ=255)
                                                                             aTTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCAc                                                                         <  2:1552603/150‑1 (MQ=255)
                                                                              tttCTTCATAGTGGGCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACg                                                                        <  1:926370/150‑1 (MQ=255)
                                                                                   tcATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCg                                                                   >  2:260028/1‑150 (MQ=255)
                                                                                   tcATAGTGGTCAGGCTTTTTATACGCTTCGCCTATCTCTTTCATTTCCCAGAATGAATCACGTTCGTCAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCg                                                                   <  1:114916/150‑1 (MQ=255)
                                                                                     aTAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGcc                                                                 <  1:1102467/150‑1 (MQ=255)
                                                                                     aTAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATATCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGcc                                                                 <  1:21236/150‑1 (MQ=255)
                                                                                       aGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAg                                                               <  2:590056/150‑1 (MQ=255)
                                                                                          ggTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCa                                                            <  2:813596/150‑1 (MQ=255)
                                                                                           ggcAGGCTGTTTATACGCTTCGCCGATCTCTTTCCTTTCCCAGAATGCATCACGTGCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAg                                                           <  1:1064101/148‑1 (MQ=255)
                                                                                             cAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGt                                                         >  2:554073/1‑150 (MQ=255)
                                                                                               ggCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGAACGGCCACCCCACGGGTCGCCAGTCAGGTCa                                                       >  2:1006697/1‑150 (MQ=255)
                                                                                                     ttATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTc                                                 >  1:631568/1‑150 (MQ=255)
                                                                                                      tataCGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGATCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGATGAGGTTGACCACTCCAGCGTACGGCCACCCCACGGGTCGACAGTCAGGTCACGGTTCt                                                >  1:1005406/1‑150 (MQ=255)
                                                                                                          cGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtc                                            >  2:863541/1‑150 (MQ=255)
                                                                                                            cTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtcgc                                          <  1:55271/150‑1 (MQ=255)
                                                                                                              tCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTAACGGTTCTggtcgcgg                                        >  2:1264733/1‑150 (MQ=255)
                                                                                                                  cTTTCTCTTTCATTTACCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTggtcgcggtcgc                                    >  2:1449019/1‑150 (MQ=255)
                                                                                                                     tctcTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGaa                                 <  1:1465021/150‑1 (MQ=255)
                                                                                                                      ctctTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGGTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAAt                                <  1:425399/150‑1 (MQ=255)
                                                                                                                       tctTTCATTTCCGAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCAGACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATa                               <  2:825842/150‑1 (MQ=255)
                                                                                                                       tctTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATa                               <  2:1386390/150‑1 (MQ=255)
                                                                                                                        ctTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAg                              <  1:1663381/150‑1 (MQ=255)
                                                                                                                        ctTTCATTTCCCAGAATGCATCACGATCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAg                              <  1:669256/150‑1 (MQ=255)
                                                                                                                             aTTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACg                         <  2:1577623/150‑1 (MQ=255)
                                                                                                                                tCCTAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTaa                      >  2:1579518/1‑149 (MQ=255)
                                                                                                                                 cccAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACa                     >  1:163917/1‑150 (MQ=255)
                                                                                                                                  ccAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGTTTCTATTCGCGACCGCGACCAGAACCGtgacc                    <  1:855473/150‑5 (MQ=38)
                                                                                                                                  ccAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTgg                                              >  1:589464/1‑124 (MQ=255)
                                                                                                                                  ccAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTgg                                              <  2:589464/124‑1 (MQ=255)
                                                                                                                                  ccAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTACCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGTTTCTATTCGCGACCGCGACCAGAACCGtgacc                    <  1:542233/150‑5 (MQ=38)
                                                                                                                                   cAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTGCTGGTCGCGGTCGCGAATAGAAACGTACATc                   <  1:1711723/150‑1 (MQ=255)
                                                                                                                                     gAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTg                 <  1:1730155/150‑1 (MQ=255)
                                                                                                                                     gAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTg                 <  1:827230/150‑1 (MQ=255)
                                                                                                                                        tGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAAt              <  1:573652/150‑1 (MQ=255)
                                                                                                                                         gccgcACGTTCACGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAAGGTACATCTGAATa             >  2:247310/5‑150 (MQ=255)
                                                                                                                                         gCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATa             <  2:1265268/150‑1 (MQ=255)
                                                                                                                                           cgcACGTTCGGGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAAc           <  1:336857/148‑1 (MQ=255)
                                                                                                                                            tCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACg          >  1:882103/1‑150 (MQ=255)
                                                                                                                                            gcACGTTCACGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACg          >  2:632917/2‑150 (MQ=255)
                                                                                                                                            gcACGTTCACGAACGTGATGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCCGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACg          >  2:293659/2‑150 (MQ=255)
                                                                                                                                             cACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATATCGa         >  2:533373/1‑150 (MQ=255)
                                                                                                                                              aCGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATTGAAACGGACATCTGAATAACGAg        >  1:1713849/1‑150 (MQ=255)
                                                                                                                                               cGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAgg       <  1:737328/150‑1 (MQ=255)
                                                                                                                                                    gTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCagag  >  1:1368321/1‑150 (MQ=255)
                                                                                                                                                    gTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCagag  <  2:1408595/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAG  >  NC_000913/448817‑449114

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: