Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 458,744 T→A 5.1% intergenic (+44/‑144) clpX → / → lon ATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease/DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913458,7440TA 5.1% 343.1 / 4.9 112intergenic (+44/‑144)clpX/lonATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease/DNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base T (58/48);  new base A (2/4);  total (60/52)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

Only 100 of 116 total aligned reads displayed.

GGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGA  >  NC_000913/458595‑458891
                                                                                                                                                     |                                                                                                                                                   
ggAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTACCATACAATTAGTTAACCAAAAAGGGGGGAtttt                                                                                                                                                     >  2:510219/1‑150 (MQ=255)
 gAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATATGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGCGGGATTTTa                                                                                                                                                    >  2:1011067/1‑150 (MQ=255)
  aaGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAt                                                                                                                                                   >  1:33470/1‑150 (MQ=255)
   aGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATTGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAtc                                                                                                                                                  <  1:180576/150‑1 (MQ=255)
   aGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCACCAGGCATCTGGTGAATAAGTACCCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAtc                                                                                                                                                  >  2:304945/1‑150 (MQ=255)
   aGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAtc                                                                                                                                                  >  2:62958/1‑150 (MQ=255)
   aGACGTAGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGACGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAtc                                                                                                                                                  <  2:24297/150‑1 (MQ=255)
    gACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAtct                                                                                                                                                 <  1:1321347/150‑1 (MQ=255)
         cGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCt                                                                                                                                            >  1:1483685/1‑150 (MQ=255)
           aaaaaGTGGTTATCGACTAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAATCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTACCCAAAAAGGGGGGATTTTATCTCCCCttt                                                                                                                                          <  1:1856090/150‑1 (MQ=255)
           aaaaaGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCttt                                                                                                                                          >  1:54179/1‑150 (MQ=255)
              aaGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCtttttg                                                                                                                                       <  2:1377781/150‑4 (MQ=255)
              aaGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCtttttg                                                                                                                                       <  2:1866718/150‑4 (MQ=255)
               aGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAAtt                                                                                                                                      >  1:814727/1‑150 (MQ=255)
                  ggTTCTCGAAGAGTCGGTAATTGATGGTCAAAGCAAACCTATGCTGATTTATGTCAAGCCGGAAGAGCAACAGGAATCTGGTGAATAATTAACCATACCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAAttttt                                                                                                                                   <  2:324524/150‑1 (MQ=255)
                  ggTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAAttttt                                                                                                                                   >  2:1624551/1‑150 (MQ=255)
                  ggTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAAttttt                                                                                                                                   >  2:47072/1‑150 (MQ=255)
                   tttttCGACGATTCTGTAATTTATTGTCAAAGCAAACGGTTGCTGATTTATGGCTAGCCGGAAGCGCAACAGGCATCTGGTGATTAATTAACCATTCCCATACTATTAGTTAACCAAAAAGGGGTGATGTTATCTCCCCTTTAATTTTTc                                                                                                                                  <  1:1753338/146‑1 (MQ=255)
                    ttATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTcc                                                                                                                                 <  1:544001/150‑1 (MQ=255)
                    ttATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTcc                                                                                                                                 >  1:1790163/1‑150 (MQ=255)
                     tATCGACGAGTCTGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCct                                                                                                                                <  1:297315/150‑1 (MQ=255)
                     tATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACGATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCct                                                                                                                                <  2:1138710/150‑1 (MQ=255)
                     tATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCct                                                                                                                                >  2:1034403/1‑150 (MQ=255)
                     tATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCct                                                                                                                                >  1:1037796/1‑150 (MQ=255)
                      aTCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCctc                                                                                                                               <  1:690143/150‑1 (MQ=255)
                        cgacgaGTCGGTAATTGATTGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTa                                                                                                                             <  1:404200/150‑1 (MQ=255)
                           cgaGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATtc                                                                                                                          <  2:1223692/150‑1 (MQ=255)
                              gTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCg                                                                                                                       >  1:1704095/1‑150 (MQ=255)
                                cGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGc                                                                                                                     >  1:818003/1‑150 (MQ=255)
                                cGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGAGAATAATTAACCATTCCCAAACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTATACTCTATTCTCGGc                                                                                                                     >  1:379757/1‑150 (MQ=255)
                                 ggTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCg                                                                                                                    <  2:1242045/150‑1 (MQ=255)
                                    aaTTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTg                                                                                                                 <  2:648685/150‑1 (MQ=255)
                                    aaTTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTg                                                                                                                 <  2:18568/150‑1 (MQ=255)
                                        gATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAAtg                                                                                                             >  1:177163/1‑150 (MQ=255)
                                            gtaaaaGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTggg                                                                                                         <  1:99722/147‑1 (MQ=255)
                                            gTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTggg                                                                                                         >  1:1560869/1‑150 (MQ=255)
                                             tCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTgggg                                                                                                        >  1:1041684/1‑150 (MQ=255)
                                               aaaGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCATATATTATCGGCTTTTAATGTGGGGGa                                                                                                      >  1:1772983/1‑150 (MQ=255)
                                                  gCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTAATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAAc                                                                                                   >  1:483680/1‑150 (MQ=255)
                                                        cGTTGCTGATTTATGGCAAGCCGGACGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATcccc                                                                                             <  2:1040328/150‑1 (MQ=255)
                                                        cGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATcccc                                                                                             <  2:1727356/150‑1 (MQ=255)
                                                            gCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCatat                                                                                         <  1:1767044/150‑1 (MQ=255)
                                                             cTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCAaata                                                                                        >  2:638092/1‑150 (MQ=255)
                                                                aTTTATGGCGCGCCGGAAGCGCAAAAGGCTTCTGGTGAATAATTAACCATTCCGATACAATTTGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGAGTTGAATGTGGGGGAAACATCCCCATATACTg                                                                                     >  1:1090229/1‑150 (MQ=255)
                                                                aTTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTg                                                                                     >  1:546450/1‑150 (MQ=255)
                                                                aTTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTg                                                                                     >  1:1119165/1‑150 (MQ=255)
                                                                 tttATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGa                                                                                    <  1:1375503/150‑1 (MQ=255)
                                                                   tATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACg                                                                                  <  2:1667551/150‑1 (MQ=255)
                                                                   tATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACg                                                                                  <  2:682036/150‑1 (MQ=255)
                                                                   tATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACg                                                                                  <  1:886195/150‑1 (MQ=255)
                                                                   tATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACg                                                                                  <  1:911243/150‑1 (MQ=255)
                                                                    aTGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGt                                                                                 >  1:1595163/1‑150 (MQ=255)
                                                                       gCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACa                                                                              <  2:385361/150‑1 (MQ=255)
                                                                         aaGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATg                                                                            <  2:884542/150‑1 (MQ=255)
                                                                          aGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGt                                                                           >  2:880036/1‑150 (MQ=255)
                                                                           gCCGGAAGCGCAACAGGCATCTTGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGAACATGtt                                                                          <  2:136019/150‑1 (MQ=255)
                                                                           gCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGtt                                                                          >  2:322368/1‑150 (MQ=255)
                                                                           gCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGtt                                                                          <  1:1434352/150‑1 (MQ=255)
                                                                                 aGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTATTTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATAAACTGACGTACATGTTAATAGa                                                                    >  2:901737/1‑150 (MQ=255)
                                                                                 aGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGa                                                                    <  1:614369/150‑1 (MQ=255)
                                                                                  gcgcATCAGGCATCTGTTGAATAATTTACCAATCCCATAAAATTAGTTAACCAAAAAGGGTGGATTTTATCTCCACTTTAATTTTTCCTCTATTCTCGGAGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGAt                                                                   <  1:476953/150‑1 (MQ=255)
                                                                                  gcgcAACAGGCATCTGGTGAATAATTTACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGAt                                                                   >  2:330205/1‑150 (MQ=255)
                                                                                    gcAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTAGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTAAATGTTAATAGATgg                                                                 >  1:178279/1‑150 (MQ=255)
                                                                                        cAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTg                                                             <  2:33470/150‑1 (MQ=255)
                                                                                                tGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTAAATAGATGGCGTGAAGCACAg                                                     >  1:1895429/1‑150 (MQ=255)
                                                                                                tGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATATTAATATATGGCTTGAAGCACAg                                                     >  2:1044639/1‑150 (MQ=255)
                                                                                                 ggTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGt                                                    <  1:1103530/150‑1 (MQ=255)
                                                                                                 ggTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCACCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGt                                                    >  1:169181/1‑150 (MQ=255)
                                                                                                         atTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCaa                                            >  2:1434664/1‑149 (MQ=255)
                                                                                                         atTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTAATCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTTACGTACATGTTAATAGATGGCGTGATGCACAGTCgtgtaaa                                            >  1:26060/1‑147 (MQ=255)
                                                                                                              ccATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGTGAAACACCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGAt                                       <  1:1011067/150‑1 (MQ=255)
                                                                                                               cATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGAtt                                      >  2:1551451/1‑150 (MQ=255)
                                                                                                                 ttCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTAc                                    <  2:1851847/150‑1 (MQ=255)
                                                                                                                     cATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTgg                                >  2:568829/1‑150 (MQ=255)
                                                                                                                     cATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTgg                                >  2:1888018/1‑150 (MQ=255)
                                                                                                                      aTACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACTTGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGc                               >  1:1378079/1‑150 (MQ=255)
                                                                                                                          aaTTAGTTTACCAAAAAGGGTGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACTTGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGaa                           <  1:504500/150‑1 (MQ=255)
                                                                                                                           aTTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTCTTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGGGAAGCACAGTCGTGTCATCTGATTACCTGGCGGaaa                          <  2:1700135/150‑1 (MQ=255)
                                                                                                                            ttAGTTAACCAAAAAGGGTGGATTTAATCTACCCTTTAATTTTTCCTCTATTCTCGGCGTTGAAGTTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAAt                         <  1:292250/150‑1 (MQ=255)
                                                                                                                            ttAGTTAACCAAAAAGGGGGGTTTTTATCTCCCCTTTAATTTTTCCTCTAATCTCGGCGTTGTATTGGGGGTAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAAt                         <  1:945454/150‑1 (MQ=255)
                                                                                                                            ttAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAAt                         <  1:1833545/150‑1 (MQ=255)
                                                                                                                              aGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTa                       <  2:718284/150‑1 (MQ=255)
                                                                                                                               gTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTaa                      >  1:544325/1‑150 (MQ=255)
                                                                                                                               gTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCGGGCGGAAATTaa                      >  1:536936/1‑150 (MQ=255)
                                                                                                                                 tAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGTAAACATCCCCATATCCTGACTTACATGTTAATAGATTTCGTGAAGCACAGTCTTGTCATCTGAATACATCGCTGAAAATaaaa                    >  1:159421/1‑149 (MQ=255)
                                                                                                                                 tAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAAc                    <  1:1034403/150‑1 (MQ=255)
                                                                                                                                 tAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAAc                    >  2:1011364/1‑150 (MQ=255)
                                                                                                                                 tAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAAc                    <  1:889743/150‑1 (MQ=255)
                                                                                                                                 aaaaattaaaGGGGAGATAAAATCCCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTTATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAATGTAAAc                    >  2:224376/8‑150 (MQ=255)
                                                                                                                                  aaCCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACt                   >  2:1617196/1‑150 (MQ=255)
                                                                                                                                  aaCCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACt                   >  2:1473963/1‑150 (MQ=255)
                                                                                                                                   aCCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTa                  >  1:456520/1‑150 (MQ=255)
                                                                                                                                   aCCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTa                  >  1:710844/1‑150 (MQ=255)
                                                                                                                                    ccAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTaa                 >  2:621135/1‑150 (MQ=255)
                                                                                                                                     cAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAag                >  2:751757/1‑150 (MQ=255)
                                                                                                                                       aaaaGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAagag              <  1:47072/150‑1 (MQ=255)
                                                                                                                                          aGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATATGATTTCCTGGCGGATATTTAACTAagagaga           >  1:756829/1‑150 (MQ=255)
                                                                                                                                              gggATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCGGGCGGAAATTAAACTAAGAGAGAGctg       >  2:1726833/1‑149 (MQ=255)
                                                                                                                                               ggATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGctct      <  2:1523571/150‑1 (MQ=255)
                                                                                                                                                   tttATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTTATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGa  >  1:1428556/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGA  >  NC_000913/458595‑458891

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