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.4% 580.8 / 11.9 196R1231L (CGT→CTT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (0/7);  new base T (83/106);  total (83/113)
Fisher's exact test for biased strand distribution p-value = 2.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-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 201 total aligned reads displayed.

AGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGAACA  >  CP009273/4178824‑4179113
                                                                                                                                                 |                                                                                                                                                
aGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGtt                                                                                                                                                <  2:1562016/148‑1 (MQ=255)
  cGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTAc                                                                                                                                              <  1:1747287/148‑1 (MQ=255)
   gATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACa                                                                                                                                             <  1:1249998/148‑1 (MQ=255)
       cGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGt                                                                                                                                         <  2:1019449/148‑1 (MQ=255)
       cGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGt                                                                                                                                         <  2:3310704/148‑1 (MQ=255)
        gTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGtt                                                                                                                                        >  2:1128391/1‑148 (MQ=255)
        gTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGtt                                                                                                                                        <  1:499502/148‑1 (MQ=255)
         tACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTa                                                                                                                                       <  2:1262970/148‑1 (MQ=255)
          aCGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTaa                                                                                                                                      >  1:2124712/1‑148 (MQ=255)
             aaGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGa                                                                                                                                   <  1:2634762/148‑1 (MQ=255)
              agagATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACTACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGaa                                                                                                                                  >  1:3341411/1‑148 (MQ=255)
              agagATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACTACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGaa                                                                                                                                  >  1:3547620/1‑148 (MQ=255)
              agagATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGaa                                                                                                                                  >  2:1330812/1‑148 (MQ=255)
               gagaTGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAg                                                                                                                                 >  2:1410387/1‑148 (MQ=255)
                 gatgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTa                                                                                                                               <  1:420642/148‑1 (MQ=255)
                  atgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTAc                                                                                                                              >  1:3481449/1‑148 (MQ=255)
                  atgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTAc                                                                                                                              >  1:3299965/1‑148 (MQ=255)
                  atgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTAc                                                                                                                              >  2:1828540/1‑148 (MQ=255)
                  atgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTAc                                                                                                                              >  1:1643339/1‑148 (MQ=255)
                     atTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAgg                                                                                                                           <  1:480095/148‑1 (MQ=255)
                      ttCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGa                                                                                                                          >  1:2738201/1‑148 (MQ=255)
                          gAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGta                                                                                                                      <  1:1527287/148‑1 (MQ=255)
                            aaTGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGtata                                                                                                                    >  2:1396539/1‑148 (MQ=255)
                            aaTGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGtata                                                                                                                    <  1:2042000/148‑1 (MQ=255)
                            aaTGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGtata                                                                                                                    <  2:1913966/148‑1 (MQ=255)
                            aaTGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGtata                                                                                                                    <  2:110924/148‑1 (MQ=255)
                             aTGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATAc                                                                                                                   >  1:814294/1‑148 (MQ=255)
                             aTGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATAc                                                                                                                   >  1:406928/1‑148 (MQ=255)
                                 cGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTc                                                                                                               >  2:1851973/1‑148 (MQ=255)
                                   tCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTg                                                                                                             >  1:226684/1‑148 (MQ=255)
                                     aGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCa                                                                                                           <  1:2585854/148‑1 (MQ=255)
                                       cTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAgg                                                                                                         <  2:3447851/148‑1 (MQ=255)
                                        tCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAggg                                                                                                        <  1:3521699/148‑1 (MQ=255)
                                         cAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGc                                                                                                       >  2:1709040/1‑148 (MQ=255)
                                           aCGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCTGGGCGt                                                                                                     >  2:3097354/1‑148 (MQ=255)
                                            cGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGtt                                                                                                    >  2:2957072/1‑148 (MQ=255)
                                              tgtTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTATCGAAGTACAGGACGTATACCGTCTGCAGGGCGTTaa                                                                                                  >  2:1111175/1‑148 (MQ=255)
                                               gtTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAg                                                                                                 <  2:3493436/148‑1 (MQ=255)
                                               gtTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAg                                                                                                 <  1:1851279/148‑1 (MQ=255)
                                                ttCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGa                                                                                                <  2:1331080/148‑1 (MQ=255)
                                                 tCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGAt                                                                                               <  2:1519225/148‑1 (MQ=255)
                                                    aaGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTaa                                                                                            <  1:302875/148‑1 (MQ=255)
                                                     aGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAAc                                                                                           >  1:3275316/1‑148 (MQ=255)
                                                      ggTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACg                                                                                          <  1:741526/148‑1 (MQ=255)
                                                          aaCGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATaa                                                                                      <  1:865419/148‑1 (MQ=255)
                                                           aCGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATaaa                                                                                     <  1:2399035/148‑1 (MQ=255)
                                                                   aaCGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGa                                                                             <  2:2751068/148‑1 (MQ=255)
                                                                   aaCGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGa                                                                             >  1:767814/1‑148 (MQ=255)
                                                                     cGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAg                                                                           >  1:3156696/1‑148 (MQ=255)
                                                                     cGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAg                                                                           >  2:2145714/1‑148 (MQ=255)
                                                                     cGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAg                                                                           >  2:260112/1‑148 (MQ=255)
                                                                      gtggtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGt                                                                          <  1:298568/148‑1 (MQ=255)
                                                                       tggtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGtt                                                                         <  2:2390810/148‑1 (MQ=255)
                                                                       tggtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGtt                                                                         <  2:669200/148‑1 (MQ=255)
                                                                        ggtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTa                                                                        <  2:1676220/148‑1 (MQ=255)
                                                                         gtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAt                                                                       >  1:2580979/1‑148 (MQ=255)
                                                                            aCGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtcgt                                                                    >  2:341795/1‑148 (MQ=255)
                                                                              gTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtcgttc                                                                  >  1:232890/1‑148 (MQ=255)
                                                                                 aTTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTc                                                               <  2:2605943/148‑1 (MQ=255)
                                                                                    tCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGa                                                            >  2:3187770/1‑148 (MQ=255)
                                                                                    tCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGa                                                            <  2:2351175/148‑1 (MQ=255)
                                                                                      cGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGAtg                                                          >  1:1484611/1‑148 (MQ=255)
                                                                                       gACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGAtgc                                                         >  2:1855788/1‑148 (MQ=255)
                                                                                            tCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGt                                                    >  2:3251519/1‑148 (MQ=255)
                                                                                                 aaGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGc                                               >  1:1018754/1‑148 (MQ=255)
                                                                                                 aaGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGc                                               <  1:1306022/148‑1 (MQ=255)
                                                                                                  aGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCt                                              >  1:722819/1‑148 (MQ=255)
                                                                                                   gcgcCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTa                                             >  2:3059514/1‑148 (MQ=255)
                                                                                                    cgccgcACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTAc                                            >  1:2784088/1‑148 (MQ=255)
                                                                                                    cgccgcACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTAc                                            <  1:2805515/148‑1 (MQ=255)
                                                                                                     gccgcACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTAcc                                           <  2:3481449/148‑1 (MQ=255)
                                                                                                         cACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCg                                       >  1:2503106/1‑148 (MQ=255)
                                                                                                         cACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCg                                       <  2:199010/148‑1 (MQ=255)
                                                                                                          aCGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGt                                      <  2:538277/148‑1 (MQ=255)
                                                                                                           cGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGtt                                     >  1:930638/1‑148 (MQ=255)
                                                                                                           cGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGtt                                     <  1:2907160/148‑1 (MQ=255)
                                                                                                           cGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGtt                                     <  1:364848/148‑1 (MQ=255)
                                                                                                             aCATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTaa                                   <  2:2918934/148‑1 (MQ=255)
                                                                                                                  ctgcgtctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGGTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTAAGATGCTGCGTAAAGATCCCATCGTTAACGCgg                              >  1:1904518/1‑148 (MQ=255)
                                                                                                                   tgcgtctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCggg                             <  1:2745712/148‑1 (MQ=255)
                                                                                                                   tgcgtctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCggg                             <  2:1000689/148‑1 (MQ=255)
                                                                                                                      gtctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAg                          <  2:575375/148‑1 (MQ=255)
                                                                                                                        ctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCt                        >  2:2458086/1‑148 (MQ=255)
                                                                                                                        ctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCt                        >  1:186734/1‑148 (MQ=255)
                                                                                                                            gtggtgTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGa                    <  1:1089112/148‑1 (MQ=255)
                                                                                                                             tggtgTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGAc                   <  2:458002/148‑1 (MQ=255)
                                                                                                                              ggtgTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACt                  >  1:1031730/1‑148 (MQ=255)
                                                                                                                                tgtTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTc                >  2:1960563/1‑148 (MQ=255)
                                                                                                                                tgtTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTc                <  1:3187770/148‑1 (MQ=255)
                                                                                                                                tgtTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTc                <  2:835548/148‑1 (MQ=255)
                                                                                                                                      tGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGaa          <  2:2503062/148‑1 (MQ=255)
                                                                                                                                      tGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGaa          >  2:1438926/1‑148 (MQ=255)
                                                                                                                                       gCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAg         <  1:660682/148‑1 (MQ=255)
                                                                                                                                         tGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGc       <  1:1551964/148‑1 (MQ=255)
                                                                                                                                          gTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCg      <  1:644326/148‑1 (MQ=255)
                                                                                                                                          gTTACTCGTTACATCGTTAACGAAGTTCAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCg      <  2:145796/148‑1 (MQ=255)
                                                                                                                                           ttACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGa     >  1:3058577/1‑148 (MQ=255)
                                                                                                                                           ttACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGa     <  1:1818030/148‑1 (MQ=255)
                                                                                                                                              ctctTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGAACa  >  2:2457476/1‑148 (MQ=255)
                                                                                                                                              ctctTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGAACa  >  1:3140845/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: ‑