Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 448,965 G→A 6.0% H559H (CAC→CAT) ‡ 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,9650GA 6.0% 369.5 / 8.1 117H559H (CAC→CAT) ‡cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base G (46/63);  new base A (4/3);  total (51/66)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

Only 100 of 121 total aligned reads displayed.

TGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAG  >  NC_000913/448816‑449114
                                                                                                                                                     |                                                                                                                                                     
tGCGAGATGATTTCGAAAATGATGATGGTGTAGAAAGCTGCTATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTGGGCTTTCTCTTTCATTTCCCAGAATGCATAACGTTCg                                                                                                                                                       <  1:103178/150‑1 (MQ=255)
tGCCAGATCATGGCGAAACCGATGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCGGTTCTCTTTCATTTCCCAGAATGCATCACGTTCg                                                                                                                                                       >  2:620859/1‑150 (MQ=255)
 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:50492/150‑1 (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)
              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)
                   cGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACgg                                                                                                                                    <  1:781809/150‑1 (MQ=255)
                   cGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACgg                                                                                                                                    <  1:441243/150‑1 (MQ=255)
                    gAAGATGGTGGAGTAAGCGGCTATGTCGATACCTGCACCGCTGTTGTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGc                                                                                                                                   <  2:350219/150‑1 (MQ=255)
                    gAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGc                                                                                                                                   <  2:745211/150‑1 (MQ=255)
                    gAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGc                                                                                                                                   <  1:729822/150‑1 (MQ=255)
                      aGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCaa                                                                                                                                 >  2:55271/1‑150 (MQ=255)
                       gATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGGTTTTCGGCATATGAATTTCTTCATAGTGGTCA‑GCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAg                                                                                                                               <  1:841305/150‑1 (MQ=255)
                        aTGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAg                                                                                                                               <  1:996963/150‑1 (MQ=255)
                        aTGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAg                                                                                                                               <  1:581691/150‑1 (MQ=255)
                             ggAGAACGCTGCAATGACGAACCCTGCACCGCTGTTTTTCGGCATATGAATTGCTTCATAGAGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTtata                                                                                                                          <  1:108925/150‑1 (MQ=255)
                               agaAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGa                                                                                                                        >  1:1265268/1‑150 (MQ=255)
                                    gCTGCAATGACGATACCTGCACCGCTGTTTTTCGTCATATGAATTTCTTAATAGTGGTCAGGCTTTTTATACGCTTAGACTTTCTATTTCATTTCACAGAATGCATCACGTTCTTGAACGTGAGGCAATACGGCAAAGTTATAGAACGGa                                                                                                                   >  2:238846/1‑150 (MQ=255)
                                     cTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAg                                                                                                                  >  2:452695/1‑150 (MQ=255)
                                       gCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGc                                                                                                                >  2:542233/1‑150 (MQ=255)
                                       gCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGc                                                                                                                >  1:1651287/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)
                                                   ccTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGtt                                                                                                    >  1:1394750/1‑150 (MQ=255)
                                                   ccTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGtt                                                                                                    >  1:1639433/1‑150 (MQ=255)
                                                    cTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTg                                                                                                   <  2:848258/150‑1 (MQ=255)
                                                      gCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGcc                                                                                                 <  2:1035572/150‑1 (MQ=255)
                                                      gCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCGTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGcc                                                                                                 <  2:669514/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)
                                                            cTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTcc                                                                                           >  2:1312455/1‑150 (MQ=255)
                                                                       cATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGcc                                                                                >  1:1013387/1‑150 (MQ=255)
                                                                       cATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATAACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGCCGGGGAAGAGGTTGCCCACTCCAGCGTACGGcc                                                                                >  2:1333695/1‑150 (MQ=255)
                                                                         tatGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCAc                                                                              >  2:1702011/1‑150 (MQ=255)
                                                                            gAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCAcccc                                                                           <  2:36754/150‑1 (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)
                                                                                         gTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGt                                                              <  2:662563/150‑1 (MQ=255)
                                                                                            ggcAGGCTGTTTATACGCTTCGCCGATCTCTTTCCTTTCCCAGAATGCATCACGTGCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAg                                                           <  1:1064101/148‑1 (MQ=255)
                                                                                            gTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGGTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAg                                                           <  1:666503/150‑1 (MQ=255)
                                                                                              cAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCCGTCAGGt                                                         >  2:1532598/1‑150 (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)
                                                                                                                        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)
                                                                                                                                ttCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCCCCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGGCGCGAATAGAAACGTa                       >  2:1600330/1‑150 (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)
                                                                                                                                       aaTGCATCAGCTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGa                <  2:764118/150‑1 (MQ=255)
                                                                                                                                         tGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAAt              <  1:573652/150‑1 (MQ=255)
                                                                                                                                          gccgcACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATa             >  1:24748/5‑150 (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)
                                                                                                                                                cGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAgg       <  1:737328/150‑1 (MQ=255)
                                                                                                                                                  ttCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGCAACGTACATCTGAATAACGAGGCa     >  1:1305145/1‑150 (MQ=255)
                                                                                                                                                     gTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCagag  >  1:1368321/1‑150 (MQ=255)
                                                                                                                                                     gTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCagag  <  2:1408595/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAG  >  NC_000913/448816‑449114

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