Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,787,100 C→T 7.2% G5S (GGC→AGC)  ppsA ← phosphoenolpyruvate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,787,1000CT 7.2% 160.2 / 4.4 56G5S (GGC→AGC) ppsAphosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base C (23/29);  new base T (2/2);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

ATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATAAAAAATAAAATTGAAGGTTCA  >  NC_000913/1786829‑1787355
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                               
aTCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCAtt                                                                                                                                                                                                                                                              >  1:63305/1‑275 (MQ=255)
       ggCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGAc                                                                                                                                                                                                                                                       <  2:63304/275‑1 (MQ=255)
         cGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACa                                                                                                                                                                                                                                                      <  1:848855/274‑1 (MQ=255)
           gatgtgctccgccGCTTTCTCAGGCTGAGTAACATCGTCAAGATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTACAGAAACTGGGTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAACAACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGCGCCATTGTTGGACATCg                                                                                                                                                                                                                                                   <  2:1286517/264‑1 (MQ=255)
             tttGCGGGCCCGCTTTCGCAAGCGTGGTAAAATCGTCAATGTCCGTTTTATCCAGCAGTTCATAAATGCGCTGGGTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGATGCATTGTTGACCCTAACCCTGGCTACATCATTCATGCCGAGTTGGTTATACCAAAGTACCAGCGGTGACGAGCCATTGTTGGACATCGaa                                                                                                                                                                                                                                                 <  2:422517/275‑1 (MQ=255)
             tttGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAAGATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGaa                                                                                                                                                                                                                                                 <  1:375493/275‑1 (MQ=255)
                                                tCAATATTCGTTTTATCCAGCAGTTAATACAAGCGCTGTTTTACGCCACATTGGTACAGAAACTGGTTAAACGCGACGGCGGTTGTGGCGAAACCATTCGGAACGGAAACTCCCATTCCGGAAAGATTAGCAATCATGTCACACAGGGAGGCATTTTTGCCCTCAACCCTGTCTACTACATTCATGCCGAGTTGGTTATACCAACACACCAGCGGTGCCGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATCAATGAACGGTTTGCGaaaa                                                                                                                                                                                                              >  2:235851/1‑274 (MQ=255)
                                                        cggtttATCCAGCAGGTCATAAATGGGCGGGTTTACGCCGCTTGGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCAGTCGGAACTGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCt                                                                                                                                                                                                      <  2:1499616/272‑1 (MQ=255)
                                                        cGTTTTGTCCAGCAGGTCATAAATGCGCTGGGTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCt                                                                                                                                                                                                      <  2:559028/275‑1 (MQ=255)
                                                        cGTTTTATCCAGCAGTTAATAAATGCGCTGGTTTACGCCGCTTTGGTCCATAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGATGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCt                                                                                                                                                                                                      <  1:1497189/275‑1 (MQ=255)
                                                        cGTTGTATCCAGCAGTTAATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGTAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATTCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCt                                                                                                                                                                                                      <  1:1564283/275‑1 (MQ=255)
                                                                    cagTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACTTCGAACAATCCTTTTGTGCTAAATGAACGGTTTGAGAAACACTTTTCTGCGCATTTATCt                                                                                                                                                                                          >  1:265310/1‑275 (MQ=255)
                                                                       ttCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAACGCACCAGCGGTGACGAGCCATTGTTGGACAt                                                                                                                                                                                                                                                     >  2:1136306/1‑213 (MQ=255)
                                                                           tAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGaa                                                                                                                                                                                                                                                 <  1:1136307/213‑1 (MQ=255)
                                                                                 cgcTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAAAGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATAATCTCACCCAGGGAGGCATTTTTGCCCCCAAACCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACCCATTTTTGGTCATTTATCTTTGCATATTTnan                                                                                                                                                                             >  2:1173173/1‑274 (MQ=255)
                                                                                      gTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGAt                                                                                                                                                                        >  1:701845/1‑275 (MQ=255)
                                                                                                cTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACAt                                                                                                                                                                                                          >  1:607104/1‑231 (MQ=255)
                                                                                                 tttGGCCCAGAAGCTGGATAAACAAGTGGGCGGTTATGGCAAAACCATTCGGAAGGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTAGATCATTCATGCCGAGTTGGTTATACAAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTGCTGCGCTTTTATCTTTGCATATTTAACCGGATGAAAAAAACg                                                                                                                                                              <  2:265309/274‑1 (MQ=255)
                                                                                                    ggTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCt                                                                                                                                                                                                      <  2:607103/231‑1 (MQ=255)
                                                                                                           aaaCTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCAATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGTTGAAAAAAACGGTAAAACGTTc                                                                                                                                                   >  2:750800/1‑275 (MQ=255)
                                                                                                                 gTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGc                                                                                                                                               >  2:646828/1‑273 (MQ=255)
                                                                                                                           tCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTATGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATAttt                                                                                                                                     >  1:279536/1‑273 (MQ=255)
                                                                                                                           tCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACAAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTATTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATGTTTACAGCAAATTTAttttt                                                                                                                                   >  2:482115/1‑275 (MQ=255)
                                                                                                                              gcggTTGTGGCGAAACCATTCGGAACTGAAACACCCATTCCGGAAAGCTTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCTAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATTAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTATCCGGATGAAACAAACGGTGATTCGTTCAAGCAAATATATTTTTTTa                                                                                                                                <  1:235852/275‑1 (MQ=255)
                                                                                                                              gcggTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTa                                                                                                                                <  1:646829/275‑1 (MQ=255)
                                                                                                                              gcagTTTTGGCCAATCCATTCGGAACGGAAACACCCATGCCGGAAAGATTATTAATCAGTTCACCCATGGGGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCGGCGGTGACGGGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACAGTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTa                                                                                                                                <  2:175627/272‑1 (MQ=255)
                                                                                                                                        cGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGa                                                                                                                      <  1:1101775/275‑1 (MQ=255)
                                                                                                                                          aaaCCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACt                                                                                                                    <  2:701844/275‑1 (MQ=255)
                                                                                                                                                                   ttCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGc                                                                                            >  2:466894/1‑274 (MQ=255)
                                                                                                                                                                    agcggacggcttAGGAATCATTTGCGCCAGGGAGGCATTTTTGACCCCAACCCTGGCTACAATATCCATGCCGCGTTGCTGATGCCCAAGCAACAGTACTGACGGGCCCTGGTTGGTGATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATATTTGCATATGTAACCGGTTGAAAAAAACGGTGAATCGTTCAAGCAAATGTATTTTTTTACTTTTTAAGCCTGATCGCAGCGTTGCGCAAATCTGct                                                                                           <  2:380422/264‑1 (MQ=255)
                                                                                                                                                                                aGTAATCATTTCACCCAGTGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAAtt                                                                              <  1:884454/275‑1 (MQ=255)
                                                                                                                                                                                       aTTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATaaa                                                                                                                                                                                                                               >  1:586345/1‑123 (MQ=255)
                                                                                                                                                                                       aTTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCACCGGTGATGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCACGCAAATATATTTTTTTACTTTTTAAGACTGATGCCAGCGTTGCGCAAATc                                                                                               >  1:380423/1‑251 (MQ=255)
                                                                                                                                                                                           cACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGaa                                                                                                                                                                                                                           <  2:586344/123‑1 (MQ=255)
                                                                                                                                                                                              gattggtatgagacgcagccccAACCCTTGTTACATGATTCATCCCGAATTGTATATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCAt                                                                <  2:1226202/257‑1 (MQ=255)
                                                                                                                                                                                              ccAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGTACAATCCTTTTGTGATAAATGAACTGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGAATGAAAAAAACGGTGAAGCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCAt                                                                <  1:876956/275‑1 (MQ=255)
                                                                                                                                                                                              ccAGGGAGGAATTTTTGCCACCAACCCTGTCTAAATCATTCAGGCCGAGTTGGGTATACCAAAGCAGCAGCCGTGACGAGCCATTGTTGGACATCTAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCGGCTCTCTTCGACAATTTCCACATAAACCAt                                                                <  2:279535/275‑1 (MQ=255)
                                                                                                                                                                                                aGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCg                                                              >  2:1571093/1‑275 (MQ=255)
                                                                                                                                                                                                                    aaCCCTGTTTACATAAGTCATGCCGAATTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATaaaa                                          <  2:352844/275‑1 (MQ=255)
                                                                                                                                                                                                                           tCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCTATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATaaa                                                                    >  1:453265/1‑242 (MQ=255)
                                                                                                                                                                                                                           tCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACCTTTCTGCGCCTTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATTTGCTCTCTTCGACAATTTCCACATaaa                                                                    >  1:1226203/1‑242 (MQ=255)
                                                                                                                                                                                                                              aCATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCAt                                              >  1:352845/1‑261 (MQ=255)
                                                                                                                                                                                                                               catcatTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCTATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCAt                                                                <  2:453264/242‑1 (MQ=255)
                                                                                                                                                                                                                                  cattcatGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTTCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGttt                            <  1:750801/275‑1 (MQ=255)
                                                                                                                                                                                                                                       atGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATa                       <  1:482116/275‑1 (MQ=255)
                                                                                                                                                                                                                                         gCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGAAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATTTGCTCTCTTCGACAATTTCCACATACACCATCGGTATAATCGAAACCATAAATATGAAATGCTGtttt                           >  2:628927/1‑269 (MQ=255)
                                                                                                                                                                                                                                          ccGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATaaaa                    >  2:296448/1‑275 (MQ=255)
                                                                                                                                                                                                                                             aGTTGGTTAAACCAAGTCACCAGCGGTGACTAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGAATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATa                       <  1:628928/269‑1 (MQ=255)
                                                                                                                                                                                                                                               ttGGTTATACCAAAGCACCAGCGGTGACGAGCTATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATccc                                                                                                              >  2:810147/1‑180 (MQ=255)
                                                                                                                                                                                                                                                  gTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGAAGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCAAGCGTTGCGCaaa                                                                                                 >  1:1087982/1‑190 (MQ=255)
                                                                                                                                                                                                                                                   ttATACCAAAGCACCAGCGGTGACGAGCTATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCg                                                                                                          <  1:810148/180‑1 (MQ=255)
                                                                                                                                                                                                                                                      taCCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTg                                                                                             <  2:1087981/190‑1 (MQ=255)
                                                                                                                                                                                                                                                            aGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATAAAAAATAAAATTGAAgg      >  2:46796/1‑271 (MQ=255)
                                                                                                                                                                                                                                                            aGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATAAAAAATAAAATTGAAGGTTc   >  1:858670/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                ccAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGAGGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATAAAAAATAAAATTGAAGGTTCa  <  1:46796/271‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATAAAAAATAAAATTGAAg       >  2:1470441/1‑254 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                               
ATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGAATCGTTCAAGCAAATATATTTTTTTACTTTTTAAGACTGATCCCAGCGTTGCGCAAATCTGCTCTCTTCGACAATTTCCACATAAACCATCGGTATAAACGAAACCATAAAAATGAAATGCTGTTTTCATAAAAAATAAAATTGAAGGTTCA  >  NC_000913/1786829‑1787355

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: