Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,932,296 G→T 52.7% 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,2960GT52.7% 43.0 / 306.7 182D131Y (GAT→TAT) pykApyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base G (50/36);  new base T (52/44);  total (102/80)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 183 total aligned reads displayed.

AAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACA  >  CP009273/1932149‑1932440
                                                                                                                                                   |                                                                                                                                                
aaaGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACg                                                                                                                                                  <  1:2325244/148‑1 (MQ=255)
  agaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTAt                                                                                                                                                <  1:1490933/148‑1 (MQ=255)
  agaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTAt                                                                                                                                                >  1:91407/1‑148 (MQ=255)
   gaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATg                                                                                                                                               <  1:3127670/148‑1 (MQ=255)
   gaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATg                                                                                                                                               >  1:1291479/1‑148 (MQ=255)
   gaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATg                                                                                                                                               >  2:3494960/1‑148 (MQ=255)
    aagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATgg                                                                                                                                              <  1:1456399/148‑1 (MQ=255)
      ggCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTc                                                                                                                                            >  2:649311/1‑148 (MQ=255)
      ggCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTc                                                                                                                                            >  1:764717/1‑148 (MQ=255)
        cAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTcgc                                                                                                                                          >  1:2572535/1‑148 (MQ=255)
        cAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTcgc                                                                                                                                          >  1:783829/1‑148 (MQ=255)
           aGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTc                                                                                                                                       >  2:698664/1‑148 (MQ=255)
           aGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTc                                                                                                                                       >  1:126161/1‑148 (MQ=255)
               ttCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGt                                                                                                                                   <  1:2200541/148‑1 (MQ=255)
                   tCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTaaa                                                                                                                               >  1:1856211/1‑148 (MQ=255)
                   tCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTaaa                                                                                                                               <  1:2424142/148‑1 (MQ=255)
                   tCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTaaa                                                                                                                               <  1:683741/148‑1 (MQ=255)
                    cAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTaaaa                                                                                                                              >  2:2510296/1‑148 (MQ=255)
                    cAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTaaaa                                                                                                                              >  2:1760849/1‑148 (MQ=255)
                      atatTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTCCAAAGGCCTGCCGGCTGACGTCGTGCCTGGTGACAGCCTGCTGCTGGCCTATGGTCGCGTCCCGTTAAATGt                                                                                                                            >  2:2287164/1‑148 (MQ=255)
                      atatTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGt                                                                                                                            >  2:2975109/1‑148 (MQ=255)
                       tatTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTa                                                                                                                           >  2:1227037/1‑148 (MQ=255)
                        atTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTAc                                                                                                                          >  1:2050717/1‑148 (MQ=255)
                        atTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTAc                                                                                                                          <  1:1318206/148‑1 (MQ=255)
                        atTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTAc                                                                                                                          <  1:1044795/148‑1 (MQ=255)
                         ttGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt                                                                                                                         >  1:1529775/1‑148 (MQ=255)
                         ttGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt                                                                                                                         >  1:1529772/1‑148 (MQ=255)
                            gggATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGa                                                                                                                      <  1:219620/148‑1 (MQ=255)
                             ggATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGaa                                                                                                                     <  1:3542407/148‑1 (MQ=255)
                              gATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAg                                                                                                                    >  2:1332975/1‑148 (MQ=255)
                              gATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAg                                                                                                                    <  2:588330/148‑1 (MQ=255)
                               aTAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGt                                                                                                                   <  2:2331491/148‑1 (MQ=255)
                               aTAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGt                                                                                                                   >  1:3277668/1‑148 (MQ=255)
                               aTAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGt                                                                                                                   <  1:3094603/148‑1 (MQ=255)
                                      ccTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGc                                                                                                            <  1:3237155/148‑1 (MQ=255)
                                       cTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCa                                                                                                           <  2:91407/148‑1 (MQ=255)
                                        tGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCAt                                                                                                          >  2:1054357/1‑148 (MQ=255)
                                        tGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCAt                                                                                                          >  2:2593925/1‑148 (MQ=255)
                                           tCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGaa                                                                                                       <  2:2728418/148‑1 (MQ=255)
                                            cGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGaaa                                                                                                      >  1:1154119/1‑148 (MQ=255)
                                               cGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAgtg                                                                                                   <  2:1425457/148‑1 (MQ=255)
                                                         ggTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCTTCTACAAAGGCCTGGCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCTCGTCCAGTTAATAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAg                                                                                         <  1:2212313/148‑1 (MQ=255)
                                                         ggTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAg                                                                                         >  2:1675407/1‑148 (MQ=255)
                                                           tAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTc                                                                                       >  1:1138990/1‑148 (MQ=255)
                                                            aaaGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCa                                                                                      >  1:106730/1‑148 (MQ=255)
                                                              aGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCAcc                                                                                    >  2:289610/1‑148 (MQ=255)
                                                              aGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCAcc                                                                                    <  2:3041878/148‑1 (MQ=255)
                                                                gTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGt                                                                                  >  1:2146622/1‑148 (MQ=255)
                                                                       cGACAAAGAAAAAGTCGGTCTCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGATGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCggtggt                                                                           >  2:3098283/1‑148 (MQ=255)
                                                                          cAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTccc                                                                        <  1:2543905/148‑1 (MQ=255)
                                                                          cAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTccc                                                                        >  2:1525939/1‑148 (MQ=255)
                                                                           aaaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTcccc                                                                       <  1:1635973/148‑1 (MQ=255)
                                                                           aaaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTcccc                                                                       >  2:3080011/1‑148 (MQ=255)
                                                                            aaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCct                                                                      <  1:1747944/148‑1 (MQ=255)
                                                                            aaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCct                                                                      <  1:1860175/148‑1 (MQ=255)
                                                                            aaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCct                                                                      >  2:1274982/1‑148 (MQ=255)
                                                                             aGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCctc                                                                     >  2:1236363/1‑148 (MQ=255)
                                                                              gAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCctct                                                                    >  2:1878497/1‑148 (MQ=255)
                                                                              gAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCctct                                                                    <  2:1264814/148‑1 (MQ=255)
                                                                               aaaaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCtctc                                                                   >  1:1688499/1‑148 (MQ=255)
                                                                               aaaaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCtctc                                                                   >  2:1169067/1‑148 (MQ=255)
                                                                                 aaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCca                                                                 >  2:2542532/1‑148 (MQ=255)
                                                                                 aaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCca                                                                 <  2:2206619/148‑1 (MQ=255)
                                                                                 aaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCca                                                                 >  1:2649541/1‑148 (MQ=255)
                                                                                 aaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCca                                                                 <  1:2327704/148‑1 (MQ=255)
                                                                                  aaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCcaa                                                                >  2:2119332/1‑148 (MQ=255)
                                                                                     tCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCcaacaa                                                             >  2:871473/1‑148 (MQ=255)
                                                                                          aTCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAg                                                        <  1:3359798/148‑1 (MQ=255)
                                                                                           tCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAgg                                                       >  1:147757/1‑148 (MQ=255)
                                                                                            cGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGt                                                      <  2:2407408/148‑1 (MQ=255)
                                                                                            cGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGt                                                      <  1:3460168/148‑1 (MQ=255)
                                                                                            cGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGt                                                      <  2:1067640/148‑1 (MQ=255)
                                                                                              aCTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGATAAAAGTACTGGAAGTACAGGGCATGAAAGTCTTCACCGAAGTCACAATCGGTGTTCACCAATCAAACaaaaaacataa                                                    >  2:2161181/1‑143 (MQ=255)
                                                                                                tACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTCTca                                                  >  2:1388532/1‑148 (MQ=255)
                                                                                                  cAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaac                                                >  2:1134850/1‑148 (MQ=255)
                                                                                                  cAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaac                                                >  1:631645/1‑148 (MQ=255)
                                                                                                    aaGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaacaa                                              >  1:2132375/1‑148 (MQ=255)
                                                                                                       gcctgcctGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACt                                           >  2:169663/1‑148 (MQ=255)
                                                                                                           gcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTggc                                       >  2:3039677/1‑148 (MQ=255)
                                                                                                               gctgACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGgcggcg                                   <  1:2447103/148‑1 (MQ=255)
                                                                                                                 tgACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGt                                 <  2:783829/148‑1 (MQ=255)
                                                                                                                 tgACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGt                                 >  1:3160174/1‑148 (MQ=255)
                                                                                                                 tgACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGt                                 <  2:2120362/148‑1 (MQ=255)
                                                                                                                  gACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGtt                                <  2:126161/148‑1 (MQ=255)
                                                                                                                    cgtcgtGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTg                              <  1:1178963/148‑1 (MQ=255)
                                                                                                                    cgtcgtGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTg                              >  1:2807963/1‑148 (MQ=255)
                                                                                                                     gtcgtGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGt                             <  1:1956865/148‑1 (MQ=255)
                                                                                                                               gTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAgc                   >  2:1605964/1‑148 (MQ=255)
                                                                                                                                tGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAgcg                  <  1:667012/148‑1 (MQ=255)
                                                                                                                                tGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAgcg                  <  2:3015572/148‑1 (MQ=255)
                                                                                                                                tGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAgcg                  <  2:2967053/148‑1 (MQ=255)
                                                                                                                                   cATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTg               >  2:2666291/1‑148 (MQ=255)
                                                                                                                                    aTCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGa              <  1:2594532/148‑1 (MQ=255)
                                                                                                                                     tCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGAc             <  1:2999429/148‑1 (MQ=255)
                                                                                                                                     tCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGAc             <  1:2999423/148‑1 (MQ=255)
                                                                                                                                          ctgctgGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGaaa        <  1:1760849/148‑1 (MQ=255)
                                                                                                                                           tgctgGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCAACGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGaaaa       <  1:432768/148‑1 (MQ=255)
                                                                                                                                                gACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACa  >  2:2431664/1‑148 (MQ=255)
                                                                                                                                                gACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACa  <  2:3145823/148‑1 (MQ=255)
                                                                                                                                                gACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACa  <  2:2689549/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACA  >  CP009273/1932149‑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: ‑