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→CAAppkpolyphosphate 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

AGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTGATGCCATTCGCGAACGCGATGTGTTGCTCTATTATCCTTA  >  NC_000913/2623604‑2624083
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
agacgcgcCGCGTCGACGCAAGCCGATGATTCTTCTGGATAAAATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGTTGACCCGCGATGCCGAATCCGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGCGCCGGTGGGTTTTGTTTATCAGCTCGATATGCCCAATg                                                                                                                                                                                                                                       >  2:73201/5‑251 (MQ=255)
aGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATACCCattg                                                                                                                                                                                                                                       >  2:269247/1‑248 (MQ=255)
  aaGCGCCGCGACGACGCAAGCCGATGATTCTTCTGGCTAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGCCGCGCTGACTGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGTCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGAGATATGCCCAATgcg                                                                                                                                                                                                                                     >  2:67902/1‑251 (MQ=255)
        cgcgTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGtt                                                                                                                                                                                                                               >  1:573592/1‑251 (MQ=255)
           gTCGACGCAAGCCGATGATTCTTATGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGGGCGTTTTGTTTATCAACGCGATATGCCCACTGCGCTGGTTGa                                                                                                                                                                                                                             >  2:170529/1‑250 (MQ=255)
                        gatgatTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGaa                                                                                                                                                                                                                >  2:285949/1‑250 (MQ=255)
                                 tctGGATAACATTCTGCGTTACTTCCTGGATGATATTTTCAAAGGCTTCTTTGATTATGACGTGCTGAATGCCTATTCAATTAAGATGTCCCGCGATGCCGAATTCGATTTAGTGCATGATATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGAAGCGTTTAACTGCTTATCCGGTGCGTTTTGTTTAACAACGCGATATGTCCTATGCGCTCGTGGAAGTGGTACGCGAAAAAcgga                                                                                                                                                                                                         <  2:381474/248‑4 (MQ=255)
                                        atcATTCTGCGTTACTTCCTTGATTATAGTTTCAAAGGCTTCTTTTATTATGCCGCGCTTAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCATTCTCAAGCATCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTccc                                                                                                                                                                                                <  2:104381/248‑1 (MQ=255)
                                               tGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGCATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGCTGTCTTACAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGACGTGTTAAGCGACAAACTGACTATTTCCCACTACGAc                                                                                                                                                                                        >  1:292006/1‑251 (MQ=255)
                                                    gtccgactgtatgatATTTTCAAATGCTTCGATGATTAAGTAGCGCTGAATTCCTATGTAAGGAAGCTGACCCGGGATGCCGAATACGATTTAGGGCATGAGATGGAAGCCAGCCTGATGTAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAACCTGTCTATTTCCCGCTACGACTCCa                                                                                                                                                                                    <  2:570401/240‑1 (MQ=255)
                                                     aCTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGCCTCCATc                                                                                                                                                                                  >  2:233496/1‑251 (MQ=255)
                                                        gCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCg                                                                                                                                                                                 <  1:195071/249‑1 (MQ=255)
                                                                          aaGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCAt                                                                                                                                                             >  2:201145/1‑251 (MQ=255)
                                                                                    tGATTATGACTCGCTGAATGCCTATTCAATGAAGATGGCCCTCGATGCCGTATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTaaa                                                                                                                                                    <  2:413858/250‑1 (MQ=255)
                                                                                      aTTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGGGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGAc                                                                                                                                                 <  1:73201/251‑1 (MQ=255)
                                                                                         aTGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACttt                                                                                                                                              >  1:439004/1‑251 (MQ=255)
                                                                                         aTGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACttt                                                                                                                                              >  1:394202/1‑251 (MQ=255)
                                                                                         aTGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCACAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAATCTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACttt                                                                                                                                              <  2:521529/251‑1 (MQ=255)
                                                                                           gACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAATACTttat                                                                                                                                            >  2:113229/1‑251 (MQ=255)
                                                                                                        ccTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAAt                                                                                                                               >  1:100222/1‑251 (MQ=255)
                                                                                                          tATTCAATGAAAACGACCCGCGATGCCGAATAAGAACTAGTGCATGACATGCAAGCCAGACTGATGGCGTTGCTGTCTTCCAGTCTCAAGCAACGTTTACCTGCTGAGCCGGGGCGTTTTGTTTCTCAACGCGATATGCCAAAAGCGCTGGTTGCAGTGTTACGCGCCAAAATGAATATTTCCCGCTACGACTCCATCGTCCCCGGCGGGCGTCACAACAATTTCAAAGAATttattatattcccacattt                                                                                                                             >  2:41508/1‑238 (MQ=255)
                                                                                                           attaaaagAAGATGATCCGGAATGCCGAATACGATTTAGTGCATGAGATGGAAGACAGCCTGATGGAGTGGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGGTTTGTTTATCAACGCGATATTCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCTATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGt                                                                                                                             <  2:160547/243‑1 (MQ=255)
                                                                                                           aTTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAAt                                                                                                                                                                                                                                        >  2:500905/1‑143 (MQ=255)
                                                                                                           aTTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAAt                                                                                                                                                                                                                                        <  1:500905/143‑1 (MQ=255)
                                                                                                           aTTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGCAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCCTCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTTTTAATTTCCCAAATGTc                                                                                                                            >  2:492778/1‑251 (MQ=255)
                                                                                                                  gAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAg                                                                                                                     <  1:6766/251‑1 (MQ=255)
                                                                                                                      aTGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCa                                                                                                                  >  1:519940/1‑250 (MQ=255)
                                                                                                                              cgATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGGTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGa                                                                                                         <  1:630875/251‑1 (MQ=255)
                                                                                                                              cgATGCCGAATACGATTAACGGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTATCAAGCAGCGTTTAACTGCTGAGCCGGTGAGTTTTGTTAATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCACGCTACGTCTCCATCGTCCCCGGCGGCCGTTATCATAATTTTAAAGACTTTTTTAATTTCCCCAATGTCGgccaagccaaaatggtt                                                                                                          >  2:98976/1‑243 (MQ=255)
                                                                                                                                           gagttAGTGCATGAGATGGAGGTCAGCCTGATGGAGTTGATGTCTTCCAGTTACAAGCTGCGTTTAACTGCTGAGCCGGTGCGTTGTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGGAAAACTGACTATTTCCCGCTACTACTCCATCTTCCCCGGCGGTCGTTATCATAATTTGAAAGACTTTATTAATTTCCCCAATGTCGGCAAATCCAATCTGGTGATCAAACCACTGcc                                                                                            <  2:394202/248‑1 (MQ=255)
                                                                                                                                            aTTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCGATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCcg                                                                                           <  1:504945/251‑1 (MQ=255)
                                                                                                                                                 gggCATGAGATGGAAGCCAGCCTGATGGATTTGATTGCTTCCAGCCTCAAGCAGTGTGTAACTGCTGATGTGGTGCGTGTTGTTTATCAACGCGATATGCCCAATTCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTTCCGCGttt                                                                                      <  2:511308/249‑1 (MQ=255)
                                                                                                                                                   gCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTAc                                                                                    <  1:201145/251‑1 (MQ=255)
                                                                                                                                                                   agccTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATaa                                                                                                                                                           >  2:165221/1‑164 (MQ=255)
                                                                                                                                                                   agccTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATaa                                                                                                                                                           <  1:165221/164‑1 (MQ=255)
                                                                                                                                                                          tGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCcg                                                                                                                                                                           >  1:650168/1‑141 (MQ=255)
                                                                                                                                                                          tGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCcg                                                                                                                                                                           <  2:650168/141‑1 (MQ=255)
                                                                                                                                                                                         ccaggatgaaGCAGCGTATAAGTGCTGGGTCGGTGCGTTTTTTGTTACAACGGGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAATCTTACTATTGCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAAGTTTAAAGCCTTTATTTATTTCTCCAGTGTCGGCAAAGCCACAGTGGGGAGCAAACCACCGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCgtgatggt                                              <  2:325407/243‑8 (MQ=255)
                                                                                                                                                                                         ccAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATgg                                               <  1:189281/250‑1 (MQ=255)
                                                                                                                                                                                           aGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTg                                                                     >  1:519146/1‑226 (MQ=255)
                                                                                                                                                                                           aGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGGTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAAGGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGttt                                            <  1:67902/251‑1 (MQ=255)
                                                                                                                                                                                           aGTCTCAAGCAGCGTTTAACTGCGGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTg                                                                     <  2:519146/226‑1 (MQ=255)
                                                                                                                                                                                              ctcAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATCAAGCCCAGTTCCGCAATGGTTTTg                                          >  1:591705/1‑250 (MQ=255)
                                                                                                                                                                                                cAAGCATCGTGTAACTGCGGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACGGCCGCGTTTACGCCATAGTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTGAt                                        <  2:275963/250‑1 (MQ=255)
                                                                                                                                                                                                  agcagcGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATaaa                                                                >  1:26651/1‑224 (MQ=255)
                                                                                                                                                                                                  agcagcGTTTAACGGCTGAGCCGGTGCGTTTTGTTTATCAACGCGATATGCCCAATGCGCTGGTTGAAGTTTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATaaa                                                                <  2:26651/224‑1 (MQ=255)
                                                                                                                                                                                                                                       aTCAACGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTGATGCCATTCGCGAACGCGATGTGTTGCTCTATTATCCTTa  <  2:642817/250‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
AGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCCCGCTACGACTCCATCGTCCCCGGCGGTCGTTATCATAATTTTAAAGACTTTATTAATTTCCCCAATGTCGGCAAAGCCAATCTGGTGAACAAACCACTGCCGCGTTTACGCCATATTTGGTTTGATAAAGCCCAGTTCCGCAATGGTTTTGATGCCATTCGCGAACGCGATGTGTTGCTCTATTATCCTTA  >  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: