Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,623,838 G→A Q265Q (CAG→CAA)  ppk → polyphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,623,8380GA95.7% 135.0 / ‑3.7 47Q265Q (CAG→CAA) ppkpolyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (22/23);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTGATGCCATTCGCGAACGCGATGTGTTGCTCTATTATCCTTA  >  NC_000913/2623604‑2624083
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
agacgcgcCGCGTCGACGCAAGCCGATGATTCTTCTGGATAAAATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGTTGACCCGCGATGCCGAATCCGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGCGCCGGTGGGTTTTGT‑TTATCAGCTCGATATGCCCAATg                                                                                                                                                                                                                                       >  2:73201/5‑251 (MQ=255)
aGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATACCCattg                                                                                                                                                                                                                                       >  2:269247/1‑248 (MQ=255)
  aaGCGCCGCGACGACGCAAGCCGATGATTCTTCTGGCTAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGCCGCGCTGACTGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGTCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGAGATATGCCCAATgcg                                                                                                                                                                                                                                     >  2:67902/1‑251 (MQ=255)
        cgcgTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGtt                                                                                                                                                                                                                               >  1:573592/1‑251 (MQ=255)
           gTCGACGCAAGCCGATGATTCTTATGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGGGCGTTTTGT‑TTATCAACGCGATATGCCCACTGCGCTGGTTGa                                                                                                                                                                                                                             >  2:170529/1‑250 (MQ=255)
                        gatgatTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGaa                                                                                                                                                                                                                >  2:285949/1‑250 (MQ=255)
                                 tctGGATAACATTCTGCGTTACTTCCTGGATGATATTTTCAAAGGCTTCTTTGATTATGACGTGCTGAATGCCTATTCAATTAAGATGTCCCGCGATGCCGAATTCGATTTAGTGCATGATATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGAAGCGTTTAACTGCTTATCCGGTGCGTTTTGT‑TTAACAACGCGATATGTCCTATGCGCTCGTGGAAGTGGTACGCGAAAAAcgga                                                                                                                                                                                                         <  2:381474/248‑4 (MQ=255)
                                        atcATTCTGCGTTACTTCCTTGATTATAGTTTCAAAGGCTTCTTTTATTATGCCGCGCTTAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCATTCTCAAGCATCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTccc                                                                                                                                                                                                <  2:104381/248‑1 (MQ=255)
                                               tGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGCATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGCTGTCTTACAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGACGTGTTAAGCGACAAACTGACTATTTCCCACTACGAc                                                                                                                                                                                        >  1:292006/1‑251 (MQ=255)
                                                    gtccgactgtatgatATTTTCAAATGCTTCGATGATTAAGTAGCGCTGAATTCCTATGTAAGGAAGCTGACCCGGGATGCCGAATACGATTTAGGGCATGAGATGGAAGCCAGCCTGATGTAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAACCTGTCTATTTCCCGCTACGACTCCa                                                                                                                                                                                    <  2:570401/240‑1 (MQ=255)
                                                     aCTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGCCTCCATc                                                                                                                                                                                  >  2:233496/1‑251 (MQ=255)
                                                        gCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCg                                                                                                                                                                                 <  1:195071/249‑1 (MQ=255)
                                                                          aaGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCAt                                                                                                                                                             >  2:201145/1‑251 (MQ=255)
                                                                                    tGATTATGACTCGCTGAATGCCTATTCAATGAAGATGGCCCTCGATGCCGTATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTaaa                                                                                                                                                    <  2:413858/250‑1 (MQ=255)
                                                                                      aTTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGGGTTTTGT‑TTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGAc                                                                                                                                                 <  1:73201/251‑1 (MQ=255)
                                                                                         aTGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACttt                                                                                                                                              >  1:439004/1‑251 (MQ=255)
                                                                                         aTGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACttt                                                                                                                                              >  1:394202/1‑251 (MQ=255)
                                                                                         aTGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCACAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAATCTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACttt                                                                                                                                              <  2:521529/251‑1 (MQ=255)
                                                                                           gACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAATACTttat                                                                                                                                            >  2:113229/1‑251 (MQ=255)
                                                                                                        ccTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAAt                                                                                                                               >  1:100222/1‑251 (MQ=255)
                                                                                                          tATTCAATGAAAACGACCCGCGATGCCGAATAAGAACTAGTGCATGACATGCAAGCCAGACTGATGGCGTTGCTGTCTTCCAGTCTCAAGCAACGTTTACCTGCTGAGCCGGGGCGTTTTGT‑TTCTCAACGCGATATGCCAAAAGCGCTGGTTGCAGTGTTACGCGCCAAAATGAATATTTCCCGCTACGACTCCATCGTCCCCGGCGGGCGTCACAACAATTTCAAAGAATttattatattcccacattt                                                                                                                             >  2:41508/1‑238 (MQ=255)
                                                                                                           attaaaagAAGATGATCCGGAATGCCGAATACGATTTAGTGCATGAGATGGAAGACAGCCTGATGGAGTGGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGGTTTGT‑TTATCAACGCGATATTCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCTATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGt                                                                                                                             <  2:160547/243‑1 (MQ=255)
                                                                                                           aTTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAAt                                                                                                                                                                                                                                        >  2:500905/1‑143 (MQ=255)
                                                                                                           aTTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAAt                                                                                                                                                                                                                                        <  1:500905/143‑1 (MQ=255)
                                                                                                           aTTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGCAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCCTCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTTTTAATTTCCCAAATGTc                                                                                                                            >  2:492778/1‑251 (MQ=255)
                                                                                                                  gAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAg                                                                                                                     <  1:6766/251‑1 (MQ=255)
                                                                                                                      aTGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCa                                                                                                                  >  1:519940/1‑250 (MQ=255)
                                                                                                                              cgATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGGTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGa                                                                                                         <  1:630875/251‑1 (MQ=255)
                                                                                                                              cgATGCCGAATACGATTAACGGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTATCAAGCAGCGTTTAACTGCTGAGCCGGTGAGTTTTGT‑TAATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCACGCTACGTCTCCATCGTCCCCGGCGGCCGTTATCATAATTTTAAAGACTTTTTTAATTTCCCCAATGTCGgccaagccaaaatggtt                                                                                                          >  2:98976/1‑243 (MQ=255)
                                                                                                                                           gagttAGTGCATGAGATGGAGGTCAGCCTGATGGAGTTGATGTCTTCCAGTTACAAGCTGCGTTTAACTGCTGAGCCGGTGCGTTGTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGGAAAACTGACTATTTCCCGCTACTACTCCATCTTCCCCGGCGGTCGTTATCATAATTTGAAAGACTTTATTAATTTCCCCAATGTCGGCAAATCCAATCTGGTGATCAAACCACTGcc                                                                                            <  2:394202/248‑1 (MQ=255)
                                                                                                                                            aTTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCGATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCcg                                                                                           <  1:504945/251‑1 (MQ=255)
                                                                                                                                                 gggCATGAGATGGAAGCCAGCCTGATGGATTTGATTGCTTCCAGCCTCAAGCAGTGTGTAACTGCTGATGTGGTGCGTGTTGT‑TTATCAACGCGATATGCCCAATTCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTTCCGCGttt                                                                                      <  2:511308/249‑1 (MQ=255)
                                                                                                                                                   gCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTAc                                                                                    <  1:201145/251‑1 (MQ=255)
                                                                                                                                                                   agccTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATaa                                                                                                                                                           >  2:165221/1‑164 (MQ=255)
                                                                                                                                                                   agccTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATaa                                                                                                                                                           <  1:165221/164‑1 (MQ=255)
                                                                                                                                                                          tGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCcg                                                                                                                                                                           >  1:650168/1‑141 (MQ=255)
                                                                                                                                                                          tGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCcg                                                                                                                                                                           <  2:650168/141‑1 (MQ=255)
                                                                                                                                                                                         ccaggatgaaGCAGCGTATAAGTGCTGGGTCGGTGCGTTTTTTGTTA‑CAACGGGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAATCTTACTATTGCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAAGTTTAAAGCCTTTATTTATTTCTCCAGTGTCGGCAAAGCCACAGTGGGGAGCAAACCACCGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCgtgatggt                                              <  2:325407/243‑8 (MQ=255)
                                                                                                                                                                                         ccAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATgg                                               <  1:189281/250‑1 (MQ=255)
                                                                                                                                                                                           aGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTg                                                                     >  1:519146/1‑226 (MQ=255)
                                                                                                                                                                                           aGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGGTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAAGGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGttt                                            <  1:67902/251‑1 (MQ=255)
                                                                                                                                                                                           aGTCTCAAGCAGCGTTTAACTGCGGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTg                                                                     <  2:519146/226‑1 (MQ=255)
                                                                                                                                                                                              ctcAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATCAAGCCCAGTTCCGCAATGGTTTTg                                          >  1:591705/1‑250 (MQ=255)
                                                                                                                                                                                                cAAGCATCGTGTAACTGCGGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACGGCCGCGTTTACGCCATAGTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTGAt                                        <  2:275963/250‑1 (MQ=255)
                                                                                                                                                                                                  agcagcGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATaaa                                                                >  1:26651/1‑224 (MQ=255)
                                                                                                                                                                                                  agcagcGTTTAACGGCTGAGCCGGTGCGTTTTGT‑TTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTTTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATaaa                                                                <  2:26651/224‑1 (MQ=255)
                                                                                                                                                                                                                                       aTCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTGATGCCATTCGCGAACGCGATGTGTTGCTCTATTATCCTTa  <  2:642817/250‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
AGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGT‑TTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTGATGCCATTCGCGAACGCGATGTGTTGCTCTATTATCCTTA  >  NC_000913/2623604‑2624083

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