Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,477,151 G→A 5.9% P160L (CCG→CTG)  dsdC ← dsd operon activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,477,1510GA 5.9% 265.3 / 4.0 84P160L (CCG→CTG) dsdCdsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base G (42/37);  new base A (3/2);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

CGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTGGGCGATAGAGGGC  >  NC_000913/2477002‑2477300
                                                                                                                                                     |                                                                                                                                                     
cGTTGCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACg                                                                                                                                                       >  1:785664/1‑150 (MQ=255)
 gTTGCTCCATGCCTGTCTGTCATGGTGCAACGTACTGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGATTTGCGCTGACgg                                                                                                                                                      <  1:901327/150‑1 (MQ=255)
 gTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTTCCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATTCTGGCAGGATTTCTTAATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACgg                                                                                                                                                      >  1:1539625/1‑150 (MQ=255)
 gTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTACTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACgg                                                                                                                                                      <  1:29800/150‑1 (MQ=255)
    gCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                                                                                                                                                   <  1:730369/150‑1 (MQ=255)
    gCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                                                                                                                                                   <  2:1679311/150‑1 (MQ=255)
    gCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                                                                                                                                                   <  2:1672838/150‑1 (MQ=255)
    gCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                                                                                                                                                   >  2:1684452/1‑150 (MQ=255)
     cTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca                                                                                                                                                  >  1:1233846/1‑150 (MQ=255)
     cTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca                                                                                                                                                  >  2:727707/1‑150 (MQ=255)
     cTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca                                                                                                                                                  >  1:1496829/1‑150 (MQ=255)
     cTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATAAGTGAGTTGCGCTGACGGCGca                                                                                                                                                  >  1:110749/1‑150 (MQ=255)
       ccATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATATCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCCGTCAGCGcaac                                                                                                                                                <  1:1409262/150‑3 (MQ=255)
         aTGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGATAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCCGTCAGCGcaactc                                                                                                                                              <  1:1480501/150‑5 (MQ=255)
          tGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc                                                                                                                                             >  2:189951/1‑150 (MQ=255)
          tGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc                                                                                                                                             >  2:1601101/1‑150 (MQ=255)
           gCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCAtcatca                                                                                                                                            <  1:651443/150‑1 (MQ=255)
           gCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCAtcatca                                                                                                                                            <  1:657423/150‑1 (MQ=255)
            ccTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa                                                                                                                                           >  2:1522295/1‑150 (MQ=255)
                tctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTa                                                                                                                                       >  2:1578715/1‑150 (MQ=255)
                     cATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCg                                                                                                                                  <  2:667922/150‑1 (MQ=255)
                     cATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCAGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCg                                                                                                                                  >  2:1082936/1‑150 (MQ=255)
                           gCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAAt                                                                                                                            >  1:1547482/1‑150 (MQ=255)
                                 tACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTTACGGCGCATCATCAAAGTAGATCGCCAAATCGATTc                                                                                                                      <  1:1516628/150‑1 (MQ=255)
                                 tACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATAAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTc                                                                                                                      <  2:110749/150‑1 (MQ=255)
                                  aCAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTcc                                                                                                                     >  2:994375/1‑150 (MQ=255)
                                  aCAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTcc                                                                                                                     <  1:937816/150‑1 (MQ=255)
                                   cAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCg                                                                                                                    >  2:712471/1‑150 (MQ=255)
                                   cAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCg                                                                                                                    >  1:956757/1‑150 (MQ=255)
                                   cAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCg                                                                                                                    >  1:1506722/1‑150 (MQ=255)
                                    aGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATTTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCgg                                                                                                                   <  2:1345231/150‑1 (MQ=255)
                                         cacaGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGt                                                                                                              >  1:1288249/1‑150 (MQ=255)
                                           caGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTg                                                                                                            >  1:230099/1‑150 (MQ=255)
                                                  aTTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTg                                                                                                     <  2:553722/150‑1 (MQ=255)
                                                  aTTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTg                                                                                                     <  1:1322189/150‑1 (MQ=255)
                                                  aTTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATAATTGCAGGATTTCTTCATCCATCAGAAATTGATGAGTGAGTTGCGCTGAAGGCGCAACATCAAAGTAGATCGCAAAATCGATTCCGGCACGTTGCAGAtta                                                                                                     >  1:1001642/1‑149 (MQ=255)
                                                      ccGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGt                                                                                                 <  1:1418940/150‑1 (MQ=255)
                                                           ttGGTTAAAGCATGTCTTTGAGCGTATTCCTGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCa                                                                                            <  1:727707/150‑1 (MQ=255)
                                                            tGGTTAAAGCATGTCTTTTAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCAt                                                                                           <  2:209660/150‑1 (MQ=255)
                                                                tAAAGCAGGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTGCTTCATCCCTCAGAAAGTGATGAGTGAGTGGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTAcc                                                                                       <  2:449554/150‑1 (MQ=255)
                                                                  aaGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAg                                                                                     <  2:1496829/150‑1 (MQ=255)
                                                                      aTGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAg                                                                                 <  1:1093507/150‑1 (MQ=255)
                                                                         tCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCAc                                                                              >  2:1641615/1‑150 (MQ=255)
                                                                          cTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACg                                                                             >  2:1210650/1‑150 (MQ=255)
                                                                             tGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTg                                                                          >  2:1719620/1‑150 (MQ=255)
                                                                              gAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGa                                                                         >  1:670977/1‑150 (MQ=255)
                                                                                 cGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCg                                                                      <  1:646687/150‑1 (MQ=255)
                                                                                  gTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGa                                                                     >  2:1703588/1‑150 (MQ=255)
                                                                                    aTTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGaaa                                                                   <  2:670977/150‑1 (MQ=255)
                                                                                     ttCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAAt                                                                  >  1:1455197/1‑150 (MQ=255)
                                                                                     ttCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCAAGTTGCAGATTGACGTTGTCATTACCAGTGACCACGGTGAGCGAAAt                                                                  >  1:965791/1‑150 (MQ=255)
                                                                                      tCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATa                                                                 <  2:230099/150‑1 (MQ=255)
                                                                                       ccGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAg                                                                >  2:1632553/1‑150 (MQ=255)
                                                                                         gggCTACTTACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGAc                                                              <  1:1036705/150‑1 (MQ=255)
                                                                                          ggCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACg                                                             <  1:1641615/150‑1 (MQ=255)
                                                                                          ggCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACg                                                             <  1:994375/150‑1 (MQ=255)
                                                                                          ggCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACg                                                             <  2:780574/150‑1 (MQ=255)
                                                                                           gCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACgg                                                            <  1:1662872/150‑1 (MQ=255)
                                                                                               cATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATAc                                                        >  2:24852/1‑150 (MQ=255)
                                                                                                 tACTGGCATGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGAATCATCAAAGTAGATCGCCAAATCGATTCCGGCTCGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGGGAGCGAAATAGACGGATACCg                                                      <  1:816146/150‑1 (MQ=255)
                                                                                                  aCTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCgg                                                     >  1:818189/1‑150 (MQ=255)
                                                                                                  aCTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCgg                                                     >  2:1257273/1‑150 (MQ=255)
                                                                                                  aCTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCgg                                                     >  2:32269/1‑150 (MQ=255)
                                                                                                  aCTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCgg                                                     >  2:32268/1‑150 (MQ=255)
                                                                                                    tGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCg                                                   <  1:631790/150‑1 (MQ=255)
                                                                                                         ggATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTaa                                              >  1:759603/1‑150 (MQ=255)
                                                                                                         ggATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGATGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTaa                                              <  1:413281/150‑1 (MQ=255)
                                                                                                              tcttcATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCa                                         >  2:20927/1‑150 (MQ=255)
                                                                                                                  catccatcAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTa                                     <  1:1601101/150‑1 (MQ=255)
                                                                                                                  caaccatcAGAAAGTGATGAGTGAGATCCGATGACGGCCGCTCATCAAAGTAGATCACCAAATCGATTCCGGCACGTTGAAGATTGACGTAGTCATTACCAGTTAGAACGGTGAGCGAAATAGACGGATACCGGCGTGAAAAGTCAccaa                                     >  1:320536/4‑148 (MQ=255)
                                                                                                                    tccatcAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAAt                                   >  2:1402464/1‑150 (MQ=255)
                                                                                                                          aGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGc                             >  1:1098038/1‑150 (MQ=255)
                                                                                                                          aGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGc                             <  1:730399/150‑1 (MQ=255)
                                                                                                                           gAATGTGAGGAGTGAGGTGCGCTGACGGCGCATCATCAACGTAGATCGCCAAATCGATTCCGGAACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAATTAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCa                            <  1:334193/150‑1 (MQ=255)
                                                                                                                             aaGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCacc                          <  2:785664/150‑1 (MQ=255)
                                                                                                                               gTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCaccaa                        >  1:956294/1‑150 (MQ=255)
                                                                                                                                   tgagtgagTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAg                    >  2:713775/1‑150 (MQ=255)
                                                                                                                                     ttgatgatGCGCCGTCAGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCCCCTTATGCGGGCCCCAACCAGCa                  >  2:711162/8‑150 (MQ=255)
                                                                                                                                     ttgatgatGCGCCGTCAGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCa                  >  2:538200/8‑150 (MQ=255)
                                                                                                                                        atgatGCGCCGTCAGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTg               >  2:579705/5‑150 (MQ=255)
                                                                                                                                          gtTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTggg             <  1:513864/150‑1 (MQ=255)
                                                                                                                                          gtTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTggg             >  2:1349027/1‑150 (MQ=255)
                                                                                                                                           ttGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTGGGc            <  1:20927/150‑1 (MQ=255)
                                                                                                                                            tGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTGGGCg           <  1:1210650/150‑1 (MQ=255)
                                                                                                                                                  gACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTTAATTCACCTAATGCGGGCACCAACCAGCATTGGGCGATagag     >  2:529807/1‑150 (MQ=255)
                                                                                                                                                  gACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTGGGCGATagag     >  2:1661379/1‑150 (MQ=255)
                                                                                                                                                    cGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGCCGGATACCGGCGTGTAAAGTCACCAAAAGCGGGCACCAACCAGCATTGGGCGATAGAggg   >  2:300696/1‑150 (MQ=255)
                                                                                                                                                    cGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTGGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTGGGCGATAGAggg   <  1:189951/150‑1 (MQ=255)
                                                                                                                                                     ggCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTGGGCGATAGAGGGc  <  1:24852/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTGGGCGATAGAGGGC  >  NC_000913/2477002‑2477300

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