Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,261,890 G→A 12.1% intergenic (‑15/+136) prs ← / ← ispE phosphoribosylpyrophosphate synthase/4‑diphosphocytidyl‑2‑C‑methylerythritol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,261,8900GA12.1% 162.1 / 13.5 58intergenic (‑15/+136)prs/ispEphosphoribosylpyrophosphate synthase/4‑diphosphocytidyl‑2‑C‑methylerythritol kinase
Reads supporting (aligned to +/- strand):  ref base G (23/28);  new base A (4/3);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01

GGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCATTCCGGGGCTTGCTCTAGCACCTGGCGGGCTTCAGACTCTGTATCAAATTCAGCAAAGACACAGGCCCCTGTCCCAGTCAGGCGCGA  >  NC_000913/1261608‑1262175
                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                             
ggATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGACATGCg                                                                                                                                                                                                                                                                                         >  1:106974/1‑289 (MQ=255)
           gTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATtggtgg                                                                                                                                                                                                                                                                              >  2:100527/1‑289 (MQ=255)
              aTACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGGCCATTGGTGGaaa                                                                                                                                                                                                                                                                           >  1:479060/1‑289 (MQ=255)
                  gACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAAt                                                                                                                                                                                                                                                                       <  2:342233/289‑1 (MQ=255)
                    ccTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCt                                                                                                                                                                                                                                                                     <  1:266489/289‑1 (MQ=255)
                        cGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCg                                                                                                                                                                                                                                                                 >  2:353392/1‑289 (MQ=255)
                              cACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGaaa                                                                                                                                                                                                                                                              >  1:113268/1‑286 (MQ=255)
                              cACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACt                                                                                                                                                                                                                                                           >  2:177246/1‑289 (MQ=255)
                                 gacgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACtgtg                                                                                                                                                                                                                                                        <  2:222604/289‑1 (MQ=255)
                                 gacgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACtgtg                                                                                                                                                                                                                                                        <  1:100527/289‑1 (MQ=255)
                                 gacgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACtgtg                                                                                                                                                                                                                                                        <  1:156415/289‑1 (MQ=255)
                                                          tAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAgg                                                                                                                                                                                                                               >  2:141382/1‑289 (MQ=255)
                                                          tAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAgg                                                                                                                                                                                                                               <  1:445702/289‑1 (MQ=255)
                                                          tAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGACATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAgg                                                                                                                                                                                                                               <  2:106974/289‑1 (MQ=255)
                                                            aTTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCg                                                                                                                                                                                                                             <  1:47435/289‑1 (MQ=255)
                                                            aTTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCg                                                                                                                                                                                                                             <  2:197649/289‑1 (MQ=255)
                                                            aTTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACATACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCg                                                                                                                                                                                                                             <  2:477991/289‑1 (MQ=255)
                                                            aTTCCATCAGGTTGTCGGTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCg                                                                                                                                                                                                                             <  1:464791/289‑1 (MQ=255)
                                                              tCCATCAGGTTGTCGGTAGTAGGGGCACAAGTGGACTGGATGATGAAAAGATCGCCACCGCGTACATTTTCAGTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCATGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGgt                                                                                                                                                                                                                           <  2:289689/289‑1 (MQ=255)
                                                                               agtagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGc                                                                                                                                                                                                                              >  2:418675/1‑269 (MQ=255)
                                                                               agtagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGc                                                                                                                                                                                                                              <  1:418675/269‑1 (MQ=255)
                                                                               agtagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGAt                                                                                                                                                                                                          <  2:113268/289‑1 (MQ=255)
                                                                                  agGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCa                                                                                                                                                                                                       <  1:177246/289‑1 (MQ=255)
                                                                                    gggCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCggg                                                                                                                                                                                                                                                    <  1:481520/242‑1 (MQ=255)
                                                                                    gggCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCggg                                                                                                                                                                                                                                                    >  2:481520/1‑242 (MQ=255)
                                                                                    gggCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAAc                                                                                                                                                                                                     >  1:298340/1‑289 (MQ=255)
                                                                                           aGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAAt                                                                                                                                                                                                                                                                       >  1:400395/1‑216 (MQ=255)
                                                                                           aGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAAt                                                                                                                                                                                                                                                                       <  2:400395/216‑1 (MQ=255)
                                                                                               gACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGcc                                                                                                                                                                                                                                                  <  1:363486/233‑1 (MQ=255)
                                                                                               gACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGcc                                                                                                                                                                                                                                                  >  2:363486/1‑233 (MQ=255)
                                                                                                    gatgatGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGAt                                                                                                                                                                                                                                         >  1:131057/1‑237 (MQ=255)
                                                                                                    gatgatGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGAt                                                                                                                                                                                                                                         <  2:131057/237‑1 (MQ=255)
                                                                                                              atatCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATTCGGCACGAGAAGAACCACAGGCATCGGTCCATTGGTGGAAAGAATCTGCTGAAAATTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGAGTATTAATGGGCTCGATGCAGCGTCGTGAACAAGGGTCCGTGGTCAcc                                                                                                                                                                           >  2:482689/1‑289 (MQ=255)
                                                                                                                 tCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGaa                                                                                                                                                                        >  2:281262/1‑289 (MQ=255)
                                                                                                                       ccGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCa                                                                                                                                                                   >  1:136967/1‑288 (MQ=255)
                                                                                                                             aCATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGcc                                                                                                                                                            <  2:298340/289‑1 (MQ=255)
                                                                                                                                     aTTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTaa                                                                                                                                                    >  2:39267/1‑289 (MQ=255)
                                                                                                                                                           ccATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCCGGCCCGAGAAGAACCTCAGACATGCGTTCATTGGTTGAAAGAATCTGCCGAAAACTGTGCGGGCCCGGCATGACCCTTTCCAGGCGGTGTCTTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAATCTCAGCTTGCCCGCCTTCACGCATGGCTCTTTGCACTgggg                                                                                                                              >  2:444429/1‑289 (MQ=255)
                                                                                                                                                            cATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGa                                                                                                                             <  1:141382/289‑1 (MQ=255)
                                                                                                                                                                            tACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTc                                                                                                             >  1:257033/1‑289 (MQ=255)
                                                                                                                                                                                  ggcgTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCcaca                                                                                                       <  1:39267/289‑1 (MQ=255)
                                                                                                                                                                                      tCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGTGAAGAACCTCAGACATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAg                                                                                                    >  2:239845/1‑288 (MQ=255)
                                                                                                                                                                                        gCCGAGTGTAGTGTCCAGGCGGTGGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACAAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCGGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCa                                                                                                 <  1:482689/289‑1 (MQ=255)
                                                                                                                                                                                               gAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGACATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAg                                                                                            <  1:239845/287‑1 (MQ=255)
                                                                                                                                                                                                             gTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGACATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCATTCCGGGGCTTg                                                                             >  2:210734/1‑288 (MQ=255)
                                                                                                                                                                                                                          tgtgCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCATTCCGGGGCTTGCTCTAGCACCTGGc                                                               >  2:479837/1‑289 (MQ=255)
                                                                                                                                                                                                                                 gTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCa                                                                                                                                                <  1:446754/201‑1 (MQ=255)
                                                                                                                                                                                                                                 gTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCa                                                                                                                                                >  2:446754/1‑201 (MQ=255)
                                                                                                                                                                                                                                   tCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGACATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCATTCCGGGGCTTGCTCTAGCACCTGGCGGGCTTCa                                                       <  1:210734/288‑1 (MQ=255)
                                                                                                                                                                                                                                      ggggTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCg                                                                                                                                                                                                                             <  1:338583/119‑1 (MQ=255)
                                                                                                                                                                                                                                      ggggTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCg                                                                                                                                                                                                                             >  2:338583/1‑119 (MQ=255)
                                                                                                                                                                                                                                      ggggTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCATTCCGGGGCTTGCTCTAGCACCTGGCGGGCTTCAGAct                                                   <  1:479837/289‑1 (MQ=255)
                                                                                                                                                                                                                                           ggCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCg                                                                                                                                                                                                                             >  1:17050/1‑114 (MQ=255)
                                                                                                                                                                                                                                           ggCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCg                                                                                                                                                                                                                             <  2:17050/114‑1 (MQ=255)
                                                                                                                                                                                                                                                                    tCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCATTCCGGGGCTTGCTCTAGCACCTGGCGGGCTTCAGACTCTGTATCAAATTCAGCAAAGACACAGGccc                     <  1:281262/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 gaaCCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCATTCCGGGGCTTGCTCTAGCACCTGGCGGGCTTCAGAct                                                   >  2:465976/1‑246 (MQ=255)
                                                                                                                                                                                                                                                                                 gaaCCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCATTCCGGGGCTTGCTCTAGCACCTGGCGGGCTTCAGAct                                                   <  1:465976/246‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        aGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCATTCCGGGGCTTGCTCTAGCACCTGGCGGGCTTCAGACTCTGTATCAAATTCAGCAAAGACACAGGCCCCTGTCCCAGTCAGGCGCGa  >  2:171352/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                             
GGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCATTCCGGGGCTTGCTCTAGCACCTGGCGGGCTTCAGACTCTGTATCAAATTCAGCAAAGACACAGGCCCCTGTCCCAGTCAGGCGCGA  >  NC_000913/1261608‑1262175

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