Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,654,553 C→G 7.5% intergenic (+26/‑130) slp → / → dctR starvation lipoprotein/putative DNA‑binding transcriptional regulator DctR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,654,5530CG 7.5% 310.4 / 12.0 93intergenic (+26/‑130)slp/dctRstarvation lipoprotein/putative DNA‑binding transcriptional regulator DctR
Reads supporting (aligned to +/- strand):  ref base C (43/42);  new base G (3/4);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GTGGTTAATACCACCTATAACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAG  >  NC_000913/3654405‑3654702
                                                                                                                                                    |                                                                                                                                                     
gTGGTTAATACCACCTATAACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTtct                                                                                                                                                      <  1:148681/150‑1 (MQ=255)
  ggTTAATACCCCCTATAACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAAAGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCa                                                                                                                                                    >  2:1467146/1‑150 (MQ=255)
  ggTTAATACCACCTATAACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCa                                                                                                                                                    >  2:1187490/1‑150 (MQ=255)
      aaTACCACCTATAACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCAt                                                                                                                                                >  2:383756/1‑150 (MQ=255)
       aTACCACCTATAACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATc                                                                                                                                               >  2:1499697/1‑150 (MQ=255)
         accaccTATAACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTg                                                                                                                                             <  1:1457972/150‑1 (MQ=255)
          ccaccTATAACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGc                                                                                                                                            >  1:853457/1‑150 (MQ=255)
           caccTATAACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGTGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTGAGAAAGCATCTGca                                                                                                                                           <  1:504703/150‑1 (MQ=255)
           caccTATAACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGca                                                                                                                                           <  1:1261069/150‑1 (MQ=255)
           caccTATAACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTGAGAAAGCATCTGca                                                                                                                                           <  1:1619712/150‑1 (MQ=255)
            accTATAACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGcat                                                                                                                                          >  2:625435/1‑150 (MQ=255)
             ccTATAAGCTGTGGGATTAGGGCTATGGTGGATTGTGGCCGGAACCGGGCTGGGGGTAGCCTTACTAGACGAATGCGGTGAGTCAGGTAACAGGTGAGCTTGTCAAATAACACCACGGAAAGATGCAGATGGTGAGAAAGCATCTGcatc                                                                                                                                         <  2:312597/150‑1 (MQ=255)
             ccTATAACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGcatc                                                                                                                                         <  2:865324/150‑1 (MQ=255)
             ccTATAACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTGAGAAAGCATCTGcatc                                                                                                                                         <  1:1449483/150‑1 (MQ=255)
               tataACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGcatcat                                                                                                                                       >  2:1019347/1‑150 (MQ=255)
                  aaCCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCa                                                                                                                                    >  1:791817/1‑150 (MQ=255)
                    ccTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCAtt                                                                                                                                  >  2:1124696/1‑150 (MQ=255)
                       gtgGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACa                                                                                                                               >  1:149677/1‑150 (MQ=255)
                          ggATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCa                                                                                                                            <  2:1265899/150‑1 (MQ=255)
                           gATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCaa                                                                                                                           <  2:1127649/150‑1 (MQ=255)
                            aTTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCaaa                                                                                                                          <  1:1531211/150‑1 (MQ=255)
                              tACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGTGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAAtt                                                                                                                        >  2:811454/1‑150 (MQ=255)
                                cGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTaa                                                                                                                      >  2:1614218/1‑150 (MQ=255)
                                       ggTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTTAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAg                                                                                                               <  2:336212/150‑1 (MQ=255)
                                       ggTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAg                                                                                                               >  2:261604/1‑150 (MQ=255)
                                               cTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTaa                                                                                                       >  1:1106975/1‑150 (MQ=255)
                                                       aaCCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTATCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACtattatt                                                                                               >  2:812763/1‑150 (MQ=255)
                                                        aCCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACtattatta                                                                                              <  2:856873/150‑1 (MQ=255)
                                                        aCCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACtattatta                                                                                              <  2:853457/150‑1 (MQ=255)
                                                         ccGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACtattattat                                                                                             <  1:1318281/150‑1 (MQ=255)
                                                         ccGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACtattattat                                                                                             >  1:180169/1‑150 (MQ=255)
                                                           gggCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACtattattatta                                                                                           <  1:353281/150‑1 (MQ=255)
                                                           gggCTGGGGTGCGCCTTACTACACCAATGCGGGGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACtattattatta                                                                                           <  2:970628/150‑1 (MQ=255)
                                                            ggCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACtattattattat                                                                                          >  1:1331678/1‑150 (MQ=255)
                                                              cTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAGTAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATaa                                                                                        >  2:1295513/1‑150 (MQ=255)
                                                                ggggTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGc                                                                                      <  2:1477143/150‑1 (MQ=255)
                                                                  ggTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGccc                                                                                    >  1:891896/1‑150 (MQ=255)
                                                                   gTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCcct                                                                                   >  1:886686/1‑150 (MQ=255)
                                                                       gcCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCcctgtcc                                                                               <  2:1366199/150‑1 (MQ=255)
                                                                       gcCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCcctgtcc                                                                               >  1:1420259/1‑150 (MQ=255)
                                                                        ccTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCcctgtcct                                                                              <  2:180169/150‑1 (MQ=255)
                                                                         cTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAATTAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCcctgtcctg                                                                             <  2:431121/150‑1 (MQ=255)
                                                                           tactacACCAATGTGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGtt                                                                           <  1:811454/150‑1 (MQ=255)
                                                                           tactacACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAGTAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGtt                                                                           <  1:1295513/150‑1 (MQ=255)
                                                                            actacACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTa                                                                          <  2:835431/150‑1 (MQ=255)
                                                                            actacACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTa                                                                          >  1:20496/1‑150 (MQ=255)
                                                                               acacCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAAtt                                                                       <  2:345676/150‑1 (MQ=255)
                                                                               acacCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCACTGTACTGTTAAtt                                                                       >  2:1559247/1‑150 (MQ=255)
                                                                               acacCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTAAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATAAAATTAATACAAAGTAAGCTAACTATTATTATTATAAGCCATGTCATGTTAAtt                                                                       >  2:713907/1‑150 (MQ=255)
                                                                                cacCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTa                                                                      >  1:1133824/1‑150 (MQ=255)
                                                                                aacCAATGCGGTGAGTCAGGTAACACCTTAGCTGTTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTa                                                                      <  1:1596391/149‑1 (MQ=255)
                                                                                      tGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTg                                                                <  2:149677/150‑1 (MQ=255)
                                                                                       gCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTgg                                                               <  2:1514170/150‑1 (MQ=255)
                                                                                          gTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCaa                                                            >  2:198686/1‑150 (MQ=255)
                                                                                            gAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAAc                                                          <  1:1499697/150‑1 (MQ=255)
                                                                                            gAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAAc                                                          >  1:1149332/1‑150 (MQ=255)
                                                                                              gTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTg                                                        >  1:230414/1‑150 (MQ=255)
                                                                                              gTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTAAATTACCTTTGGCAAACTg                                                        >  2:272207/1‑150 (MQ=255)
                                                                                                cAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGAt                                                      <  2:791817/150‑1 (MQ=255)
                                                                                                cAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGAt                                                      >  1:727050/1‑150 (MQ=255)
                                                                                                       cacCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAAt                                                                        >  2:46142/1‑125 (MQ=255)
                                                                                                       cacCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAAt                                                                        <  1:46142/125‑1 (MQ=255)
                                                                                                         ccTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGtt                                             >  1:290345/1‑150 (MQ=255)
                                                                                                          cTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTa                                            >  1:1433907/1‑150 (MQ=255)
                                                                                                                tGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTAAGCATATGCATCATGCATTACATCAAATTAATACACATTAATCTAACTATTATTATTATAAGACCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGtta                                      <  1:1559247/150‑3 (MQ=255)
                                                                                                                 ggTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCg                                     <  2:842650/150‑1 (MQ=255)
                                                                                                                 ggTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCg                                     >  1:1444873/1‑150 (MQ=255)
                                                                                                                   tCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCa                                   >  1:924644/1‑150 (MQ=255)
                                                                                                                     aaaTAACACCACCTAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCAcc                                 <  1:625435/150‑1 (MQ=255)
                                                                                                                     aaaTAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCGTTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAGTTAATGTACccacca                                >  1:219932/1‑150 (MQ=255)
                                                                                                                      aaTAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCa                                <  1:1321062/150‑1 (MQ=255)
                                                                                                                       aTAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAg                               <  1:1124696/150‑1 (MQ=255)
                                                                                                                          acacCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAg                            <  1:55461/150‑1 (MQ=255)
                                                                                                                           caccacCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGt                           >  2:1361152/1‑150 (MQ=255)
                                                                                                                           caccacCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGt                           >  2:1311637/1‑150 (MQ=255)
                                                                                                                               accGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTTATTATAAAGTTAATGTCCGCACCAGGAGTCGGt                       <  1:639951/150‑1 (MQ=255)
                                                                                                                               accGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGt                       <  1:1019347/150‑1 (MQ=255)
                                                                                                                               accGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGt                       <  1:695022/150‑1 (MQ=255)
                                                                                                                               accGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGt                       <  1:294028/150‑1 (MQ=255)
                                                                                                                                ccGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGtt                      >  2:894571/1‑150 (MQ=255)
                                                                                                                                  gcatgatgcagatgcTGAGAAAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTAt                    >  2:1335409/5‑150 (MQ=255)
                                                                                                                                  gAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTAt                    >  2:1391117/1‑150 (MQ=255)
                                                                                                                                  gAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTAt                    >  1:848025/1‑150 (MQ=255)
                                                                                                                                       atgcagatgcTGAGAAAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGttta               >  2:189837/1‑149 (MQ=255)
                                                                                                                                       atgcagatgcTGAGAAAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGttta               >  1:1433428/1‑149 (MQ=255)
                                                                                                                                        tgcagatgcTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCt              >  1:955709/1‑150 (MQ=255)
                                                                                                                                         gcagatgcTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCtt             <  2:963027/150‑1 (MQ=255)
                                                                                                                                          cagatgcTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAAt                                          >  1:852040/1‑120 (MQ=255)
                                                                                                                                          cagatgcTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAAt                                          <  2:852040/120‑1 (MQ=255)
                                                                                                                                           agatgcTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTtat           <  2:955709/150‑1 (MQ=255)
                                                                                                                                            gatgcTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTtata          >  1:1169559/1‑150 (MQ=255)
                                                                                                                                            gatgcTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTtata          <  2:568997/150‑1 (MQ=255)
                                                                                                                                             atgcTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATaa         <  2:1433907/150‑1 (MQ=255)
                                                                                                                                                  ttCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTAcc    <  1:261604/150‑1 (MQ=255)
                                                                                                                                                    ctcAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTAATATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGAGTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAg  <  2:1438353/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGGTTAATACCACCTATAACCTGTGGGATTACGGCTATGGTGCATTCTGGCCGGAACCGGGCTGGGGTGCGCCTTACTACACCAATGCGGTGAGTCAGGTAACACCTGAGCTGGTCAAATAACACCACCGAAAGATGCAGATGCTTTCTCAGCATCTGCATCATGCATTACATCAAATTAATACACAGTAAGCTAACTATTATTATTATAAGCCCTGTCCTGTTAATTACCTTTGGCAAACTGATTATAAAGTTAATGTCCGCACCAGGAGTCGGTTATGTTTCTTATAATTACCAG  >  NC_000913/3654405‑3654702

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