Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,720,305 G→A 8.5% intergenic (+44/+143) cspA → / ← hokA RNA chaperone and antiterminator, cold‑inducible/toxic polypeptide, small

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,720,3050GA 8.5% 390.9 / 21.2 130intergenic (+44/+143)cspA/hokARNA chaperone and antiterminator, cold‑inducible/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base G (52/67);  new base A (2/9);  total (54/76)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

Only 100 of 137 total aligned reads displayed.

TAAAGGTAATACACTATGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAATTCACACAGTGAATCTCTTATCATCCAGGTGAAAAATAAAAGCGTGAAACAAATCACT  >  NC_000913/3720034‑3720568
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                       
taAAGGTAATACACTATGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTATTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCgggg                                                                                                                                                                                                                                                                      >  4:64057/1‑275 (MQ=255)
          aacatATGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCAGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAt                                                                                                                                                                                                                                                            <  4:524891/271‑1 (MQ=255)
              tATGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGACCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCCTTTGGCGGGGATTTTTTTATTTGt                                                                                                                                                                                                                                                        >  2:32557/1‑275 (MQ=255)
               aTGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCCAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtt                                                                                                                                                                                                                                                       >  2:598615/1‑275 (MQ=255)
                tGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAACGGCTTCGCCTTCATCCCTCCTGACGATGGCTCTAAAGATGTGTTAGTACCCTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGAACGAAGGTCCGAAAGTGTCCTTCACCATCTAACGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTCACCAGCCTGTAATCTCTGCTTAAACGCGCAGAATCTAAGTCCCCTGCCTTTCGCAGGTGTTTTTTTTATtttttt                                                                                                                                                                                                                                                      >  4:271142/1‑275 (MQ=255)
                                        tAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatctaaccgt                                                                                                                                                                                                                              >  4:709512/1‑268 (MQ=255)
                                                        tGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCtattat                                                                                                                                                                                                              <  2:709514/275‑1 (MQ=255)
                                                         gACAAAGGCTTCGGCTTCCTCACTCCTCACGGTGGATCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGACGGCCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCTGCAGCTAGTAACGTAACCAGCCTGTAATCTATGCTTAAACGCACAGATGCTAAGATCCCTGCCAATGGGCGGGGATTTTTTTAGTTGGTTTAAGGAAATAAATAATCGATCCGTTAATAAATTCtattttt                                                                                                                                                                                                             >  4:621265/1‑275 (MQ=255)
                                                                                 ccTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTggg                                                                                                                                                                                     <  2:64058/275‑1 (MQ=255)
                                                                                                         ttCGTAAACTTCTCGGCGTTCCAGAACGATGGTTACAAAGCTCTGGACGAAGGTCAGAAAGTGTCCTTAACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCg                                                                                                                                                              <  4:598613/274‑1 (MQ=255)
                                                                                                                            tCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtttt                                                                                                                                                                                                                                                     >  4:755963/1‑168 (MQ=255)
                                                                                                                            tCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGt                                                                                                                                                                                    >  1:60034/1‑233 (MQ=255)
                                                                                                                               aGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGa                                                                                                                                                                                                                                                                     >  8:35610/1‑149 (MQ=255)
                                                                                                                               aGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCaggg                                                                                                                                                                                                                                                                      <  5:838679/148‑1 (MQ=255)
                                                                                                                                gAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAgg                                                                                                                                                                                                                                                 <  2:755965/168‑1 (MQ=255)
                                                                                                                                 aaCGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGAtt                                                                                                                                                                                                                                                                   <  7:891868/149‑1 (MQ=255)
                                                                                                                                 aaCGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTTAAGAATAAATTTGGGTGCAATGAAAATGGGCAACGCCGTAAATCAGGCGGGAATAATTTCccg                                                                                                                                     >  4:326854/1‑275 (MQ=255)
                                                                                                                                 aaCGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAt                                                                                                                                                                                                                                                                    <  6:328268/148‑1 (MQ=255)
                                                                                                                                   cGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtct                                                                                                                                                                                                                                                                  <  5:271930/148‑3 (MQ=255)
                                                                                                                                   cGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtct                                                                                                                                                                                                                                                                  <  6:549260/148‑3 (MQ=255)
                                                                                                                                     aTGGTTACAAATCTCTGGGCGAAGGTCAGAAAATGTCCTTCAACCACGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGAttttt                                                                                                                                                                                                                                                                <  5:350372/148‑1 (MQ=255)
                                                                                                                                      tGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAataaataa                                                                                                                                                                                                                                       >  2:275587/1‑172 (MQ=255)
                                                                                                                                      tGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGAACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCTATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCTCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGaa                                                                                                                                >  1:142687/1‑275 (MQ=255)
                                                                                                                                          tACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAttt                                                                                                                                                                                                                                                          <  7:361517/149‑1 (MQ=255)
                                                                                                                                          tACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcga                                                                                                                                                                                                                                   <  4:275586/172‑1 (MQ=255)
                                                                                                                                           aCAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtcttagattct                                                                                                                                                                                                                                                          <  5:595031/148‑8 (MQ=255)
                                                                                                                                             aaaTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtc                                                                                                                                                                                                                                                                   <  8:841120/137‑2 (MQ=255)
                                                                                                                                             aaaTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtc                                                                                                                                                                                                                                                                   >  7:841120/1‑136 (MQ=255)
                                                                                                                                              aaTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAt                                                                                                                                                                                                                                                                    <  6:402187/135‑1 (MQ=255)
                                                                                                                                              aaTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAt                                                                                                                                                                                                                                                                    >  5:402187/1‑135 (MQ=255)
                                                                                                                                                 ctctGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTc                                                                                                                                                                                                                                                    <  6:138168/148‑1 (MQ=255)
                                                                                                                                                 ctctGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTc                                                                                                                                                                                                                                                    <  6:122998/148‑1 (MQ=255)
                                                                                                                                                    tGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                                                                                                                                                                                                                                                   >  4:433331/1‑146 (MQ=255)
                                                                                                                                                        cGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaa                                                                                                                                                                                                                                               <  2:433333/146‑1 (MQ=255)
                                                                                                                                                               cAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatc                                                                                                                                                                                                                                     <  8:262301/149‑1 (MQ=255)
                                                                                                                                                               cAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatc                                                                                                                                                                                                                                     >  7:140401/1‑149 (MQ=255)
                                                                                                                                                               cAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatc                                                                                                                                                                                                                                     <  7:35610/149‑1 (MQ=255)
                                                                                                                                                                aGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCg                                                                                                                                                                                                                                                                         <  5:573940/112‑1 (MQ=255)
                                                                                                                                                                aGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaa                                                                                                                                                                                                                                              >  8:673176/1‑139 (MQ=255)
                                                                                                                                                                aGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaa                                                                                                                                                                                                                                              <  7:673176/139‑1 (MQ=255)
                                                                                                                                                                aGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGtaat                                                                                                                                                                                                                          <  2:423915/159‑1 (MQ=255)
                                                                                                                                                                    aGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcgatc                                                                                                                                                                                                                                 >  5:137138/1‑148 (MQ=255)
                                                                                                                                                                          cTTCACCATCGGAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGtaat                                                                                                                                                                                                                          >  7:528565/1‑149 (MQ=255)
                                                                                                                                                                              aCCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATaaaa                                                                                                                                                                                                                      >  7:376590/1‑149 (MQ=255)
                                                                                                                                                                                cATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                                                                                                                                                                                                                                                   <  7:328374/118‑1 (MQ=255)
                                                                                                                                                                                cATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                                                                                                                                                                                                                                                   >  8:328374/1‑118 (MQ=255)
                                                                                                                                                                                  tCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtt                                                                                                                                                                                                                                                       <  7:80809/112‑1 (MQ=255)
                                                                                                                                                                                  tCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtt                                                                                                                                                                                                                                                       >  8:80809/1‑112 (MQ=255)
                                                                                                                                                                                  tCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTAtttt                                                                                                                                                                                                        >  2:816710/1‑159 (MQ=255)
                                                                                                                                                                                   cGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCta                                                                                                                                                                                                                  <  5:499802/148‑1 (MQ=255)
                                                                                                                                                                                     agaggggtgcTAAGGGCCGGGGAGCGGGNAAGGTAACCAGCCTGGTATTTGTGCNTAAAAACACAGATTCGAAGTGCCCTGCCATGTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATAGAGCGGGTNATAAAATCTAGTATTAGTTTTGTGAAGAATAAATTTGGGTGAAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCAt                                                                                  <  3:142686/266‑1 (MQ=255)
                                                                                                                                                                                      aaGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtgtg                                                                                                                                                                                                    <  4:816708/159‑1 (MQ=255)
                                                                                                                                                                                       aGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCtattatt                                                                                                                                                                                                             <  7:221888/149‑1 (MQ=255)
                                                                                                                                                                                           gcgcTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTc                                                                                                                                                                                                                                                    >  6:231677/1‑106 (MQ=255)
                                                                                                                                                                                               tAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtgt                                                                                                                                                                                                     <  7:879226/149‑1 (MQ=255)
                                                                                                                                                                                               tAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtg                                                                                                                                                                                                      <  6:477455/148‑1 (MQ=255)
                                                                                                                                                                                                  aGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTTCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGTAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTga                                                                                                                                                                                                   <  6:286961/148‑1 (MQ=255)
                                                                                                                                                                                                  aGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAaataaata                                                                                                                                                                                                                                        >  5:63287/1‑111 (MQ=255)
                                                                                                                                                                                                  aGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAaataaata                                                                                                                                                                                                                                        <  6:63287/111‑1 (MQ=255)
                                                                                                                                                                                                                gtaacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                                                                                                                                                                                                                                                   >  6:281603/1‑86 (MQ=255)
                                                                                                                                                                                                                gtaacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                                                                                                                                                                                                                                                   <  5:281603/86‑1 (MQ=255)
                                                                                                                                                                                                                gtaacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAATTcaca                                                      <  4:359251/275‑1 (MQ=255)
                                                                                                                                                                                                                     gtaacCAGCCTTTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAATTCACACAGTg                                                 <  2:326856/275‑1 (MQ=255)
                                                                                                                                                                                                                     gtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAAGAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCaa                                                                                                                                                                                <  5:337392/148‑1 (MQ=255)
                                                                                                                                                                                                                      taacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAAt                                                                                                                                                                               >  6:754807/1‑148 (MQ=255)
                                                                                                                                                                                                                        acCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCaa                                                                                                                                                                                <  5:52040/145‑1 (MQ=255)
                                                                                                                                                                                                                        acCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCaa                                                                                                                                                                                >  6:52040/1‑145 (MQ=255)
                                                                                                                                                                                                                           aGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAttt                                                                                                                                                                                                                                                          <  6:505483/68‑1 (MQ=255)
                                                                                                                                                                                                                           aGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGaa                                                                                                                                                                          >  6:344773/1‑148 (MQ=255)
                                                                                                                                                                                                                           aGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGaa                                                                                                                                                                          <  5:344773/148‑1 (MQ=255)
                                                                                                                                                                                                                           aGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAAt                                                                                                                                                                         <  7:791207/149‑1 (MQ=255)
                                                                                                                                                                                                                             ccTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATg                                                                                                                                                                        <  6:339894/148‑1 (MQ=255)
                                                                                                                                                                                                                             ccTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATg                                                                                                                                                                        <  5:754807/148‑1 (MQ=255)
                                                                                                                                                                                                                               tGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCaa                                                                                                                                                                                <  5:19084/138‑1 (MQ=255)
                                                                                                                                                                                                                                 tAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTACTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCaa                                                                                                                                                                   <  8:140401/149‑1 (MQ=255)
                                                                                                                                                                                                                                      tctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcgat                                                                                                                                                                                                                                  >  7:233580/1‑81 (MQ=255)
                                                                                                                                                                                                                                      tctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcgat                                                                                                                                                                                                                                  <  8:233580/81‑1 (MQ=255)
                                                                                                                                                                                                                                       ctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGaa                                                                                                                                                                          >  6:368100/1‑136 (MQ=255)
                                                                                                                                                                                                                                       ctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGaa                                                                                                                                                                          <  5:368100/136‑1 (MQ=255)
                                                                                                                                                                                                                                        tGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAATTCACACAGTGAATCTCTTATCATCCAGGt                              >  4:338337/1‑275 (MQ=255)
                                                                                                                                                                                                                                         gCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtgt                                                                                                                                                                                                     <  8:742367/107‑1 (MQ=255)
                                                                                                                                                                                                                                         gCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtgt                                                                                                                                                                                                     >  7:742367/1‑107 (MQ=255)
                                                                                                                                                                                                                                         gCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgtaa                                                                                                                                                           <  8:376590/149‑1 (MQ=255)
                                                                                                                                                                                                                                         gCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgta                                                                                                                                                            <  6:137138/148‑1 (MQ=255)
                                                                                                                                                                                                                                           ttAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatc                                                                                                                                                                                                                                     >  7:55361/1‑73 (MQ=255)
                                                                                                                                                                                                                                           ttAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatc                                                                                                                                                                                                                                     <  8:55361/73‑1 (MQ=255)
                                                                                                                                                                                                                                            tAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATaaa                                                                                                                                                                                           >  5:622809/1‑114 (MQ=255)
                                                                                                                                                                                                                                            tAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATaaa                                                                                                                                                                                           <  6:622809/114‑1 (MQ=255)
                                                                                                                                                                                                                                               aaGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGtaagtaag                                                                                                                                                      >  5:301928/1‑148 (MQ=255)
                                                                                                                                                                                                                                                aGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAgg                                                                                                                                                     >  5:258438/1‑148 (MQ=255)
                                                                                                                                                                                                                                                 gCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATcgc                                                                                                                                                                                                                               <  2:546962/73‑1 (MQ=255)
                                                                                                                                                                                                                                                          cTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTg                                                                                                                                                                                       <  6:850609/104‑1 (MQ=255)
                                                                                                                                                                                                                                                          cTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTg                                                                                                                                                                                       >  5:850609/1‑104 (MQ=255)
                                                                                                                                                                                                                                                            aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTg                                                                                                                                                                                       >  5:485452/1‑102 (MQ=255)
                                                                                                                                                                                                                                                            aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAAttt                                                                                                                                         <  6:258438/148‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAAttt                                                                                                                                         >  6:590875/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                 cccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCANAACACTCCTGTTTAATATGCAATTCACACAGTGAATCCCTTATCATCCAGGGACAAAAAAAAAACGGGAaaaaaaaa     >  3:170694/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                    tGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAAtt                                                          >  2:491143/1‑219 (MQ=255)
                                                                                                                                                                                                                                                                     gCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAATTCACACAGTGAATCTCTTATCATCCAGGTGAAAAATAAAAGCGTGAAACAAATCACt  <  1:170695/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aTTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAATTcaca                                                      <  4:491141/219‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                       
TAAAGGTAATACACTATGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAATTCACACAGTGAATCTCTTATCATCCAGGTGAAAAATAAAAGCGTGAAACAAATCACT  >  NC_000913/3720034‑3720568

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: