Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 458,738 G→A 8.9% intergenic (+38/‑150) 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,7380GA 8.9% 299.3 / 12.2 103intergenic (+38/‑150)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 G (37/55);  new base A (2/7);  total (40/63)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

Only 100 of 107 total aligned reads displayed.

GTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCT  >  NC_000913/458589‑458887
                                                                                                                                                     |                                                                                                                                                     
gTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGGGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAggggg                                                                                                                                                       <  1:1046000/150‑1 (MQ=255)
gTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAggggg                                                                                                                                                       <  2:745043/150‑1 (MQ=255)
   cATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGAt                                                                                                                                                    >  2:1698542/1‑150 (MQ=255)
    aTGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGAtt                                                                                                                                                   <  1:1714313/150‑1 (MQ=255)
    aTGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGAtt                                                                                                                                                   >  1:326708/1‑150 (MQ=255)
     tGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGAttt                                                                                                                                                  <  1:755860/150‑1 (MQ=255)
       gAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACGGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTa                                                                                                                                                >  2:944211/1‑150 (MQ=255)
           aCGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAATCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGTGGGGATTTTAtctt                                                                                                                                            >  2:1033819/1‑149 (MQ=255)
             gTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTACAATACAATTAGTTAACCAAAAAGGGGGGATTTTATCGccc                                                                                                                                          >  1:282094/1‑150 (MQ=255)
              tCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTcccc                                                                                                                                         <  1:960001/150‑1 (MQ=255)
              tCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTcccc                                                                                                                                         >  2:702875/1‑150 (MQ=255)
               ggAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTTTGGCAAGCCGGAAGCGCAACATGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCt                                                                                                                                        <  2:1144575/149‑1 (MQ=255)
               cGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCATAGCAAACCGTTGCTGATTTATGGCAAGCCGGTAGCTCAACTGTAATCTGGTGAATAATTCACCATTCCGATACAATTAGTTTTCCAAAAAGGGGAGATAAAATCCCCCCt                                                                                                                                        <  1:1164029/150‑1 (MQ=255)
               cGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCACGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTGCCCt                                                                                                                                        >  1:107967/1‑150 (MQ=255)
               cGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCt                                                                                                                                        >  1:1679730/1‑150 (MQ=255)
               cGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCtt                                                                                                                                       <  1:1671129/150‑1 (MQ=255)
                 aaaaaGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAACCCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCttt                                                                                                                                      <  2:344428/150‑1 (MQ=255)
                  aaaaGTGTTTATCGACGAGTCGGTAATTGATGGTCAACGCAAACCGTTGCTGATTAATGGCAAGCCGGAAGCGCAACAGGCATCTTGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTa                                                                                                                                     <  1:285093/150‑1 (MQ=255)
                   aaaGTGGTTATCGTCGAGTAGGATATTGTTGGTCATAGCAAATCGTTGCTGTTTTATGGCAAGCCGGATGCGGTACAGTCATCTGTTGAATAATTAACCTTTACGTTACAATTAGTTAACCAAAAAGGGGAGATTAAATCCCCCCttttt                                                                                                                                    <  1:1112574/150‑3 (MQ=255)
                    aaGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAAt                                                                                                                                   <  2:225674/150‑1 (MQ=255)
                    aaGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAAt                                                                                                                                   <  1:791450/150‑1 (MQ=255)
                    aaGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAAt                                                                                                                                   <  2:424310/150‑1 (MQ=255)
                    aaGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCtttttg                                                                                                                                   <  2:249371/150‑4 (MQ=255)
                    aaGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCtttttg                                                                                                                                   <  1:1386684/150‑4 (MQ=255)
                     aGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGCATTTTATCTCCCCTTTAAtt                                                                                                                                  >  2:1451725/1‑150 (MQ=255)
                      gTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAAttt                                                                                                                                 <  1:925553/150‑1 (MQ=255)
                        ggTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGTCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCtttttggtta                                                                                                                               <  1:605630/150‑8 (MQ=255)
                           tATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTcca                                                                                                                            >  1:1680355/1‑149 (MQ=255)
                           tATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCct                                                                                                                            >  1:665316/1‑150 (MQ=255)
                            aTCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCctc                                                                                                                           >  1:1405282/1‑150 (MQ=255)
                            aTCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCctc                                                                                                                           >  2:100614/1‑150 (MQ=255)
                            aTCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTCctct                                                                                                                          <  2:627417/150‑1 (MQ=255)
                             tCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCctct                                                                                                                          >  2:1265021/1‑150 (MQ=255)
                              cgacgaGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTa                                                                                                                         >  2:552476/1‑150 (MQ=255)
                                 cgaGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATtc                                                                                                                      >  2:917665/1‑150 (MQ=255)
                                     tCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCgg                                                                                                                  <  2:661566/150‑1 (MQ=255)
                                     tCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCgg                                                                                                                  <  2:638009/150‑1 (MQ=255)
                                      cGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGc                                                                                                                 >  1:1316768/1‑150 (MQ=255)
                                      cGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGc                                                                                                                 <  2:953952/150‑1 (MQ=255)
                                       ggTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCg                                                                                                                <  2:326708/150‑1 (MQ=255)
                                       ggTAATTGATGGTCAAAGCAAACAGTAGCTCATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTACTTAACCAAAAACGGGCGATTTTATATCACCTATAATTATTCCTATATTATCAGCg                                                                                                                >  2:765348/1‑150 (MQ=255)
                                        gTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGt                                                                                                               <  1:569657/150‑1 (MQ=255)
                                        gTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGt                                                                                                               <  2:47309/150‑1 (MQ=255)
                                          aaTTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTg                                                                                                             >  2:411206/1‑150 (MQ=255)
                                          aaTTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTg                                                                                                             >  2:404861/1‑150 (MQ=255)
                                                tGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTgg                                                                                                      <  1:668549/150‑1 (MQ=255)
                                                  ggaaaaTCAACCCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATAAAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTggg                                                                                                     <  2:257405/147‑1 (MQ=255)
                                                  gTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTggg                                                                                                     <  2:193087/150‑1 (MQ=255)
                                                   tCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATAAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTgggg                                                                                                    <  2:1070275/150‑1 (MQ=255)
                                                    cAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTggggg                                                                                                   <  2:222913/150‑1 (MQ=255)
                                                      aaGCAAACCGTTGCTGATTTATTGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAATAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGaa                                                                                                 <  1:1502979/150‑1 (MQ=255)
                                                      aaGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGaa                                                                                                 <  1:8013/150‑1 (MQ=255)
                                                      aaGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGaa                                                                                                 >  2:1256590/1‑150 (MQ=255)
                                                       aGCAAACCGTTGCGGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGaaa                                                                                                >  2:1254153/1‑150 (MQ=255)
                                                         cAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACa                                                                                              >  2:805162/1‑150 (MQ=255)
                                                         cAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCACATTAATTTTTCCTCTATTCTCGTCGTTTAATGTGGGGGAAACa                                                                                              >  2:1426657/1‑150 (MQ=255)
                                                          aaaCCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACaa                                                                                             >  1:185346/1‑149 (MQ=255)
                                                              cGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATcccc                                                                                         <  2:665316/150‑1 (MQ=255)
                                                                  gCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCatat                                                                                     >  1:528856/1‑150 (MQ=255)
                                                                       tttATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTTACCATTCCCATACACTTAGTTAACCAAAAAGGGGTGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTATAATGTGGGGGAAACATCCCCATATACTGa                                                                                <  2:102480/150‑1 (MQ=255)
                                                                        ttATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGAc                                                                               <  1:196396/150‑1 (MQ=255)
                                                                         tATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACg                                                                              >  1:963950/1‑150 (MQ=255)
                                                                            ggCAAGCCGGAAGCGCAACTGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTCACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTAc                                                                           <  2:828078/150‑1 (MQ=255)
                                                                            ggCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTAc                                                                           <  2:1136684/150‑1 (MQ=255)
                                                                             gCAAGCCGGAAGCGCAACATGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACa                                                                          <  2:1625275/150‑1 (MQ=255)
                                                                             gCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATTCTGACGTACa                                                                          >  1:318876/1‑150 (MQ=255)
                                                                             gCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACa                                                                          <  1:1254153/150‑1 (MQ=255)
                                                                                    ggAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACACCCCCATATACTGACGTACATGTTAAt                                                                   <  2:1304149/150‑1 (MQ=255)
                                                                                            aaCAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCg                                                           <  2:1316768/150‑1 (MQ=255)
                                                                                              aaGGCATCTTGTGAATAAATAACCATTCCCATACAATTAGTTAACCAAAAATTTGGGATTTTATCTCCCCTTAAATTTTTCCTCTATTCTCGGCGTTGATTGTGGGCGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCTTg                                                         <  1:241394/149‑1 (MQ=255)
                                                                                               gggCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGa                                                        <  1:944211/149‑1 (MQ=255)
                                                                                                 gCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAg                                                      <  2:524377/150‑1 (MQ=255)
                                                                                                 gCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAg                                                      >  1:493958/1‑150 (MQ=255)
                                                                                                  cATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGc                                                     >  1:640812/1‑150 (MQ=255)
                                                                                                   aTCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGca                                                    <  2:459654/150‑1 (MQ=255)
                                                                                                     cTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATATCCCCTTAAATTTTTCATCTATTCTCGGCGTTGAATTTGGGGGAAACATCCCAATATAATGACGAACATGTTAATAGATGGCGTGAAGaaca                                                  >  2:1574507/1‑150 (MQ=255)
                                                                                                       ggTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTTGTGTAAAAATCCCCATATACTGACGTAAATGTTAATAGATGTCGTGAAGCACAGt                                                >  1:888522/1‑150 (MQ=255)
                                                                                                       ggTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGt                                                <  2:1680355/150‑1 (MQ=255)
                                                                                                         tGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCg                                              >  1:770695/1‑150 (MQ=255)
                                                                                                           aataatTAACCATTCCTATACAATTAGTTAACAAATAAGGGGGGATTTTATCTTCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgtg                                            <  1:1033819/150‑1 (MQ=255)
                                                                                                           aataatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGCGAAGCACAGTCgtg                                            <  2:359113/150‑1 (MQ=255)
                                                                                                                 tAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCt                                      >  2:1496213/1‑150 (MQ=255)
                                                                                                                     cATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGAtt                                  <  2:1288658/150‑1 (MQ=255)
                                                                                                                       ttCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTAc                                <  2:1450868/150‑1 (MQ=255)
                                                                                                                        tCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTAcc                               <  1:468752/150‑1 (MQ=255)
                                                                                                                          ccttACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAAGGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTg                             <  1:100614/147‑1 (MQ=255)
                                                                                                                           cATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTgg                            >  2:1251877/1‑150 (MQ=255)
                                                                                                                            aTACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGc                           <  1:411206/150‑1 (MQ=255)
                                                                                                                             tACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTTCCTGGCg                          >  1:72627/1‑150 (MQ=255)
                                                                                                                             tACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCg                          >  1:1226238/1‑150 (MQ=255)
                                                                                                                             tACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCg                          <  1:1256590/150‑1 (MQ=255)
                                                                                                                               cAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGa                        <  2:963950/150‑1 (MQ=255)
                                                                                                                                    aGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTa                   <  1:1363295/150‑1 (MQ=255)
                                                                                                                                          ccAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTaa             <  2:1405282/150‑1 (MQ=255)
                                                                                                                                           cAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACaaat            >  2:1097118/1‑149 (MQ=255)
                                                                                                                                           attaaaGGGGAGATAAAATCCCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAag            >  1:1243160/4‑150 (MQ=255)
                                                                                                                                           attaaaGGGGAGATAAAATCCCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAag            >  1:383921/4‑150 (MQ=255)
                                                                                                                                             aaaaGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAagag          <  1:1251877/150‑1 (MQ=255)
                                                                                                                                                   ggggATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGct    >  2:872277/1‑150 (MQ=255)
                                                                                                                                                     ggATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTAATATGATTACCTGGCGGAAATTAAACAAATAGAGAGctct  >  1:174974/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCT  >  NC_000913/458589‑458887

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