Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,178,969 G→T 100% R1231L (CGT→CTT)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,178,9690GT96.5% 514.1 / 9.3 172R1231L (CGT→CTT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base T (76/90);  total (79/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 173 total aligned reads displayed.

AGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGA  >  CP009273/4178824‑4179110
                                                                                                                                                 |                                                                                                                                             
aGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCttt                                                                                                                                             >  2:63734/1‑148 (MQ=255)
   gATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACa                                                                                                                                          <  2:1236821/148‑1 (MQ=255)
     tCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATc                                                                                                                                        >  2:750067/1‑148 (MQ=255)
      ccGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGTCGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCg                                                                                                                                       >  2:1648817/1‑148 (MQ=255)
      ccGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCg                                                                                                                                       >  2:1648548/1‑148 (MQ=255)
      ccGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCg                                                                                                                                       <  2:2404332/148‑1 (MQ=255)
       cGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGt                                                                                                                                      >  2:2253742/1‑148 (MQ=255)
         tACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTa                                                                                                                                    <  2:1386690/148‑1 (MQ=255)
          aCGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTaa                                                                                                                                   <  2:1512958/148‑1 (MQ=255)
          aCGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTaa                                                                                                                                   <  2:749336/148‑1 (MQ=255)
          aCGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTaa                                                                                                                                   >  1:1763356/1‑148 (MQ=255)
             aaGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGa                                                                                                                                >  1:2725294/1‑148 (MQ=255)
              agagATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGaa                                                                                                                               >  1:1731347/1‑148 (MQ=255)
                agaTGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGt                                                                                                                             <  1:3504765/148‑1 (MQ=255)
                 gatgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTACTTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTa                                                                                                                            >  2:1424836/1‑148 (MQ=255)
                 gatgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTa                                                                                                                            <  2:1635330/148‑1 (MQ=255)
                    gatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAg                                                                                                                         >  1:1457770/1‑148 (MQ=255)
                    gatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAg                                                                                                                         >  1:1332997/1‑148 (MQ=255)
                     atTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAgg                                                                                                                        <  2:3097803/148‑1 (MQ=255)
                        ccGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACg                                                                                                                     >  1:2087101/1‑148 (MQ=255)
                           aaaTGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGtat                                                                                                                  <  2:439819/148‑1 (MQ=255)
                           aaaTGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGtat                                                                                                                  <  1:626118/148‑1 (MQ=255)
                           aaaTGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTGTACATCGTTAACGAAGTACAGGACGtat                                                                                                                  <  1:2953614/148‑1 (MQ=255)
                              tGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATAcc                                                                                                               >  1:397327/1‑148 (MQ=255)
                                 cGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTc                                                                                                            <  1:2177622/148‑1 (MQ=255)
                                 cGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTc                                                                                                            <  2:2038231/148‑1 (MQ=255)
                                 cGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTc                                                                                                            >  1:1905015/1‑148 (MQ=255)
                                   tCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTg                                                                                                          <  1:1449382/148‑1 (MQ=255)
                                      gCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAg                                                                                                       <  2:2809130/148‑1 (MQ=255)
                                       cTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAgg                                                                                                      <  1:1640821/148‑1 (MQ=255)
                                         cAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGc                                                                                                    >  1:2000063/1‑148 (MQ=255)
                                         cAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGc                                                                                                    <  2:2523223/148‑1 (MQ=255)
                                          aaCGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCg                                                                                                   <  1:754904/148‑1 (MQ=255)
                                           aCGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGt                                                                                                  >  1:53788/1‑148 (MQ=255)
                                           aCGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGt                                                                                                  <  2:3529608/148‑1 (MQ=255)
                                           aCGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGt                                                                                                  >  1:115057/1‑148 (MQ=255)
                                            cGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGtt                                                                                                 >  2:205678/1‑148 (MQ=255)
                                              tgtTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTaa                                                                                               >  1:1200718/1‑148 (MQ=255)
                                                 tCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGAt                                                                                            >  1:3110428/1‑148 (MQ=255)
                                                 tCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGAt                                                                                            <  2:342356/148‑1 (MQ=255)
                                                   gAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTa                                                                                          <  1:1695952/148‑1 (MQ=255)
                                                    aaGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTaa                                                                                         >  2:991079/1‑148 (MQ=255)
                                                    aaGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTaa                                                                                         <  2:3515867/148‑1 (MQ=255)
                                                          aaCGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATaa                                                                                   <  2:295310/148‑1 (MQ=255)
                                                          aaCGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATaa                                                                                   >  2:1211830/1‑148 (MQ=255)
                                                           aCGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATaaa                                                                                  <  2:332909/148‑1 (MQ=255)
                                                               gtAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAAcaca                                                                              <  1:2742421/148‑1 (MQ=255)
                                                                tAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACAt                                                                             >  2:3276557/1‑148 (MQ=255)
                                                                tAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACAt                                                                             <  1:1014877/148‑1 (MQ=255)
                                                                  gAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCg                                                                           >  1:2314012/1‑148 (MQ=255)
                                                                   aaCGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGa                                                                          >  1:542275/1‑148 (MQ=255)
                                                                   aaCGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGa                                                                          <  2:2780202/148‑1 (MQ=255)
                                                                     cGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAg                                                                        >  1:2574527/1‑148 (MQ=255)
                                                                       tggtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGtt                                                                      <  2:1058468/148‑1 (MQ=255)
                                                                       tggtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGtt                                                                      <  1:782687/148‑1 (MQ=255)
                                                                         gtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAg                                                                        <  1:99760/144‑1 (MQ=255)
                                                                         gtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAg                                                                        <  1:100157/144‑1 (MQ=255)
                                                                          tgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTTGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtc                                                                   <  1:2282541/148‑1 (MQ=255)
                                                                          tgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtc                                                                   <  2:3284638/148‑1 (MQ=255)
                                                                             cGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtcgtt                                                                >  2:2110007/1‑148 (MQ=255)
                                                                                aaTTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtcgttcgt                                                             <  1:1423225/148‑1 (MQ=255)
                                                                                aaTTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtcgttcgt                                                             >  1:3552149/1‑148 (MQ=255)
                                                                                   ttCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAg                                                          <  2:827087/148‑1 (MQ=255)
                                                                                   ttCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAg                                                          <  2:2724138/148‑1 (MQ=255)
                                                                                       gACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGAtgc                                                      <  1:3387768/148‑1 (MQ=255)
                                                                                       gACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGAtgc                                                      >  2:2628602/1‑148 (MQ=255)
                                                                                       gACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGAtgc                                                      >  2:2171838/1‑148 (MQ=255)
                                                                                       gACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGAcgc                                                      <  1:799738/148‑1 (MQ=255)
                                                                                        aCGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGGTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGAtgct                                                     >  2:2472839/1‑148 (MQ=255)
                                                                                             ccGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTa                                                >  2:44920/1‑148 (MQ=255)
                                                                                                gAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAg                                             <  2:371430/148‑1 (MQ=255)
                                                                                                gAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAg                                             >  1:2045643/1‑148 (MQ=255)
                                                                                                 aaGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATTCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGc                                            <  1:422911/148‑1 (MQ=255)
                                                                                                 aaGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGc                                            <  2:1746013/148‑1 (MQ=255)
                                                                                                 aaGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGc                                            >  2:2607059/1‑148 (MQ=255)
                                                                                                 aaGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGc                                            <  2:780100/148‑1 (MQ=255)
                                                                                                  aGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCt                                           <  2:1457770/148‑1 (MQ=255)
                                                                                                      ccgcACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCa                                       <  1:3540096/148‑1 (MQ=255)
                                                                                                          aCGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGt                                   <  1:3350870/148‑1 (MQ=255)
                                                                                                          aCGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGt                                   <  1:3365390/148‑1 (MQ=255)
                                                                                                            gACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTa                                 <  2:1763356/148‑1 (MQ=255)
                                                                                                                 tctgcgtctgcgTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAAcgcg                            <  2:1664309/148‑1 (MQ=255)
                                                                                                                    gcgtctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGt                         >  2:816450/1‑148 (MQ=255)
                                                                                                                       tctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGc                      <  1:1267819/148‑1 (MQ=255)
                                                                                                                        ctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCt                     >  1:3327178/1‑148 (MQ=255)
                                                                                                                            gtggtgTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGa                 <  1:2253742/148‑1 (MQ=255)
                                                                                                                            gtggtgTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGa                 <  1:2522864/148‑1 (MQ=255)
                                                                                                                            gtggtgTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGa                 <  1:2522868/148‑1 (MQ=255)
                                                                                                                             tggtgTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGAc                <  2:966038/148‑1 (MQ=255)
                                                                                                                              ggtgTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACt               <  1:205678/148‑1 (MQ=255)
                                                                                                                                tgtTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTc             >  1:2262741/1‑148 (MQ=255)
                                                                                                                                  ttCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCt           <  1:1999022/148‑1 (MQ=255)
                                                                                                                                  ttCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCt           <  2:663626/148‑1 (MQ=255)
                                                                                                                                   tCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTg          >  1:283574/1‑148 (MQ=255)
                                                                                                                                    cATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTgg         >  1:2441546/1‑148 (MQ=255)
                                                                                                                                      tGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGaa       <  1:2575693/148‑1 (MQ=255)
                                                                                                                                          gTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCg   >  1:3176883/1‑148 (MQ=255)
                                                                                                                                           ttACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGa  <  2:542275/148‑1 (MQ=255)
                                                                                                                                           ttACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGa  <  2:1828290/148‑1 (MQ=255)
                                                                                                                                           ttACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGa  <  2:1905015/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
AGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGA  >  CP009273/4178824‑4179110

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: