Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,405,933 C→G 7.7% A237P (GCA→CCA)  lrhA ← DNA‑binding transcriptional dual regulator LrhA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,405,9330CG 7.7% 297.2 / 10.7 89A237P (GCA→CCA) lrhADNA‑binding transcriptional dual regulator LrhA
Reads supporting (aligned to +/- strand):  ref base C (47/35);  new base G (5/2);  total (52/37)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGG  >  NC_000913/2405785‑2406081
                                                                                                                                                    |                                                                                                                                                    
ctctCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGcc                                                                                                                                                     >  2:701825/1‑150 (MQ=255)
ctctCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGcc                                                                                                                                                     >  2:695159/1‑150 (MQ=255)
   tcCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTc                                                                                                                                                  <  2:475306/150‑1 (MQ=255)
    ccATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCa                                                                                                                                                 >  2:563836/1‑150 (MQ=255)
           cTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAg                                                                                                                                          <  1:1110160/150‑1 (MQ=255)
             gATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGAc                                                                                                                                        <  1:1860917/150‑1 (MQ=255)
                aGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCgg                                                                                                                                     <  1:1277116/150‑1 (MQ=255)
                 gATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGc                                                                                                                                    >  2:279478/1‑150 (MQ=255)
                 gATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAACAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGc                                                                                                                                    >  1:1815865/1‑150 (MQ=255)
                  aTCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCt                                                                                                                                   >  1:110762/1‑150 (MQ=255)
                  aTCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCt                                                                                                                                   <  1:1690946/150‑1 (MQ=255)
                   tCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGTCTTGGCGTGACGGCAAGGCCggttg                                                                                                                                  <  1:206272/150‑4 (MQ=255)
                   tCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCtt                                                                                                                                  <  1:1027604/150‑1 (MQ=255)
                    cACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCttt                                                                                                                                 >  1:224886/1‑150 (MQ=255)
                     aCCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTc                                                                                                                                >  2:123235/1‑150 (MQ=255)
                     aCCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTc                                                                                                                                >  2:721167/1‑150 (MQ=255)
                          ggCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGc                                                                                                                           <  2:1768787/150‑1 (MQ=255)
                          ggCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGc                                                                                                                           >  1:241275/1‑150 (MQ=255)
                          ggCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGc                                                                                                                           <  1:1574197/150‑1 (MQ=255)
                                cTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGAGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACg                                                                                                                     >  1:1094464/1‑150 (MQ=255)
                                  cATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGaa                                                                                                                   >  1:525660/1‑150 (MQ=255)
                                   attattACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAAc                                                                                                                  >  1:24313/1‑150 (MQ=255)
                                      attaCTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGc                                                                                                               <  1:1878963/150‑1 (MQ=255)
                                       ttaCTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCACCCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGcc                                                                                                              <  2:1188273/150‑1 (MQ=255)
                                            ggACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAg                                                                                                         >  2:1660566/1‑150 (MQ=255)
                                              aCGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCg                                                                                                       <  1:1718853/150‑1 (MQ=255)
                                                     aTAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGc                                                                                                <  2:960996/150‑1 (MQ=255)
                                                      tAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCg                                                                                               >  1:1508438/1‑150 (MQ=255)
                                                        gCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGTCTTGGCGTGACGGCAAGGCCggttg                                                                                                                                  >  2:1106466/1‑110 (MQ=38)
                                                        gCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGTCTTGGCGTGACGGCAAGGCCggttg                                                                                                                                  <  1:1106466/113‑4 (MQ=38)
                                                         cacaGCAAATATTCGGTATCAGGTAGCGGCTGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTAACGCCAAGACCGGCTTTCACTGCCTCACGAAAGGCCGGATGCGTCGAGGCGACa                                                                                            >  2:1848190/1‑150 (MQ=255)
                                                           caGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATa                                                                                          >  2:544048/1‑150 (MQ=255)
                                                            aGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATaa                                                                                         >  1:1953981/1‑150 (MQ=255)
                                                             gCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAg                                                                                        >  1:1270470/1‑150 (MQ=255)
                                                                   atTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTc                                                                                  <  1:1198419/150‑1 (MQ=255)
                                                                   atTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTc                                                                                  >  2:1062773/1‑150 (MQ=255)
                                                                   atTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTc                                                                                  >  1:1330475/1‑150 (MQ=255)
                                                                   atTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTATCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTc                                                                                  >  2:695277/1‑150 (MQ=255)
                                                                    ttCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCATGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCg                                                                                 >  1:1026193/1‑150 (MQ=255)
                                                                    ttCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCg                                                                                 >  1:1418307/1‑150 (MQ=255)
                                                                      cGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGcc                                                                               >  1:907005/1‑150 (MQ=255)
                                                                       ggTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCCGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCa                                                                              >  2:739304/1‑150 (MQ=255)
                                                                              ggTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAata                                                                       >  1:1563680/1‑150 (MQ=255)
                                                                                    ggcggTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCt                                                                 >  1:1467025/1‑150 (MQ=255)
                                                                                     gcggTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCtt                                                                <  2:491269/150‑1 (MQ=255)
                                                                                      cggTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCttt                                                               <  1:1961195/150‑1 (MQ=255)
                                                                                       ggTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTtta                                                              <  1:563836/150‑1 (MQ=255)
                                                                                       ggTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTtta                                                              >  2:1701902/1‑150 (MQ=255)
                                                                                       ggTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTtta                                                              <  2:24313/150‑1 (MQ=255)
                                                                                        gTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTttat                                                             <  1:1701902/150‑1 (MQ=255)
                                                                                          aaGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGACATGGAATATCTGCTttatta                                                           >  2:1812303/1‑150 (MQ=255)
                                                                                              cATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCg                                                       <  1:23588/150‑1 (MQ=255)
                                                                                                  tACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCg                                                    <  2:1518535/149‑1 (MQ=255)
                                                                                                  tACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCg                                                    >  1:1518535/1‑149 (MQ=255)
                                                                                                  tACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGc                                                   >  1:1741401/1‑150 (MQ=255)
                                                                                                   aCTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGcc                                                  <  1:1062773/150‑1 (MQ=255)
                                                                                                   aCTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGcc                                                  <  2:1815865/150‑1 (MQ=255)
                                                                                                    cTCCACTCAGCACACGCAGGTCCGTGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCa                                                 <  2:1094464/150‑1 (MQ=255)
                                                                                                    cTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCATGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCa                                                 <  2:1026193/150‑1 (MQ=255)
                                                                                                      ccACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGc                                               >  1:153246/1‑150 (MQ=255)
                                                                                                      ccACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTCTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGc                                               <  1:695277/150‑1 (MQ=255)
                                                                                                           cAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCat                                          >  1:1632984/1‑150 (MQ=255)
                                                                                                             gCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCatat                                        >  1:1071715/1‑150 (MQ=255)
                                                                                                                 acGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATcgcg                                    <  2:110762/150‑1 (MQ=255)
                                                                                                                 acGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATcgcg                                    >  2:1109076/1‑150 (MQ=255)
                                                                                                                  cGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGa                                   >  2:764760/1‑150 (MQ=255)
                                                                                                                     aGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                                <  2:1270470/150‑1 (MQ=255)
                                                                                                                      ggTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACg                               >  1:701041/1‑150 (MQ=255)
                                                                                                                        tCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACggg                             >  1:352630/1‑150 (MQ=255)
                                                                                                                            ggCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAg                         >  1:1455804/1‑150 (MQ=255)
                                                                                                                              cTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGa                       <  1:544048/150‑1 (MQ=255)
                                                                                                                                catcatCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGatc                     >  2:1653531/1‑150 (MQ=255)
                                                                                                                                catcatCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGatc                     >  2:826742/1‑150 (MQ=255)
                                                                                                                                 atcatcTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGatca                    >  1:838993/1‑150 (MQ=255)
                                                                                                                                 atcatcTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGatca                    >  1:1461570/1‑150 (MQ=255)
                                                                                                                                 atcatcTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGatca                    <  2:330876/150‑1 (MQ=255)
                                                                                                                                  tcatcTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGatcct                   >  2:1700500/1‑148 (MQ=255)
                                                                                                                                    atcTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATcc                 <  1:488499/150‑1 (MQ=255)
                                                                                                                                    atcTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATcc                 >  2:327721/1‑150 (MQ=255)
                                                                                                                                      ctcAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGatca                    <  2:801715/145‑1 (MQ=255)
                                                                                                                                      ctcAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGatca                    >  1:801715/1‑145 (MQ=255)
                                                                                                                                          aCCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAAt           <  2:153246/150‑1 (MQ=255)
                                                                                                                                           ccGGTCTTGGCGTGACGGCAAGGCCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATa          >  1:1690602/1‑150 (MQ=255)
                                                                                                                                           ccGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATa          <  1:1835297/150‑1 (MQ=255)
                                                                                                                                            cGGTCTTGGCGTGACGGCAAGGCCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATAc         >  1:1580629/1‑150 (MQ=255)
                                                                                                                                             ggtcTTGGCGTGACGGCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACa        >  1:1533883/4‑150 (MQ=255)
                                                                                                                                             ggtcTTGGCGTGACGGCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACa        >  2:1526108/4‑150 (MQ=255)
                                                                                                                                                cTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAaga     >  2:1077760/1‑150 (MQ=255)
                                                                                                                                                  tGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAgg   >  1:742281/1‑150 (MQ=255)
                                                                                                                                                   gCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAggg  <  2:1632984/150‑1 (MQ=255)
                                                                                                                                                   gCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAggg  <  2:525660/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCTCCATCGCCTGATAGATCACCTGGGCCAGCTCATTATTACTGGACGGATCATAGCACAGCAAATATTCGGTATCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGG  >  NC_000913/2405785‑2406081

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: