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,9690GT97.6% 643.9 / 6.1 212R1231L (CGT→CTT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base T (84/123);  total (85/127)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 215 total aligned reads displayed.

AGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGAACA  >  CP009273/4178824‑4179113
                                                                                                                                                 |                                                                                                                                                
aGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCttt                                                                                                                                                >  1:2537779/1‑148 (MQ=255)
  cGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTAc                                                                                                                                              >  2:899166/1‑148 (MQ=255)
             aaGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGa                                                                                                                                   <  2:2410532/148‑1 (MQ=255)
              agagATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGaa                                                                                                                                  <  1:1897530/148‑1 (MQ=255)
                  atgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTATTTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTAc                                                                                                                              <  1:1374941/148‑1 (MQ=255)
                  atgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTAc                                                                                                                              <  1:3171643/148‑1 (MQ=255)
                  atgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTAc                                                                                                                              <  1:932380/148‑1 (MQ=255)
                  atgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTAc                                                                                                                              <  1:3174483/148‑1 (MQ=255)
                  atgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTAc                                                                                                                              >  1:255020/1‑148 (MQ=255)
                     atTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAgg                                                                                                                           <  2:2956742/148‑1 (MQ=255)
                     atTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAgg                                                                                                                           <  2:546853/148‑1 (MQ=255)
                     atTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAgg                                                                                                                           <  2:1050588/148‑1 (MQ=255)
                      ttCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGa                                                                                                                          <  1:2984237/148‑1 (MQ=255)
                      ttCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGa                                                                                                                          <  1:897485/148‑1 (MQ=255)
                          gAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGta                                                                                                                      <  2:1804734/148‑1 (MQ=255)
                              tGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATAcc                                                                                                                  >  1:2914514/1‑148 (MQ=255)
                                gCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGt                                                                                                                <  2:1423615/148‑1 (MQ=255)
                                    cAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGc                                                                                                            >  2:2968297/1‑148 (MQ=255)
                                    cAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGc                                                                                                            >  2:731875/1‑148 (MQ=255)
                                        tCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAggg                                                                                                        <  2:1705967/148‑1 (MQ=255)
                                        tCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAggg                                                                                                        <  1:1215701/148‑1 (MQ=255)
                                            cGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGtt                                                                                                    >  1:3569239/1‑148 (MQ=255)
                                             gtgtTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTa                                                                                                   <  2:3615789/148‑1 (MQ=255)
                                               gtTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAg                                                                                                 <  1:608560/148‑1 (MQ=255)
                                               gtTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAg                                                                                                 >  1:3698925/1‑148 (MQ=255)
                                               gtTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAg                                                                                                 <  2:3121969/148‑1 (MQ=255)
                                                ttCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGa                                                                                                <  1:1413490/148‑1 (MQ=255)
                                                ttCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGa                                                                                                <  1:1412870/148‑1 (MQ=255)
                                                ttCGAAGGTGAACGTGTAGAACGTGGTGACGTAAGTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGa                                                                                                <  2:2200643/148‑1 (MQ=255)
                                                 tCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGAt                                                                                               <  2:2952483/148‑1 (MQ=255)
                                                 tCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGAt                                                                                               <  1:619751/148‑1 (MQ=255)
                                                  cGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGAtt                                                                                              <  1:534625/148‑1 (MQ=255)
                                                   gAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTa                                                                                             <  1:3631785/148‑1 (MQ=255)
                                                   gAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTa                                                                                             >  2:2067816/1‑148 (MQ=255)
                                                      ggTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACg                                                                                          <  1:3403877/148‑1 (MQ=255)
                                                      ggTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACg                                                                                          <  1:3469655/148‑1 (MQ=255)
                                                            cGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAac                                                                                    >  2:2661184/1‑148 (MQ=255)
                                                            cGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAac                                                                                    <  2:1509234/148‑1 (MQ=255)
                                                              tgtAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAacac                                                                                  >  1:1916297/1‑148 (MQ=255)
                                                               gtAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAAcaca                                                                                 <  1:3211581/148‑1 (MQ=255)
                                                                tAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACAt                                                                                <  1:2420488/148‑1 (MQ=255)
                                                                   aaCGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGa                                                                             >  2:2352415/1‑148 (MQ=255)
                                                                    aCGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGaa                                                                            >  2:3505878/1‑148 (MQ=255)
                                                                    aCGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGaa                                                                            <  1:429334/148‑1 (MQ=255)
                                                                     cGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAg                                                                           <  1:3629673/148‑1 (MQ=255)
                                                                       tggtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGtt                                                                         <  1:3766998/148‑1 (MQ=255)
                                                                        ggtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTa                                                                        <  1:3682456/148‑1 (MQ=255)
                                                                        ggtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTa                                                                        <  1:1579789/148‑1 (MQ=255)
                                                                         gtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAt                                                                       <  2:3528305/148‑1 (MQ=255)
                                                                             cGTACTGTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtcgtt                                                                   <  1:3347028/148‑1 (MQ=255)
                                                                             cGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtcgtt                                                                   <  1:3346006/148‑1 (MQ=255)
                                                                                aaTTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtcgttcgt                                                                >  1:1774492/1‑148 (MQ=255)
                                                                                aaTTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtcgttcgt                                                                <  2:554344/148‑1 (MQ=255)
                                                                                 aTTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTc                                                               >  2:1629070/1‑148 (MQ=255)
                                                                                   ttCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAg                                                             >  2:3450927/1‑148 (MQ=255)
                                                                                            tCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGt                                                    >  2:1105935/1‑148 (MQ=255)
                                                                                            tCCGGAAGCGACGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGt                                                    >  2:2308651/1‑148 (MQ=255)
                                                                                            tCCGGAAGCGACGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGt                                                    >  2:598083/1‑148 (MQ=255)
                                                                                            tCCAGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGt                                                    >  1:2862004/1‑148 (MQ=255)
                                                                                             ccGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTa                                                   >  2:1404435/1‑148 (MQ=255)
                                                                                             ccGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTa                                                   >  1:174878/1‑148 (MQ=255)
                                                                                              cGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTaa                                                  <  1:1287298/148‑1 (MQ=255)
                                                                                                gAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAg                                                <  1:2847213/148‑1 (MQ=255)
                                                                                                 aaGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGc                                               >  2:297608/1‑148 (MQ=255)
                                                                                                    cgccgcACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTAc                                            >  1:1746564/1‑148 (MQ=255)
                                                                                                    cgccgcACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTAc                                            >  1:308503/1‑148 (MQ=255)
                                                                                                         cACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCg                                       <  2:3713278/148‑1 (MQ=255)
                                                                                                         cACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCg                                       <  1:2076849/148‑1 (MQ=255)
                                                                                                         cACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCg                                       >  1:3668736/1‑148 (MQ=255)
                                                                                                         cACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCg                                       <  2:448278/148‑1 (MQ=255)
                                                                                                           cGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGtt                                     >  2:793850/1‑148 (MQ=255)
                                                                                                             aCATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTaa                                   <  2:1520160/148‑1 (MQ=255)
                                                                                                             aCATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTaa                                   <  1:2472919/148‑1 (MQ=255)
                                                                                                              cATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAAc                                  >  2:1184442/1‑148 (MQ=255)
                                                                                                                ttCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAAcgc                                >  1:315629/1‑148 (MQ=255)
                                                                                                                  ctgcgtctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCgg                              <  2:1916297/148‑1 (MQ=255)
                                                                                                                    gcgtctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGNTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGt                            <  1:1931971/148‑1 (MQ=255)
                                                                                                                      gtctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAg                          <  2:255020/148‑1 (MQ=255)
                                                                                                                      gtctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAg                          <  2:1568237/148‑1 (MQ=255)
                                                                                                                       tctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGc                         >  2:82682/1‑148 (MQ=255)
                                                                                                                       tctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGc                         >  2:219969/1‑148 (MQ=255)
                                                                                                                       tctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGc                         >  1:1080688/1‑148 (MQ=255)
                                                                                                                       tctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGc                         >  2:217951/1‑148 (MQ=255)
                                                                                                                        ctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCt                        <  1:1946193/148‑1 (MQ=255)
                                                                                                                        ctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCt                        <  1:2231399/148‑1 (MQ=255)
                                                                                                                        ctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCt                        >  2:2189409/1‑148 (MQ=255)
                                                                                                                             tggtgTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGAc                   <  1:2384080/148‑1 (MQ=255)
                                                                                                                                   tCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTg             >  1:2635831/1‑148 (MQ=255)
                                                                                                                                     aTGCTGTTACTCTGTACATCGTTAACTAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGa           <  2:1230594/148‑1 (MQ=255)
                                                                                                                                      tGCTGTTACTCTTTACATCGTTAACGAAGTACAGTACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGaa          <  1:251145/148‑1 (MQ=255)
                                                                                                                                      tGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGaa          >  2:387991/1‑148 (MQ=255)
                                                                                                                                      tGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTACGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGaa          >  2:3602420/1‑148 (MQ=255)
                                                                                                                                       gCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAg         >  1:1824984/1‑148 (MQ=255)
                                                                                                                                       gCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAg         <  1:3269182/148‑1 (MQ=255)
                                                                                                                                        cTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAgg        <  1:1521938/148‑1 (MQ=255)
                                                                                                                                        cTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAgg        <  2:917830/148‑1 (MQ=255)
                                                                                                                                        cTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAgg        <  2:92765/148‑1 (MQ=255)
                                                                                                                                          gTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCg      >  1:1347730/1‑148 (MQ=255)
                                                                                                                                           ttACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGa     >  1:1603011/1‑148 (MQ=255)
                                                                                                                                            tACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGaa    >  2:1584143/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGAACA  >  CP009273/4178824‑4179113

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: