Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,114,527 A→C 8.0% intergenic (‑175/+23) pppA ← / ← yghJ prepilin peptidase/putative lipoprotein YghJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,114,5270AC 8.0% 291.4 / 17.1 100intergenic (‑175/+23)pppA/yghJprepilin peptidase/putative lipoprotein YghJ
Reads supporting (aligned to +/- strand):  ref base A (48/44);  new base C (4/4);  total (52/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 = 4.31e-01

GGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAG  >  NC_000913/3114381‑3114670
                                                                                                                                                  |                                                                                                                                               
ggATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAgcg                                                                                                                                              >  1:478543/1‑150 (MQ=255)
  aTGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGcc                                                                                                                                            <  1:648773/150‑1 (MQ=255)
   tGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGTCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCt                                                                                                                                           <  2:450925/150‑1 (MQ=255)
   tGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCt                                                                                                                                           <  1:731764/150‑1 (MQ=255)
    gcggcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCtt                                                                                                                                          >  2:43919/1‑150 (MQ=255)
       gcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                                       <  2:1091619/150‑1 (MQ=255)
        cgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                                      <  1:43919/150‑1 (MQ=255)
          tAAAACGCCTTATCCGTCCTCCGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGGTAACCCGGTGCGCCTTATTTCATGCCGGATAAGGCGCTCGCGCCGCATCCgg                                                                                                                                    <  1:1261973/150‑1 (MQ=255)
           aaaaCGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                                                                                                   <  1:1188671/150‑1 (MQ=255)
            aaaCGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                                  >  2:87085/1‑150 (MQ=255)
            aaaCGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                                  <  1:1627425/150‑1 (MQ=255)
              aCGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATAAGGCGCTCGCGCCGCATCCGGcatg                                                                                                                                <  1:1320999/150‑4 (MQ=255)
              aCGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATAAGGCGCTCGCGCCGCATCCGGcatg                                                                                                                                <  2:687108/150‑4 (MQ=255)
                gCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATAAGGCGCTCGCGCCGCATCCGGcatgaa                                                                                                                              <  1:930328/150‑6 (MQ=255)
                     aTCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCtt                                                                                                                         <  2:1101186/150‑1 (MQ=255)
                       ccGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTAc                                                                                                                       <  1:744195/150‑1 (MQ=255)
                       ccGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTAc                                                                                                                       >  2:78999/1‑150 (MQ=255)
                       ccGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTAc                                                                                                                       <  1:1586821/150‑1 (MQ=255)
                        cGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACt                                                                                                                      <  1:1315907/150‑1 (MQ=255)
                        cGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACt                                                                                                                      >  1:694830/1‑150 (MQ=255)
                        cGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACt                                                                                                                      >  1:60055/1‑150 (MQ=255)
                        cGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACt                                                                                                                      >  2:118778/1‑150 (MQ=255)
                         gTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTc                                                                                                                     <  1:1374935/150‑1 (MQ=255)
                             tACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCa                                                                                                                 >  1:857760/1‑150 (MQ=255)
                             tACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCa                                                                                                                 >  2:1309120/1‑150 (MQ=255)
                              aCGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAg                                                                                                                >  2:865643/1‑150 (MQ=255)
                              aCGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCAGCCTACGGGCTTACTCGGCAg                                                                                                                <  2:1265230/150‑1 (MQ=255)
                                gggTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGAc                                                                                                              <  2:1546547/150‑1 (MQ=255)
                                  gtgtCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACAt                                                                                                            >  1:731912/1‑150 (MQ=255)
                                      cTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTa                                                                                                        >  1:1291497/1‑150 (MQ=255)
                                      cTGCCAGCGCAAATAATCGCGGCTTTCCACGTCTGTAGGCCCGATAAGCAGGCTCATCGGGCAAATTTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTa                                                                                                        >  2:318242/1‑150 (MQ=255)
                                       tGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTAt                                                                                                       >  1:1541166/1‑150 (MQ=255)
                                       tGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTAt                                                                                                       >  2:361830/1‑150 (MQ=255)
                                             cgcAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTc                                                                                                                     >  2:834741/1‑130 (MQ=255)
                                             cgcAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTc                                                                                                                     <  1:834741/130‑1 (MQ=255)
                                             cgcAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCgg                                                                                                 <  1:949779/150‑1 (MQ=255)
                                               cAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTa                                                                                               >  2:1485570/1‑150 (MQ=255)
                                               cAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTa                                                                                               <  2:731912/150‑1 (MQ=255)
                                                  ataatCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCTGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAAcc                                                                                            <  1:922484/150‑1 (MQ=255)
                                                       cgcgGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGAtt                                                                                       <  1:1170635/150‑1 (MQ=255)
                                                        gcgGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTa                                                                                      >  1:1233856/1‑150 (MQ=255)
                                                           gCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATg                                                                                   <  1:799144/150‑1 (MQ=255)
                                                                ccccGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTc                                                                              >  2:1204481/1‑150 (MQ=255)
                                                                 cccGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTcc                                                                             >  2:1668446/1‑150 (MQ=255)
                                                                 cccGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTcc                                                                             >  1:1587627/1‑150 (MQ=255)
                                                                 cccGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTcc                                                                             >  1:75526/1‑150 (MQ=255)
                                                                 cccGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTcc                                                                             <  1:1668446/150‑1 (MQ=255)
                                                                           ggCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGtt                                                                   >  1:565296/1‑150 (MQ=255)
                                                                           ggCCCGATAAGCAGGCGCATCGGACAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGtt                                                                   >  2:428904/1‑150 (MQ=255)
                                                                               cGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCgg                                                               <  1:1624110/150‑1 (MQ=255)
                                                                                 aTAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGtt                                                             <  2:1415938/150‑1 (MQ=255)
                                                                                  tAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGttt                                                            <  1:78999/150‑1 (MQ=255)
                                                                                    aGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCg                                                          <  1:1204481/150‑1 (MQ=255)
                                                                                     gCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATg                                                                                   >  1:1211719/1‑124 (MQ=255)
                                                                                     gCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATg                                                                                   <  2:1211719/124‑1 (MQ=255)
                                                                                     gCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCgg                                                         >  1:894909/1‑150 (MQ=255)
                                                                                     gCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCgg                                                         >  1:1533697/1‑150 (MQ=255)
                                                                                     gCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCgg                                                         <  2:1233856/150‑1 (MQ=255)
                                                                                      cAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGc                                                        <  2:1587627/150‑1 (MQ=255)
                                                                                      cAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGc                                                        <  2:618323/150‑1 (MQ=255)
                                                                                      cAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGc                                                        <  1:43764/150‑1 (MQ=255)
                                                                                            cATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCg                                                  <  1:1223575/150‑1 (MQ=255)
                                                                                              tCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAg                                                >  1:85267/1‑150 (MQ=255)
                                                                                                gggCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCg                                                          >  1:1098976/1‑138 (MQ=255)
                                                                                                gggCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCg                                                          <  2:1098976/138‑1 (MQ=255)
                                                                                                 ggCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGa                                             <  1:1184075/150‑1 (MQ=255)
                                                                                                          gtTAACCCGGTGCGCCTTTTTTCATGCCGGATGCCGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCg                                                  <  2:644177/136‑1 (MQ=255)
                                                                                                          gtTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCgt                                    >  2:1651961/1‑150 (MQ=255)
                                                                                                          gtTAACCCGGTGCGCCTTATTTCATGCCGGATGCCGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCg                                                  >  1:644177/1‑136 (MQ=255)
                                                                                                           ttAACCCGGTGCGCCTTATTTCATGTCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCgtg                                   <  1:1553584/150‑1 (MQ=255)
                                                                                                           ttAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCgtg                                   >  2:1197923/1‑150 (MQ=255)
                                                                                                           ttAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCgtg                                   <  2:1533697/150‑1 (MQ=255)
                                                                                                             aaCCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGtgtt                                 >  2:337340/1‑150 (MQ=255)
                                                                                                              aCCGGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGAGATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGtgttg                                <  1:743782/150‑1 (MQ=255)
                                                                                                              aCCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCTAGCGtgttg                                <  2:60055/150‑1 (MQ=255)
                                                                                                               cccGGTGCGCCTTATTTCATGCCGGATGCGGTGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGtgttgt                               >  2:840196/1‑150 (MQ=255)
                                                                                                                 cGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTaa                             >  1:912407/1‑150 (MQ=255)
                                                                                                                   gTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCTAGCGTTTTGTAAGc                           <  2:478543/150‑1 (MQ=255)
                                                                                                                      cgcCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGa                        >  1:1304407/1‑150 (MQ=255)
                                                                                                                       gcCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGGGACGCGAGCGTGTTGTAAGCCGAc                       >  1:109052/1‑150 (MQ=255)
                                                                                                                        ccTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACt                      >  2:863084/1‑150 (MQ=255)
                                                                                                                        ccTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACt                      >  2:429241/1‑150 (MQ=255)
                                                                                                                         cTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTg                     <  2:1669284/150‑1 (MQ=255)
                                                                                                                           tATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGc                   <  2:694830/150‑1 (MQ=255)
                                                                                                                           tATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGc                   <  2:75526/150‑1 (MQ=255)
                                                                                                                            cccgtaggCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGGTCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCt                  >  2:413843/8‑150 (MQ=38)
                                                                                                                            aTTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCt                  >  1:601553/1‑150 (MQ=255)
                                                                                                                            aTTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCt                  >  1:135292/1‑150 (MQ=255)
                                                                                                                             ccgtaggCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTc                 >  1:633753/7‑150 (MQ=38)
                                                                                                                              ttCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTca                <  1:428904/150‑1 (MQ=255)
                                                                                                                               tCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTcac               >  2:1669285/1‑150 (MQ=255)
                                                                                                                                cATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTcaca              >  1:157895/1‑150 (MQ=255)
                                                                                                                                   gCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCg           >  1:1287042/1‑150 (MQ=255)
                                                                                                                                   gCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCg           >  1:258752/1‑150 (MQ=255)
                                                                                                                                    ccGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGc          >  2:111322/1‑150 (MQ=255)
                                                                                                                                      ggataaggcgCTCGCGCCGCATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGccc        >  2:565338/7‑150 (MQ=255)
                                                                                                                                      ggATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGccc        >  2:1585476/1‑150 (MQ=255)
                                                                                                                                        aTGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTc      <  2:1291497/150‑1 (MQ=255)
                                                                                                                                          gcggcgCGAGCGCCTTATCCGGCGTACGGGCTTACTCGGCAGAGATCTTATGCTCGGTAACCTGAGTAATGGTTTCCGGTCGCTGTTCCGGTTTCGGCAGATCGAGTGACGCGCGCGTGTTGTAAGCCGCCTGGCTCACACCGCCCTCGa    >  2:1306036/1‑150 (MQ=255)
                                                                                                                                            ggcgCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAg  >  1:24304/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAG  >  NC_000913/3114381‑3114670

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