Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,178,969 G→T 85.5% 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,9690GT85.5% 609.3 / 121.3 277R1231L (CGT→CTT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (21/19);  new base T (101/135);  total (122/155)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

Only 100 of 280 total aligned reads displayed.

GTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGAACAGG  >  CP009273/4178822‑4179115
                                                                                                                                                   |                                                                                                                                                  
gTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGt                                                                                                                                                   >  2:312441/1‑149 (MQ=255)
    cGATCCGTCCGCAGAGAAGAGTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACa                                                                                                                                               <  2:4786862/149‑1 (MQ=255)
      aTCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATc                                                                                                                                             >  2:391825/1‑149 (MQ=255)
         cgttggAAGAGTTGATTCCGATATGGCTTCAGCTGAACGTGTTCGATGGTGATCGTGTAGAACGTGTTGAAGTTTTTTCCGAATGTCCGGATGAGTAGCAAGACATTCTGCGTCTGCGTGGTGTTCATGATGTTAATCTTTACATCGtt                                                                                                                                          <  1:3140747/144‑1 (MQ=255)
         cGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGtt                                                                                                                                          >  2:4011776/1‑149 (MQ=255)
          gTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTa                                                                                                                                         <  1:2888482/149‑1 (MQ=255)
           tACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTaa                                                                                                                                        <  1:1956996/149‑1 (MQ=255)
               aaGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGaa                                                                                                                                    <  1:3370136/149‑1 (MQ=255)
               aaGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGaa                                                                                                                                    >  1:4336369/1‑149 (MQ=255)
                   gatgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTAc                                                                                                                                >  1:2020225/1‑149 (MQ=255)
                   gatgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTAc                                                                                                                                <  2:2970508/149‑1 (MQ=255)
                         tCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACg                                                                                                                          <  1:907088/149‑1 (MQ=255)
                            gAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGtat                                                                                                                       >  2:41490/1‑149 (MQ=255)
                              aaTGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATAc                                                                                                                     >  2:1675556/1‑149 (MQ=255)
                                 ggCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGt                                                                                                                  <  2:3254136/149‑1 (MQ=255)
                                 ggCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGt                                                                                                                  <  2:4415758/149‑1 (MQ=255)
                                  gCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTc                                                                                                                 <  2:2549994/149‑1 (MQ=255)
                                  gCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTc                                                                                                                 <  1:871556/149‑1 (MQ=255)
                                    gTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTg                                                                                                               >  1:1836713/1‑149 (MQ=255)
                                     tCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGc                                                                                                              <  1:2794708/149‑1 (MQ=255)
                                     tCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGc                                                                                                              <  2:3668912/149‑1 (MQ=255)
                                     tCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATGTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGc                                                                                                              <  2:4512475/149‑1 (MQ=255)
                                           cAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCg                                                                                                        >  2:3315594/1‑149 (MQ=255)
                                           cAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCg                                                                                                        >  2:400890/1‑149 (MQ=255)
                                             aCGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGtt                                                                                                      <  2:1987735/149‑1 (MQ=255)
                                             aCGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGtt                                                                                                      >  2:2464804/1‑149 (MQ=255)
                                               gtgtTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTaa                                                                                                    >  2:2696773/1‑149 (MQ=255)
                                                 gtTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGa                                                                                                  >  1:226456/1‑149 (MQ=255)
                                                 gtTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGa                                                                                                  >  2:2143958/1‑149 (MQ=255)
                                                     gAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTaa                                                                                              >  1:3490956/1‑149 (MQ=255)
                                                      aaGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAAc                                                                                             >  1:1717235/1‑149 (MQ=255)
                                                        ggTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGa                                                                                           <  2:2010937/149‑1 (MQ=255)
                                                         gTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGAt                                                                                          <  1:418961/149‑1 (MQ=255)
                                                          tGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATa                                                                                         >  2:4353620/1‑149 (MQ=255)
                                                          tGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATa                                                                                         <  2:3891339/149‑1 (MQ=255)
                                                             aCGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAac                                                                                      >  2:2091785/1‑149 (MQ=255)
                                                             aCGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAac                                                                                      >  2:2774841/1‑149 (MQ=255)
                                                               gtgtAGAACGTGGTGACGTAATTTCCGCCGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCATTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAacac                                                                                    <  1:3433652/149‑1 (MQ=255)
                                                                tgtAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAAcaca                                                                                   <  1:1299208/149‑1 (MQ=255)
                                                                tgtAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAAcaca                                                                                   <  1:3446846/149‑1 (MQ=255)
                                                                tgtAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAAcaca                                                                                   >  1:2503862/1‑149 (MQ=255)
                                                                tgtAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAAcaca                                                                                   <  2:1708975/149‑1 (MQ=255)
                                                                tgtAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAAcaca                                                                                   <  1:438404/149‑1 (MQ=255)
                                                                 gtAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACAt                                                                                  <  1:366584/149‑1 (MQ=255)
                                                                  tAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATc                                                                                 <  2:3332325/149‑1 (MQ=255)
                                                                    gAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGa                                                                               <  2:120510/149‑1 (MQ=255)
                                                                    gAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGa                                                                               >  1:599000/1‑149 (MQ=255)
                                                                     aaCGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGaa                                                                              >  1:2540124/1‑149 (MQ=255)
                                                                      aCGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAg                                                                             <  2:1836141/149‑1 (MQ=255)
                                                                      aCGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAg                                                                             <  1:94544/149‑1 (MQ=255)
                                                                        gtggtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGtt                                                                           >  2:4066659/1‑149 (MQ=255)
                                                                        gtggtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGtt                                                                           >  2:4318193/1‑149 (MQ=255)
                                                                           gtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtc                                                                        <  1:235636/149‑1 (MQ=255)
                                                                           gtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtc                                                                        >  1:4568188/1‑149 (MQ=255)
                                                                              aCGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtcgtt                                                                     >  2:4664346/1‑149 (MQ=255)
                                                                                gTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtcgttcg                                                                   <  2:3132351/149‑1 (MQ=255)
                                                                                   aTTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCa                                                                >  2:2651969/1‑149 (MQ=255)
                                                                                   aTTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCa                                                                <  1:3333463/149‑1 (MQ=255)
                                                                                    tttCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAg                                                               <  2:3574114/149‑1 (MQ=255)
                                                                                     ttCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGATACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGa                                                              >  2:156478/1‑149 (MQ=255)
                                                                                      tCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGAt                                                             <  2:2196812/149‑1 (MQ=255)
                                                                                      tCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGACTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGAt                                                             <  2:560852/149‑1 (MQ=255)
                                                                                         gACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGAtgct                                                          <  1:159053/149‑1 (MQ=255)
                                                                                         gACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGAtgct                                                          <  2:3729324/149‑1 (MQ=255)
                                                                                          aCGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGAtgctg                                                         >  1:4475659/1‑149 (MQ=255)
                                                                                           cGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGAtgctgc                                                        >  2:484421/1‑149 (MQ=255)
                                                                                            ggTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCg                                                       >  2:1429080/1‑149 (MQ=255)
                                                                                                  gAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGc                                                 <  1:717422/149‑1 (MQ=255)
                                                                                                   aaGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTACCGAAGTACATGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACTCATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCt                                                <  1:2204343/149‑1 (MQ=255)
                                                                                                   aaGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCt                                                <  2:3161386/149‑1 (MQ=255)
                                                                                                   aaGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCt                                                <  1:4602255/149‑1 (MQ=255)
                                                                                                     gcgcCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTAc                                              <  1:3781262/149‑1 (MQ=255)
                                                                                                     gcgcCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTAc                                              <  2:3838564/149‑1 (MQ=255)
                                                                                                       gccgcACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCa                                            <  1:3604911/149‑1 (MQ=255)
                                                                                                       gccgcACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTGAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCa                                            <  1:106354/149‑1 (MQ=255)
                                                                                                         cgcACGACATTCTGCGTCTGCGTGGTGTTCCTGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATc                                          <  1:4637037/149‑1 (MQ=255)
                                                                                                         cgcACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATc                                          <  1:3897250/149‑1 (MQ=255)
                                                                                                           cACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGt                                        <  2:3712051/149‑1 (MQ=255)
                                                                                                              gACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTaa                                     >  2:1294253/1‑149 (MQ=255)
                                                                                                                     tgcgtctgcgtGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGt                              >  2:1880032/1‑149 (MQ=255)
                                                                                                                     tgcgtctgcgtGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGt                              >  2:4269332/1‑149 (MQ=255)
                                                                                                                        gtctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGc                           <  2:3205784/149‑1 (MQ=255)
                                                                                                                        gtctgcgtGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGc                           >  2:2820152/1‑149 (MQ=255)
                                                                                                                         tctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCt                          <  1:512919/149‑1 (MQ=255)
                                                                                                                           tgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTcc                        <  1:259289/149‑1 (MQ=255)
                                                                                                                            gcgtGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCg                       >  1:1105000/1‑149 (MQ=255)
                                                                                                                             cgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGa                      <  1:1843714/149‑1 (MQ=255)
                                                                                                                             cgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGa                      <  1:1843652/149‑1 (MQ=255)
                                                                                                                             cgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGa                      <  2:2246920/149‑1 (MQ=255)
                                                                                                                               tggtgTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACt                    <  2:2819911/149‑1 (MQ=255)
                                                                                                                               tggtgTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACt                    >  1:424742/1‑149 (MQ=255)
                                                                                                                                   gtTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCt                <  1:946688/149‑1 (MQ=255)
                                                                                                                                    ttCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTg               >  2:4275180/1‑149 (MQ=255)
                                                                                                                                    ttCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTg               >  1:3289825/1‑149 (MQ=255)
                                                                                                                                     tCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTgg              >  2:3451595/1‑149 (MQ=255)
                                                                                                                                       aTGCTGTTACTCTTTACATCGTTCACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGaa            >  2:2864655/1‑149 (MQ=255)
                                                                                                                                         gCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAgg          <  2:431615/149‑1 (MQ=255)
                                                                                                                                           tGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCg        >  1:4103240/1‑149 (MQ=255)
                                                                                                                                                cTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGAACAg   <  1:3974926/149‑1 (MQ=255)
                                                                                                                                                 tCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGAACAgg  >  1:603660/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGAACAGG  >  CP009273/4178822‑4179115

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: