Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,881,358 T→A 7.4% D151V (GAC→GTC)  dmlR ← transcriptional activator of dmlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,881,3580TA 7.4% 255.7 / 8.4 94D151V (GAC→GTC) dmlRtranscriptional activator of dmlA
Reads supporting (aligned to +/- strand):  ref base T (43/44);  new base A (3/4);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

TGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCT  >  NC_000913/1881211‑1881506
                                                                                                                                                   |                                                                                                                                                    
tgatCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACa                                                                                                                                                    >  1:646689/1‑150 (MQ=255)
 gatCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCAt                                                                                                                                                   <  2:574154/150‑1 (MQ=255)
       gcTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGcccc                                                                                                                                             <  1:1555229/150‑3 (MQ=255)
        ctctTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAAt                                                                                                                                            >  2:1165719/1‑150 (MQ=255)
         tctTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATg                                                                                                                                           >  2:957955/1‑150 (MQ=255)
         tctTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATg                                                                                                                                           <  1:1376537/150‑1 (MQ=255)
         tctTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGccccta                                                                                                                                           <  2:646689/150‑5 (MQ=255)
           tttAATGACCAGGCACGGAAATGCAGAAAGATCGGTTATATTTTTCGGCGCGCCATGTTGAGCAATAATTTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGcccctaat                                                                                                                                         <  1:351639/150‑7 (MQ=255)
            ttAATGCCCAGCCACGGCAATGCAGAAAGATCGGGTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGcccctaatt                                                                                                                                        <  2:858883/150‑8 (MQ=255)
                 gACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCTATAAATTCCGGCGACGCGCAGAGAATACGGTAAATGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATc                                                                                                                                   >  2:1372332/1‑150 (MQ=255)
                 gACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATc                                                                                                                                   >  2:1379151/1‑150 (MQ=255)
                 gACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATc                                                                                                                                   >  1:406682/1‑150 (MQ=255)
                       gCACGGCAATGCAGAAAGAACGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATc                                                                                                                             <  2:1290653/150‑1 (MQ=255)
                          cGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGac                                                                                                                          >  2:1673379/1‑150 (MQ=255)
                           ggCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGTGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGaca                                                                                                                         <  1:1716062/150‑1 (MQ=255)
                            gCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCGGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGacac                                                                                                                        >  2:673120/1‑150 (MQ=255)
                            gCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGacac                                                                                                                        <  2:597030/150‑1 (MQ=255)
                              aaTGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGTCTATGTCATCGCCAATGCGAATATCGAGATCGACACCt                                                                                                                      <  2:1128695/150‑1 (MQ=255)
                                tGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGaa                                                                                                                                       >  2:74685/1‑131 (MQ=255)
                                tGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGaa                                                                                                                                       <  1:74685/131‑1 (MQ=255)
                                    taattaTCGGTTAAATGTATCGGCGCGCCATGTTGCGCAATAAAATCAGTCGACGCGCATAGAATACGGTAATTTGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTg                                                                                                                <  1:819992/145‑1 (MQ=255)
                                    gAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCTATGTCATAGCCAATGCGAATATCGAGATCGACACCTTCATTg                                                                                                                <  2:478868/150‑1 (MQ=255)
                                    gAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTg                                                                                                                <  1:1254530/150‑1 (MQ=255)
                                     aaaGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGa                                                                                                               >  2:1240215/1‑150 (MQ=255)
                                        gATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATATGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCa                                                                                                            <  1:1338640/150‑1 (MQ=255)
                                         aTCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCTCCAATGCGAATATCGAGATCGACACCTTCATTGACCaa                                                                                                           >  1:246241/1‑150 (MQ=255)
                                         aTCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCaa                                                                                                           >  1:1225386/1‑150 (MQ=255)
                                           cGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAAt                                                                                                         >  2:701297/1‑150 (MQ=255)
                                            ggTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTCACCAAATc                                                                                                        >  2:1254248/1‑150 (MQ=255)
                                               tAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGAc                                                                                                     <  2:411507/150‑1 (MQ=255)
                                                aaaTGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGCTCGACACCTTCATTGACCAAATCGAcc                                                                                                    >  1:1339312/1‑150 (MQ=255)
                                                aaaTGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGAcc                                                                                                    >  2:1243249/1‑150 (MQ=255)
                                                aaaTGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGAcc                                                                                                    <  1:1685884/150‑1 (MQ=255)
                                                aaaTGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAAAATCGAGATCGACACCTTCATTGACCAAATCGAcc                                                                                                    >  1:686841/1‑150 (MQ=255)
                                                 aaTGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATATGGTAATTGGTCGCCAGTTTGTGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCa                                                                                                   >  2:934067/1‑150 (MQ=255)
                                                      ttCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCgg                                                                                              <  1:888228/150‑1 (MQ=255)
                                                         ggCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCt                                                                                           >  2:496323/1‑150 (MQ=255)
                                                         ggCGCGCCATGTTGCGCAATAAATTCATGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCTACACCTTCATTGACCAAATCGACCAGTCGGTCt                                                                                           <  2:1385917/150‑1 (MQ=255)
                                                            gcgcCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTcc                                                                                        <  1:1699757/150‑1 (MQ=255)
                                                            gcgcCATGGTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTATGTGGTATTTCATAGCCAAGGCGAATATCGAGATCTACACCTTCATTGACCAAATCGACCAGTCGGTCTTcc                                                                                        <  2:19332/150‑1 (MQ=255)
                                                             cgcCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCt                                                                                           <  2:462572/146‑1 (MQ=255)
                                                             cgcCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCt                                                                                           >  1:462572/1‑146 (MQ=255)
                                                             cgcCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCca                                                                                       >  1:1276866/1‑150 (MQ=255)
                                                               ccATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCcaca                                                                                     >  1:1180674/1‑150 (MQ=255)
                                                                      gcgcAATAAATTCCGGCGACGCGCATAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACTCCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAAc                                                                              <  1:169468/150‑1 (MQ=255)
                                                                       cgcAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACg                                                                             >  2:1447135/1‑150 (MQ=255)
                                                                        gcAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGt                                                                            >  2:957312/1‑150 (MQ=255)
                                                                         cAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTa                                                                           <  2:690836/150‑1 (MQ=255)
                                                                         ataTAAATTCCGGCGACGCGCTGAGAATACGGTAATTGGTCGCCTGTTTGCGGGCGATTAAATTAGGGTCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTa                                                                           <  2:756216/148‑1 (MQ=255)
                                                                           aTAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAAt                                                                         >  1:107637/1‑150 (MQ=255)
                                                                            tAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAAtt                                                                        <  1:1243867/150‑1 (MQ=255)
                                                                             aaaTTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTc                                                                       <  2:687414/150‑1 (MQ=255)
                                                                               aTTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAAAGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCACTCGGTCTTCCACATCAAAACGTAATTCCa                                                                     >  2:885771/1‑150 (MQ=255)
                                                                                 tCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGt                                                                   <  2:1339312/150‑1 (MQ=255)
                                                                                  ccGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACTCCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCAAGtt                                                                  >  2:1656098/1‑150 (MQ=255)
                                                                                      cGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCgg                                                              <  2:362304/150‑1 (MQ=255)
                                                                                         cgcgCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATa                                                           <  1:1673379/150‑1 (MQ=255)
                                                                                          gcgcAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATAc                                                          <  1:504314/150‑1 (MQ=255)
                                                                                            gcAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCTATTAAATTAGGGGCTATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGTTACGc                                                        <  1:775526/150‑1 (MQ=255)
                                                                                             cAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGcc                                                       <  2:406682/150‑1 (MQ=255)
                                                                                               gagaATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCtt                                                     >  1:172691/1‑150 (MQ=255)
                                                                                                   aTACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGcc                                                 >  1:1331016/1‑150 (MQ=255)
                                                                                                     aCGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAg                                               <  1:957312/150‑1 (MQ=255)
                                                                                                       ggTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAgcg                                             <  2:686841/150‑1 (MQ=255)
                                                                                                            ttGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCtaat                                        >  2:771320/1‑150 (MQ=255)
                                                                                                             tGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCtaata                                       <  1:1165719/150‑1 (MQ=255)
                                                                                                               gTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAAc                                     >  1:1186910/1‑150 (MQ=255)
                                                                                                                 cGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGc                                   >  1:1152408/1‑150 (MQ=255)
                                                                                                                 cGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGc                                   <  1:369649/150‑1 (MQ=255)
                                                                                                                 cGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGc                                   >  2:1343164/1‑150 (MQ=255)
                                                                                                                   ccAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCg                                 >  1:359051/1‑150 (MQ=255)
                                                                                                                     aGTTTGCGGGCGATTAAATTAGGGGCTATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCAATGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCATCGATAATAACTCCGGa                               >  1:44944/1‑150 (MQ=255)
                                                                                                                       tttGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGc                             <  2:107637/150‑1 (MQ=255)
                                                                                                                         tGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGccc                           >  1:1490012/1‑149 (MQ=255)
                                                                                                                          gCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGccac                          <  1:957955/150‑1 (MQ=255)
                                                                                                                           cGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGccacc                         >  2:1382895/1‑150 (MQ=255)
                                                                                                                              gCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCaccacc                      >  2:1632419/1‑150 (MQ=255)
                                                                                                                              gCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCaccacc                      >  2:1625742/1‑150 (MQ=255)
                                                                                                                              gCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCaccacc                      <  2:1432107/150‑1 (MQ=255)
                                                                                                                                 tttAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTgt                   <  2:1652074/149‑1 (MQ=255)
                                                                                                                                 aTTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATCCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTgt                   <  1:701297/150‑1 (MQ=255)
                                                                                                                                  ttAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCCCCACCTgtc                  >  1:1411524/1‑150 (MQ=255)
                                                                                                                                     aaTTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTgtcgtc               <  2:1225386/150‑1 (MQ=255)
                                                                                                                                     aaTTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTgtcgtc               <  2:342116/150‑1 (MQ=255)
                                                                                                                                      aTTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGCCCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTcc              >  1:1050669/1‑150 (MQ=255)
                                                                                                                                       gcattggCGATGACATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCa             >  2:338932/6‑150 (MQ=255)
                                                                                                                                        cattggCGATGACATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCAAATAACGCCGGAGCCCCCACCTGTCGTCCaa            >  2:602272/5‑150 (MQ=255)
                                                                                                                                           gggCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAAc         <  1:1632419/150‑1 (MQ=255)
                                                                                                                                           gggCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAAc         <  1:1625742/150‑1 (MQ=255)
                                                                                                                                           gggCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAAc         <  1:1343164/150‑1 (MQ=255)
                                                                                                                                           gggCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTGCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGGCCaaacc         >  1:180652/1‑148 (MQ=255)
                                                                                                                                             gCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCa       >  2:1583879/1‑150 (MQ=255)
                                                                                                                                                 tGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAgc   <  1:1447135/150‑1 (MQ=255)
                                                                                                                                                  gTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAgct  >  2:1244320/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCT  >  NC_000913/1881211‑1881506

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: