Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,888,274 C→A 100% P88Q (CCG→CAG)  chrR → chromate reductase, Class I, flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,888,2740CA98.6% 690.9 / 0.8 233P88Q (CCG→CAG) chrRchromate reductase, Class I, flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (128/100);  total (130/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 233 total aligned reads displayed.

TTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTT  >  CP009273/3888129‑3888418
                                                                                                                                                 |                                                                                                                                                
ttACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAgg                                                                                                                                                >  2:3144876/1‑148 (MQ=255)
  aCCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGt                                                                                                                                              >  1:1837326/1‑148 (MQ=255)
  aCCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCAAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGAAGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTAAAGGGt                                                                                                                                              >  2:2634911/1‑148 (MQ=255)
     atccatTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTggg                                                                                                                                           <  2:834147/148‑1 (MQ=255)
     atccatTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTggg                                                                                                                                           <  2:636777/148‑1 (MQ=255)
     atccatTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTggg                                                                                                                                           <  2:561124/148‑1 (MQ=255)
      tccatTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGGAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGc                                                                                                                                          >  1:2161509/1‑148 (MQ=255)
      tccatTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGTTCCGTCAGGCTGACGTTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGc                                                                                                                                          <  2:2963977/148‑1 (MQ=255)
        catTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGATGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTg                                                                                                                                        >  1:2429563/1‑148 (MQ=255)
         atTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGa                                                                                                                                       >  1:1336391/1‑148 (MQ=255)
         atTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGa                                                                                                                                       >  1:1336434/1‑148 (MQ=255)
           tGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGaaa                                                                                                                                     >  1:1976687/1‑148 (MQ=255)
            gCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGaaaa                                                                                                                                    <  2:3196863/148‑1 (MQ=255)
            gCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGaaaa                                                                                                                                    <  2:2557097/148‑1 (MQ=255)
              cGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAAt                                                                                                                                  <  1:1641562/148‑1 (MQ=255)
               gACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATg                                                                                                                                 >  1:1871870/1‑148 (MQ=255)
               gACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATg                                                                                                                                 >  1:3324453/1‑148 (MQ=255)
                 cagtCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGcc                                                                                                                               <  1:2047107/145‑1 (MQ=255)
                  aTTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCa                                                                                                                              <  2:1478074/148‑1 (MQ=255)
                    tCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATc                                                                                                                            >  1:2307178/1‑148 (MQ=255)
                      ccTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGa                                                                                                                          >  1:1665344/1‑148 (MQ=255)
                      ccTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGa                                                                                                                          <  2:1365772/148‑1 (MQ=255)
                       cTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGAc                                                                                                                         <  1:1397799/148‑1 (MQ=255)
                       cTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTAATCGGTACAGGGTGGGCTGAAAAATGCCATCGAc                                                                                                                         >  2:2286409/1‑148 (MQ=255)
                        ttGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACt                                                                                                                        >  1:9788/1‑148 (MQ=255)
                             tGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCtt                                                                                                                   >  1:2516016/1‑148 (MQ=255)
                               aCGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTc                                                                                                                 >  1:2143922/1‑148 (MQ=255)
                                   tGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATCCCATCGACTGGCTTTCCCGc                                                                                                             <  2:2125349/148‑1 (MQ=255)
                                    gACGTACAGCAGTAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGTAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGcc                                                                                                            <  1:2597060/148‑1 (MQ=255)
                                    gACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGcc                                                                                                            >  1:554455/1‑148 (MQ=255)
                                     aCGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCt                                                                                                           >  1:2567809/1‑148 (MQ=255)
                                          ccgcagGAAGAAGGTTGTCCAGCACCGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGAGGGTGTGGTGATCGTCACGCCGTAATATAACTACACGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCgg                                                                                                      <  1:2221996/146‑1 (MQ=255)
                                          cagcagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGACTATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCgg                                                                                                      >  2:2404189/1‑148 (MQ=255)
                                          cagcagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCgg                                                                                                      >  2:2404136/1‑148 (MQ=255)
                                             cagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATc                                                                                                   >  1:641814/1‑148 (MQ=255)
                                              agtaagaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCTGGCTGACGGTGTGGTGATCGTCACGCCGTAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCa                                                                                                  <  1:1058865/145‑1 (MQ=255)
                                               ggAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCaa                                                                                                 >  2:3420125/1‑148 (MQ=255)
                                                gcagaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAc                                                                                                >  2:1513863/3‑148 (MQ=255)
                                                gaagaaGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAc                                                                                                >  1:2962893/1‑148 (MQ=255)
                                                 aagaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAcc                                                                                               >  1:3424414/1‑148 (MQ=255)
                                                 aagaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAcc                                                                                               >  1:923130/1‑148 (MQ=255)
                                                    aagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCt                                                                                            >  1:97920/1‑148 (MQ=255)
                                                      ggTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTgg                                                                                          >  1:1211235/1‑148 (MQ=255)
                                                          ttCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCCGGCTGACGGTGTGGTGATCGGCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCgg                                                                                      >  1:435450/1‑148 (MQ=255)
                                                                cAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAAcc                                                                                >  2:1960914/1‑148 (MQ=255)
                                                                cAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAAcc                                                                                >  1:1105371/1‑148 (MQ=255)
                                                                  aCGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTTGGTTGATCGTCACGCCGGAATATAACTACTCGGTACAGGTTGGGCTGAAAAATGCAATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCgg                                                                              <  2:640286/148‑1 (MQ=255)
                                                                  aCGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCgg                                                                              >  2:2098117/1‑148 (MQ=255)
                                                                  aCGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCgg                                                                              >  2:2345152/1‑148 (MQ=255)
                                                                  aCGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCgg                                                                              <  2:138494/148‑1 (MQ=255)
                                                                    ggTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTa                                                                            <  2:2224073/148‑1 (MQ=255)
                                                                     gTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTAt                                                                           >  2:380972/1‑148 (MQ=255)
                                                                     gTTGAAGCTCTGGCGGAACAGAGCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTAt                                                                           <  2:2289259/148‑1 (MQ=255)
                                                                      ttGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTCAACCGGTAtt                                                                          >  1:3075961/1‑148 (MQ=255)
                                                                         aaGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGAt                                                                       <  2:3475399/148‑1 (MQ=255)
                                                                          aGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGAtt                                                                      >  2:2562721/1‑148 (MQ=255)
                                                                            ctctGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCa                                                                    >  1:3360923/1‑148 (MQ=255)
                                                                             tctGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAg                                                                   <  2:1715817/148‑1 (MQ=255)
                                                                              attgCGGAACAGATACGTCAGGCTGACGGTGTTGTGATCGTCACGCCGGAATATAAATACTCGGTACAGGGTGGGATGAAAAATGCCATCGACTGTCTTTCACGCCTGCAGGATCAACCGCTGGCCGGTAAACCGTTATTGATTCAGa                                                                  >  2:1617087/4‑148 (MQ=255)
                                                                                ggCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGAcc                                                                <  2:1976672/148‑1 (MQ=255)
                                                                                 gCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCa                                                               >  2:513482/1‑148 (MQ=255)
                                                                                 gCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCa                                                               <  2:1722638/148‑1 (MQ=255)
                                                                                 gCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCa                                                               >  1:911096/1‑148 (MQ=255)
                                                                                       cAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCaa                                                         <  2:2209542/148‑1 (MQ=255)
                                                                                       cAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCaa                                                         <  2:1279982/148‑1 (MQ=255)
                                                                                         gATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATg                                                       >  1:746925/1‑148 (MQ=255)
                                                                                          aTCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATgg                                                      >  1:3026477/1‑148 (MQ=255)
                                                                                            ccGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGc                                                    >  2:1611493/1‑148 (MQ=255)
                                                                                               tCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTg                                                 >  1:2397729/1‑148 (MQ=255)
                                                                                                cAGGCTGACGGTGGGGTGATCGTCACGCCGGAATATAACGACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGa                                                <  1:94128/148‑1 (MQ=255)
                                                                                                 aGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTCAACCGGTATGGAGTCAGACCAGCTCAATGGGCGTGag                                               >  1:2496260/1‑147 (MQ=255)
                                                                                                  ggCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGAtt                                              <  2:3112194/148‑1 (MQ=255)
                                                                                                           tgtgGTGATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGgcgcg                                     <  1:1236876/148‑1 (MQ=255)
                                                                                                                 gATCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGt                               >  2:2812245/1‑148 (MQ=255)
                                                                                                                  aTCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTc                              >  2:208024/1‑148 (MQ=255)
                                                                                                                  aTCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTc                              >  2:818009/1‑148 (MQ=255)
                                                                                                                  aTCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTc                              >  1:2038725/1‑148 (MQ=255)
                                                                                                                  aTCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTc                              >  2:3498590/1‑148 (MQ=255)
                                                                                                                  aTCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCGGTc                              <  2:876163/148‑1 (MQ=255)
                                                                                                                   tCGTCACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCa                             >  2:1996995/1‑148 (MQ=255)
                                                                                                                     gacccGCCGGAATATAACTACGCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGt                           <  1:300211/144‑1 (MQ=255)
                                                                                                                     gTCACGCCGTAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGt                           <  1:1545757/148‑1 (MQ=255)
                                                                                                                       cACGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTAt                         <  2:1409967/148‑1 (MQ=255)
                                                                                                                         cGCCGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCa                       <  2:1336434/148‑1 (MQ=255)
                                                                                                                           ccGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCAcc                     <  2:2429563/148‑1 (MQ=255)
                                                                                                                            cGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCt                    >  1:67358/1‑148 (MQ=255)
                                                                                                                            cGGAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCt                    >  1:3124121/1‑148 (MQ=255)
                                                                                                                              gAATATAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTAAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTgc                  >  2:403412/1‑148 (MQ=255)
                                                                                                                                atatAACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTgcgc                >  1:1775859/1‑148 (MQ=255)
                                                                                                                                  ataACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCa              <  1:3351103/148‑1 (MQ=255)
                                                                                                                                   taACTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAg             <  2:3324453/148‑1 (MQ=255)
                                                                                                                                    aaCTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGa            <  2:3307551/148‑1 (MQ=255)
                                                                                                                                    aaCTACTCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGa            <  2:609210/148‑1 (MQ=255)
                                                                                                                                     actactCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGAt           <  2:1306968/148‑1 (MQ=255)
                                                                                                                                     actactCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGAt           >  1:2014203/1‑148 (MQ=255)
                                                                                                                                      ctactCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGAtt          <  2:2567809/148‑1 (MQ=255)
                                                                                                                                         ctCGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTg       >  2:1221927/1‑148 (MQ=255)
                                                                                                                                           cGGTACAGGGTGGGTTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGt     >  1:3230633/1‑148 (MQ=255)
                                                                                                                                           cGGTACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGt     >  1:2241808/1‑148 (MQ=255)
                                                                                                                                              tACAGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGtttt  <  2:1105371/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTACCATCCATTGCCGACATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTT  >  CP009273/3888129‑3888418

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