Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,932,296 G→T 100% D131Y (GAT→TAT)  pykA → pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,932,2960GT96.1% 541.3 / 10.5 179D131Y (GAT→TAT) pykApyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base G (6/1);  new base T (87/85);  total (93/86)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 179 total aligned reads displayed.

AGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACA  >  CP009273/1932151‑1932440
                                                                                                                                                 |                                                                                                                                                
agaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGAGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTAt                                                                                                                                                <  2:801022/148‑1 (MQ=255)
  aagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATgg                                                                                                                                              <  2:1858992/148‑1 (MQ=255)
  aagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATgg                                                                                                                                              <  2:1263040/148‑1 (MQ=255)
       aaaGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTcgcg                                                                                                                                         >  1:32459/1‑148 (MQ=255)
           ttttCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCa                                                                                                                                     <  2:1534936/148‑1 (MQ=255)
           ttttCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCa                                                                                                                                     <  2:1765251/148‑1 (MQ=255)
               ccTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTa                                                                                                                                 >  2:455349/1‑148 (MQ=255)
               ccTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTa                                                                                                                                 >  2:1313326/1‑148 (MQ=255)
                cTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTaa                                                                                                                                >  2:2067183/1‑148 (MQ=255)
                cTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTaa                                                                                                                                <  1:1495290/148‑1 (MQ=255)
                  cAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTaaaa                                                                                                                              >  1:2054686/1‑148 (MQ=255)
                  cAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTaaaa                                                                                                                              <  1:3157388/148‑1 (MQ=255)
                  cAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTaaaa                                                                                                                              >  1:766795/1‑148 (MQ=255)
                  cAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTcaaa                                                                                                                              >  1:3132887/1‑148 (MQ=255)
                    atatTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGt                                                                                                                            <  1:2133521/148‑1 (MQ=255)
                       ttGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACt                                                                                                                         >  1:143915/1‑148 (MQ=255)
                        tGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTg                                                                                                                        >  1:1733985/1‑148 (MQ=255)
                        tGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTg                                                                                                                        >  1:2567561/1‑148 (MQ=255)
                        tGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTg                                                                                                                        >  1:2392690/1‑148 (MQ=255)
                          gggATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGa                                                                                                                      >  2:2166049/1‑148 (MQ=255)
                            gATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAg                                                                                                                    >  1:2886001/1‑148 (MQ=255)
                            gATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAg                                                                                                                    <  1:1238789/148‑1 (MQ=255)
                            gATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAg                                                                                                                    <  2:1114134/148‑1 (MQ=255)
                             aTAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCTTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGt                                                                                                                   <  2:40416/148‑1 (MQ=255)
                             aTAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGt                                                                                                                   >  2:1085325/1‑148 (MQ=255)
                             aTAAAGTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGt                                                                                                                   <  1:624124/148‑1 (MQ=255)
                               aaaTTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTc                                                                                                                 <  1:2640167/148‑1 (MQ=255)
                                  ttCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAgg                                                                                                              >  2:278804/1‑148 (MQ=255)
                                     cTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACACCGAAACAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGATAAAAGTACTGGAAGTTCAGGGCa                                                                                                           <  2:1841730/148‑1 (MQ=255)
                                        cTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGa                                                                                                        >  1:2692612/1‑148 (MQ=255)
                                         tCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGaa                                                                                                       <  1:1110215/148‑1 (MQ=255)
                                           gACGCCAACCTGTGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAg                                                                                                     <  1:2729032/148‑1 (MQ=255)
                                            aCGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAgt                                                                                                    >  1:398731/1‑148 (MQ=255)
                                             cGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCTGGGCATGAAAgtg                                                                                                   <  2:2095193/148‑1 (MQ=255)
                                              gCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAgtgt                                                                                                  <  2:2097304/148‑1 (MQ=255)
                                                cAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTc                                                                                                <  2:766795/148‑1 (MQ=255)
                                                cAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTc                                                                                                <  2:3015259/148‑1 (MQ=255)
                                                  aCCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCAc                                                                                              >  1:499426/1‑148 (MQ=255)
                                                   ccTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCAcc                                                                                             >  2:1423692/1‑148 (MQ=255)
                                                     tGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGa                                                                                           >  2:2457554/1‑148 (MQ=255)
                                                       ggTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAg                                                                                         <  1:1852206/148‑1 (MQ=255)
                                                       ggTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAg                                                                                         <  1:3403683/148‑1 (MQ=255)
                                                       ggTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAg                                                                                         >  1:1308577/1‑148 (MQ=255)
                                                        gTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGt                                                                                        <  1:3329699/148‑1 (MQ=255)
                                                          aaaGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCa                                                                                      <  2:2981644/148‑1 (MQ=255)
                                                           aaGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCAc                                                                                     >  1:875039/1‑148 (MQ=255)
                                                               tGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGACGTCACCGTc                                                                                 >  1:2129165/1‑148 (MQ=255)
                                                               tGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTc                                                                                 >  2:2608867/1‑148 (MQ=255)
                                                                   ggCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCggtg                                                                             >  1:1401879/1‑148 (MQ=255)
                                                                      gACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTc                                                                          <  1:2136412/148‑1 (MQ=255)
                                                                        cAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTccc                                                                        <  1:1090950/148‑1 (MQ=255)
                                                                          aaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCct                                                                      <  2:1204844/148‑1 (MQ=255)
                                                                            gAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCctct                                                                    <  1:430736/148‑1 (MQ=255)
                                                                             aaaaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCGGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCtctc                                                                   >  2:1628430/1‑148 (MQ=255)
                                                                             aaaaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCtctc                                                                   <  1:1250792/148‑1 (MQ=255)
                                                                             aaaaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCtctc                                                                   >  1:2580696/1‑148 (MQ=255)
                                                                             aaaaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCtctc                                                                   >  1:2580872/1‑148 (MQ=255)
                                                                             aaaaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCtctc                                                                   >  1:2666823/1‑148 (MQ=255)
                                                                                  gTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCcaaca                                                              <  2:1385636/148‑1 (MQ=255)
                                                                                    cGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAAt                                                            >  1:1470804/1‑148 (MQ=255)
                                                                                       tATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATaaa                                                         >  2:1382012/1‑148 (MQ=255)
                                                                                        aTCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAg                                                        <  2:928329/148‑1 (MQ=255)
                                                                                        aTCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAg                                                        >  1:484516/1‑148 (MQ=255)
                                                                                        aTCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAg                                                        >  1:2240209/1‑148 (MQ=255)
                                                                                         tCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAgg                                                       <  2:2908045/148‑1 (MQ=255)
                                                                                          cGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGt                                                      <  1:1951337/148‑1 (MQ=255)
                                                                                          cGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGt                                                      <  2:3278327/148‑1 (MQ=255)
                                                                                            aCTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTAt                                                    >  1:2123005/1‑148 (MQ=255)
                                                                                            aCTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTAt                                                    <  2:1175687/148‑1 (MQ=255)
                                                                                               aCAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaa                                                 >  1:2535174/1‑148 (MQ=255)
                                                                                                  aaGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaacaa                                              >  2:740116/1‑148 (MQ=255)
                                                                                                  aaGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaacaa                                              <  1:2067183/148‑1 (MQ=255)
                                                                                                   aGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACaaa                                             >  1:2593432/1‑148 (MQ=255)
                                                                                                   aGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACaaa                                             >  1:1022065/1‑148 (MQ=255)
                                                                                                    ggCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAAc                                            >  1:1551811/1‑148 (MQ=255)
                                                                                                     gcctgcctGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACt                                           >  2:638419/1‑148 (MQ=255)
                                                                                                     gcctgcctGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACt                                           >  2:2085909/1‑148 (MQ=255)
                                                                                                        tgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTgg                                        >  1:3009161/1‑148 (MQ=255)
                                                                                                        tgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTgg                                        <  1:1144873/148‑1 (MQ=255)
                                                                                                        tgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTgg                                        <  1:1926868/148‑1 (MQ=255)
                                                                                                           ctgctgACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTggcgg                                     >  2:2358912/1‑148 (MQ=255)
                                                                                                               tgACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGt                                 <  2:1202272/148‑1 (MQ=255)
                                                                                                                    tcgtGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTc                            >  2:556124/1‑148 (MQ=255)
                                                                                                                             gTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAgc                   >  2:3226452/1‑148 (MQ=255)
                                                                                                                              tGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAgcg                  <  1:3041427/148‑1 (MQ=255)
                                                                                                                              tGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAgcg                  >  2:2732587/1‑148 (MQ=255)
                                                                                                                               gACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAgcgc                 <  1:2406129/148‑1 (MQ=255)
                                                                                                                                aCATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCt                >  2:2350767/1‑148 (MQ=255)
                                                                                                                                 cATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTg               >  2:2589335/1‑148 (MQ=255)
                                                                                                                                 cATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTg               <  2:517634/148‑1 (MQ=255)
                                                                                                                                 cATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTg               <  2:517927/148‑1 (MQ=255)
                                                                                                                                 cATCCTGCTGCAGGACTATGGTCTCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTg               <  2:3386297/148‑1 (MQ=255)
                                                                                                                                  aTCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGa              <  1:2931798/148‑1 (MQ=255)
                                                                                                                                   tCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGAc             >  1:2947242/1‑148 (MQ=255)
                                                                                                                                   tCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGAc             >  1:2947237/1‑148 (MQ=255)
                                                                                                                                    ccTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGAcc            >  1:2839753/1‑148 (MQ=255)
                                                                                                                                    ccTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGAcc            >  1:1232763/1‑148 (MQ=255)
                                                                                                                                      tgctgctgGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGa          <  2:1990902/148‑1 (MQ=255)
                                                                                                                                         tgctgGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGaaaa       <  2:2567561/148‑1 (MQ=255)
                                                                                                                                           ctgGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAg     <  2:1838107/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACA  >  CP009273/1932151‑1932440

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