Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,181,417 T→G 16.5% C108G (TGT→GGT)  yegV → putative kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,181,4170TG16.5% 163.2 / 5.3 76C108G (TGT→GGT) yegVputative kinase
Reads supporting (aligned to +/- strand):  ref base T (27/34);  new base G (13/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 3.67e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGA  >  NC_000913/2181190‑2181652
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
aCGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTAAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCAATAACCCCGAAGGTTAAAAAGGCTGGTGTCttgggggg                                                                                                                                                                                                                                       >  2:274494/1‑232 (MQ=255)
   agggacagcTCCCCGGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGtt                                                                                                                                                                                                                                    <  1:216971/231‑1 (MQ=255)
            tCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGAGGGCAAAAGAGGGCTTAATAAGCCTGATCGATAACGCCGAAGGGGTTAACGGCTGGGGCCTGGGGTGGGTTGGGCCGGag                                                                                                                                                                                                                           >  1:918818/1‑238 (MQ=255)
              cccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCCGGGCGTGTGGGCGGAGATGATTCGCAACCGGCTGGCAAAAGAGGGCTTAATCAGCCTGACCGTTAACGCCGAAGGGGATAACGGCGGGGGTCTGGCGCTGGttaggcggggggg                                                                                                                                                                                                                         >  2:1059354/1‑228 (MQ=255)
              cccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGGGATAACGGCTGGGGTCTGGGGCTGGGTGGGCCgggggg                                                                                                                                                                                                                         >  2:899090/1‑235 (MQ=255)
              cccTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCCCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAAGCAGCCTGATCGATAACGCCGAAGGGGATAACGGCTGGGGTCTGGGCCCGGTGGGGCCGGAggg                                                                                                                                                                                                                         >  2:632779/1‑239 (MQ=255)
                  ggCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGGGAAAAGGGCTGGGGTCGGGGGGTGGTTGAGCCggggggggga                                                                                                                                                                                                                     >  2:135196/1‑231 (MQ=255)
                   ataaaggacagATATCGAACTGAAACAGCAGGGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAAc                                                                                                                                                                                                                    <  1:547768/229‑1 (MQ=255)
                     gTGGGTGCGATTTCGAACCGGAACAGGAGGGGGTTAACGTTGGCGGCTGCGCCCGGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACg                                                                                                                                                                                                                   <  2:766769/238‑1 (MQ=255)
                     gTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACg                                                                                                                                                                                                                   <  1:899090/238‑1 (MQ=255)
                      tGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGGGTTACCGGCTGGGGGCTGGGGGCGGTTGAGGCGGAGGGctggtgtt                                                                                                                                                                                                                 >  2:216971/1‑232 (MQ=255)
                      tGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGGGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGc                                                                                                                                                                                                                  >  1:606975/1‑238 (MQ=255)
                        ggTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGGTAACGCCGAAGGTGATAACGGCTGGGGTCTGGCGCTGGTTGAGCCGGATGGGGAACGCACt                                                                                                                                                                                                               >  1:450534/1‑239 (MQ=255)
                        ggTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGCTAACGCCGAAGGTGATAACGGCTGGTCGCTGGCGATGTTTGAGCCCGATGGCGAACGCACt                                                                                                                                                                                                               >  1:713656/1‑239 (MQ=255)
                             gATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGACCGGTAACGCCGAAGGGGATAACGGCTGGGGTTTGGGCCTGGTTGGGCCGGAAGGGGAACCCACtttttt                                                                                                                                                                                                          >  1:11017/1‑237 (MQ=255)
                                    aaCTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGACGGTGATAACGGCTGGTGTCTGGCCCTGGGTGAGCCGGGTGGCGAAGACGATTTTATGTGATTc                                                                                                                                                                                                   >  2:679248/1‑239 (MQ=255)
                                      cTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCCCTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAg                                                                                                                                                                                                 <  1:767293/239‑1 (MQ=255)
                                            cagcagAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCCCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGa                                                                                                                                                                                                                      <  1:160436/212‑1 (MQ=255)
                                                       tAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGt                                                                                                                                                                                <  1:135196/239‑1 (MQ=255)
                                                       tAACGTGGGCGGCTGGGCACTGAATATGGCCGGGGCGTTAAAGCGCCTCGGCTCCGAAGCGGGTAATGCCGTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAg                                                                                                                                                                                 <  2:736934/238‑1 (MQ=255)
                                                           gTTGGGGGCTGTGCACTTATTGTTTCCGGGGGGTTAAGGGGCCTGGGCAACGAAGCGGGGAAGGCCTTGCCGGTCGGTCAGGGGGTGTGGGCGGGGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGaa                                                                                                                                                                            <  2:448808/239‑1 (MQ=255)
                                                              ggcggcTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGCTAACGCCGAAGGTGATAACGGCTGGGGTCTGGCGCTGGTTGAGCCGGATGGCGAACCGACTTTTATGTCATTCAGCGGGGTTGAAAATCAGTGGAAta                                                                                                                                                                          >  1:736934/1‑237 (MQ=255)
                                                              ggcggcTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTCAGCCGGATGGCGAACGCCGTTTTATGTCATTCAGCGGGGTTGATAATAAGTgggact                                                                                                                                                                          >  1:160566/1‑234 (MQ=255)
                                                              ggcggcTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGGGTCTGGGGCTGGGTTGGCCGGAGGGCGAACGCACCTTTTTGTCATTCAGCGGGGTTGAAAATAAgggggaag                                                                                                                                                                          >  2:721808/1‑234 (MQ=255)
                                                              ggcggcTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGGGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCCCTTTTATGTCATTCAGCGGGGTTGAAAATCAGTGGAATc                                                                                                                                                                          >  1:625352/1‑238 (MQ=255)
                                                                             aaTTTTGCCGGGGGTTTAAAGGCCCTGGGGATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGcgcg                                                                                                                                                         <  2:606975/239‑1 (MQ=255)
                                                                             aaTATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAAt                                                                                                                                                                                    <  1:362611/213‑1 (MQ=255)
                                                                                       tGGGGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGgcg                                                                                                                                                           <  2:154948/227‑1 (MQ=255)
                                                                                       tGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGgcg                                                                                                                                                           >  1:154948/1‑227 (MQ=255)
                                                                                       tGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGGGTTGAAAATCATGGGAATCGCCCGTGGCTGGGGCGATTAACCGtt                                                                                                                                              >  1:1081650/1‑239 (MQ=255)
                                                                                       tGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGGGTCTGGCGCGGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAAAAGCCAGGGGCTGGGGCGGGTAACCGtt                                                                                                                                              >  2:1005432/1‑239 (MQ=255)
                                                                                       tGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAAAGGCTGGTGTCTGGCGCTGGTCGAACCGGATGGCGAACGCCCTTTTATGCAATTCAGCGGTGTTGAAAATCAGTGGAATCGCCCGTGGTCAGCGCGGTTAAACGtt                                                                                                                                              >  1:1056103/1‑239 (MQ=255)
                                                                                       tGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGGGATAACGGCTGGTGTCTGGGGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGTAAAACAAGGAAAACGCCAGTGGCTGGGGCGGTTAAACGtt                                                                                                                                              >  2:726876/1‑239 (MQ=255)
                                                                                       tGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGGGATAAAGGCTGGTGGCTGGGGCTGGGTGAGCCGGCAGGCGAACGCACTTTTATGTCATTCATGGGGGTTGAAAATCAGGGGAATCGCCAGGGGCTGGCGCGGTTAACCGtt                                                                                                                                              >  2:462191/1‑239 (MQ=255)
                                                                                       tGGCGTTAAAGCGCCCCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGtt                                                                                                                                              <  1:875957/239‑1 (MQ=255)
                                                                                                    ccTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACCTTTATGTCCTTCACCGGTGTTGAAAATCAGGGGAAGCGGCCGTGGCGGGGGGGGTTTAACCGTGGCGccggcggg                                                                                                                                   >  2:468802/1‑232 (MQ=255)
                                                                                                    ccTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGGGGCTGGCGCTGGGTGAGCCGGATGGGGAAAGGCCTTTTAGTTCATTCAACGGGGTTTGAAATCAGGGGAAACGCCACGGGCTGGCGCGCTTAAACCGTTGGGCGTggacgc                                                                                                                                  >  2:987027/1‑235 (MQ=255)
                                                                                                                   gggTAATGCCTTGCCGCTCGGTCAGGTCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTc                                                                                                                    <  1:846166/237‑1 (MQ=255)
                                                                                                                   gggTAATGCCTTGCCGCTCGGTAAGGGGGTGTGGGGGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTc                                                                                                                    <  1:245727/237‑1 (MQ=255)
                                                                                                                       aaTTCCTTGCCGCTGGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGtt                                                                                                               <  1:1005432/238‑1 (MQ=255)
                                                                                                                            cTTGCCGCCCGGTCAGGGCGGGTGGGCGGAGATGATTCGGAACCGGATGGCAAAAGAGGGCTTAAGCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAAc                                                                                                         <  1:721808/239‑1 (MQ=255)
                                                                                                                              tGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGGTGAGCCGGATGGCGGACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCCGGCAGCCTGCGCTATTTTTTCGGGTATAAAATg                                                                                                       >  2:858210/1‑239 (MQ=255)
                                                                                                                               gCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAACTCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATTAACTg                                                                                                       >  1:1048963/1‑238 (MQ=255)
                                                                                                                                                tgtgGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGGTTAACCGTTGCGCCTGGCAGCCTGGCCTATTTTTCCGGTTTTCAACTTGCCTCGCCCTGCGGGGag                                                                                    >  1:913702/1‑238 (MQ=255)
                                                                                                                                                gttgGGGGGGGATGATTCGCAACGGGATGGAAAAAGAGGGTTTAAGCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGTTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGa                                                                                     <  2:1083606/236‑1 (MQ=255)
                                                                                                                                                  tgGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGGGGAATCGCCAGTGGCTGGGCCGATTAACCGTTGCGCCTGGCAGCCCGCCCTATTTTTCCGGTTTTCAACTGGccccccccccgggggaaat                                                                                  >  2:630461/1‑228 (MQ=255)
                                                                                                                                                  ggggCCGAGATGATTCGCAAACGTATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCTCTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAAt                                                                                   <  1:702191/237‑1 (MQ=255)
                                                                                                                                                                   cAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGGGGCGAATTGTTAGTGGAATGGCTgg                                                                >  1:739721/1‑239 (MQ=255)
                                                                                                                                                                              cAAAAGGGGTGTTAAGCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCa                                                      <  1:1030138/238‑1 (MQ=255)
                                                                                                                                                                              cAAAAGGGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCa                                                      <  1:987027/238‑1 (MQ=255)
                                                                                                                                                                              cAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCa                                                      <  1:679248/238‑1 (MQ=255)
                                                                                                                                                                                      ggCTTAATCAGGCTGATGGATAACGCCGAAGGGTATGACGGCTGGTGGCTGGCGCTGGTTGGGCCTGATGGCGAACGCTCTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGa                                             <  1:288499/238‑1 (MQ=255)
                                                                                                                                                                                       gCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGAc                                            >  1:875283/1‑238 (MQ=255)
                                                                                                                                                                                        cTTAATCCGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGACg                                           <  1:632779/238‑1 (MQ=255)
                                                                                                                                                                                        cTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGGTGGCGAACGCACTTTTATGTCATTCAGCGGGGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTTCGCCTGGCAGCCTGCTCTATTTTTCCGGTTTTCAACTGGGCTCGCCCTCGGGCGGATTTTTGGTGGAATGGGTGGAAAAGCGGCAAAACGGGAGGc                                          >  2:639984/1‑239 (MQ=255)
                                                                                                                                                                                                 ccTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCCAGACGTGACGCCGTTTATCg                                 >  1:704647/1‑239 (MQ=255)
                                                                                                                                                                                                     aTCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAAACGGTGCGCCTTGCACCCCGGTCTATTTTTCCCGTTAACAACCGGGCCCGCCCCCGCGGGAATTTTTAGTGGGATGGCGGGGAAAGCGCCAACGCCGTGCGACGGTTATCGAttt                             >  2:181951/1‑239 (MQ=255)
                                                                                                                                                                                                     aTCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGAGGAATCGCCAGTGGCTGGCGCGATTCACCGTTGCGCCTGGCAGCCTGCTTCATTTTTCCGGTTATCACCTGGCCTCGCCCCGCGGCGCATTGCTTGTGGAATGGCTGGGAGAGCTGCAAGAACGGACGCCGTTTTTCGAttt                             >  1:119859/1‑239 (MQ=255)
                                                                                                                                                                                                           agcGCCGGAGGTCATAACGGGTGGTGTCGGGCGGTGGTTGGGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGGCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGACGCCGTTTATCGATTTTGGccc                       <  2:913702/237‑1 (MQ=255)
                                                                                                                                                                                                              gCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTg                                                                                                                                       <  1:654977/127‑1 (MQ=255)
                                                                                                                                                                                                                cGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGGGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGAATAACCGGTGCGCCTGGCAGCCTGCTCTATTTTTCCGGGTATCAACTGGCCTCCCCCTGCGGCGGATTGGTAGTGGGATGGGTGGAAAGGCTGCAAGACGTGAAGCCGGTTTACGCTTTTTGGCCACgga                  >  2:671437/1‑237 (MQ=255)
                                                                                                                                                                                                                 ggagggtgtAACGGTTGGTGTCTGGGGTGGTTTGAGCCGGATGGCGAACGCATTTTTATGTCATTCAGCGGTGTTGAAAATCAGGTGATTCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGACGCCGTTTATCGATTTTGGCCCACGTAt                 <  2:450534/231‑1 (MQ=255)
                                                                                                                                                                                                                 gAAGGTGATAACGGTTGGTTGTTGGCGTTGGTTGAGGCGGATGGGGAACGCACTTTTATGTCATGCAGCGGTTTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGACGCCGTTTATCGATTTTGGCCCACGTAt                 <  2:704647/239‑1 (MQ=255)
                                                                                                                                                                                                                 gAAGGTAATAACGGCTGGGGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGACGCCGTTTATCGATTTTGGCCCACGTAt                 <  1:671437/239‑1 (MQ=255)
                                                                                                                                                                                                                  aaGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGGAGCCTGCTCTATTTTTCCGGGTATCCAATGGCCTCCCCCTCCGGCGAATTGTTTGGGGAATGGCGGGGAGGGGCGCCAAACGGGAGCCCCTTTTTTGTTTTTGGccccccgatt                >  2:553210/1‑232 (MQ=255)
                                                                                                                                                                                                                  aaGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGGTCTATTTTTCCGGTTATCAACTGGCCCCGCCCTGCGGGGAATTGTTTGTGGAATGGCTGGGAAGGGCTGCAAAAGGTGGCCCGGTTTTCCGTTTTGGccccccgaat                >  2:745734/1‑232 (MQ=255)
                                                                                                                                                                                                                    ggTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCAGACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGACGCCGTTTATCGATTTTGGCCCACGTATTg               <  1:1059354/238‑1 (MQ=255)
                                                                                                                                                                                                                      tGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGACGCCGTTTATCGATTTTGGCCCACGTATTGGc             <  1:630461/238‑1 (MQ=255)
                                                                                                                                                                                                                       gATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGGTGGAAGAGCTGCAAACCGGGGCGCCGTTTTTCGGTTTTGGCCCCCGGAATggggg           >  2:975738/1‑236 (MQ=255)
                                                                                                                                                                                                                        aTAACGGCTGGTTTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGACGCCGTTTATCGATTTTGGCCCACGTATTGGCGat          <  1:181951/239‑1 (MQ=255)
                                                                                                                                                                                                                        aTAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGCTTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTTCGGTTTTCAACTGGCCTCGCCCTGCGGCGACTTGTTATTGGGATGGGCTGAAGAGCGGCAACACGGGAAGCCCGTTTTCGGTTTTGGGGCACGGATTGGCGat          >  2:397143/1‑239 (MQ=255)
                                                                                                                                                                                                                        aTAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGACGCCGTTTATCGATTTTGGCCCACGTATTGGCGat          >  1:185392/1‑239 (MQ=255)
                                                                                                                                                                                                                             cgCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATccc     <  1:101184/238‑1 (MQ=255)
                                                                                                                                                                                                                               cTGGGGTCGGGGGCTGGTTGGGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATccc     <  1:639984/237‑1 (MQ=255)
                                                                                                                                                                                                                                 ggTGTTTGGCGCTGGTTGAGCCGGATGGCGAACCCCCTTTTATTGCATTCAGCGGTTTTGAAAATCAGTGGAATCGCCAGGGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGa  <  1:453738/238‑1 (MQ=255)
                                                                                                                                                                                                                                 ggTGTCGTGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTTTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGa  <  1:858210/238‑1 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
ACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGGCATCGAAGCGGGTAATGCCTTGCCGCTCGGTCAGGGCGTGTGGGCGGAGATGATTCGCAACCGGATGGCAAAAGAGGGCTTAATCAGCCTGATCGATAACGCCGAAGGTGATAACGGCTGGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTGACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGA  >  NC_000913/2181190‑2181652

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