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,9690GT100.0% 646.5 / NA 201R1231L (CGT→CTT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (77/124);  total (77/124)
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 202 total aligned reads displayed.

GCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGAACA  >  CP009273/4178825‑4179113
                                                                                                                                                |                                                                                                                                                
gCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGGTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTa                                                                                                                                               >  2:1332680/1‑148 (MQ=255)
 cGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGACCGTGGTGNCGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTAc                                                                                                                                              <  1:1785686/148‑1 (MQ=255)
   aTCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACAt                                                                                                                                            >  2:663521/1‑148 (MQ=255)
   aTCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACAt                                                                                                                                            <  1:3048958/148‑1 (MQ=255)
     ccGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCg                                                                                                                                          >  2:1814332/1‑148 (MQ=255)
     ccGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCg                                                                                                                                          <  2:1981604/148‑1 (MQ=255)
       gTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGtt                                                                                                                                        >  1:590021/1‑148 (MQ=255)
        tACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTa                                                                                                                                       >  2:2950998/1‑148 (MQ=255)
         aCGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTaa                                                                                                                                      >  2:1991234/1‑148 (MQ=255)
         aCGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTaa                                                                                                                                      >  2:217024/1‑148 (MQ=255)
           gAAGAGATGATTCCGAATTGGCGTCAGCTCATCGTGTTAGAAGGTGAACGTGTTGATCGTGGTGACGTAATTTCGGACGGTCCGGAAGCTCCGAACGAAATTCTTCGTCTGCGTGGTGTTCATGCTGTTACTATTTACATCGTTAACg                                                                                                                                    <  1:362635/148‑1 (MQ=255)
            aaGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGa                                                                                                                                   <  1:1456743/148‑1 (MQ=255)
             agagATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGaa                                                                                                                                  <  2:2620712/148‑1 (MQ=255)
              gagaTGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAg                                                                                                                                 <  2:1119188/148‑1 (MQ=255)
               agaTGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGt                                                                                                                                >  2:1317505/1‑148 (MQ=255)
                gatgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTa                                                                                                                               <  2:2548529/148‑1 (MQ=255)
                 atgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTAc                                                                                                                              <  2:2189416/148‑1 (MQ=255)
                 atgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTAc                                                                                                                              <  1:596589/148‑1 (MQ=255)
                 atgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTAc                                                                                                                              <  1:139820/148‑1 (MQ=255)
                  tgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACa                                                                                                                             <  2:1190263/148‑1 (MQ=255)
                  tgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACa                                                                                                                             <  1:826968/148‑1 (MQ=255)
                   gatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAg                                                                                                                            >  2:2548906/1‑148 (MQ=255)
                    atTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAgg                                                                                                                           >  1:1310371/1‑148 (MQ=255)
                       ccGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACg                                                                                                                        <  1:3144679/148‑1 (MQ=255)
                        cGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGt                                                                                                                       <  2:849610/148‑1 (MQ=255)
                            aTGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATAc                                                                                                                   >  2:2778626/1‑148 (MQ=255)
                            aTGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATAc                                                                                                                   >  1:1506400/1‑148 (MQ=255)
                             tGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATAcc                                                                                                                  >  2:2273709/1‑148 (MQ=255)
                               gCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTGTACCGt                                                                                                                >  2:2113852/1‑148 (MQ=255)
                               gCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGt                                                                                                                <  2:2626174/148‑1 (MQ=255)
                                    aGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTTATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCa                                                                                                           <  2:304011/148‑1 (MQ=255)
                                    aGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCa                                                                                                           <  1:2823191/148‑1 (MQ=255)
                                    aGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCa                                                                                                           <  2:2071641/148‑1 (MQ=255)
                                    aGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCa                                                                                                           <  1:412087/148‑1 (MQ=255)
                                     gCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAg                                                                                                          <  1:1665068/148‑1 (MQ=255)
                                     gCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAg                                                                                                          <  1:3403772/148‑1 (MQ=255)
                                      cTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAgg                                                                                                         >  2:3467991/1‑148 (MQ=255)
                                      cTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAgg                                                                                                         >  2:3467962/1‑148 (MQ=255)
                                         aaCGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCg                                                                                                      <  2:1393027/148‑1 (MQ=255)
                                         aaCGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCg                                                                                                      <  1:2321722/148‑1 (MQ=255)
                                          aCGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGt                                                                                                     <  2:1078853/148‑1 (MQ=255)
                                           cGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGtt                                                                                                    <  2:2412606/148‑1 (MQ=255)
                                            gtgtTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTa                                                                                                   <  1:2117206/148‑1 (MQ=255)
                                               ttCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGa                                                                                                <  2:278305/148‑1 (MQ=255)
                                               ttCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGa                                                                                                <  2:278865/148‑1 (MQ=255)
                                                tCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGAt                                                                                               <  2:1504608/148‑1 (MQ=255)
                                                  gAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTa                                                                                             <  2:1462175/148‑1 (MQ=255)
                                                  gAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTa                                                                                             <  2:173110/148‑1 (MQ=255)
                                                  gAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTa                                                                                             <  2:172657/148‑1 (MQ=255)
                                                   aaGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTaa                                                                                            <  1:3139809/148‑1 (MQ=255)
                                                      gTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGa                                                                                         >  1:1935228/1‑148 (MQ=255)
                                                         aaCGTGTTGAACGTGGTGACGTAATTTCCGACGTTCCGGAAGCGCCGCACTACTTTCTTCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATTCAGTCTGCAGTGCGTTAAGATTATCGATat                                                                                      <  2:485240/148‑2 (MQ=255)
                                                         aaCGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATaa                                                                                      <  1:270922/148‑1 (MQ=255)
                                                          aCGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATaaa                                                                                     <  2:3215545/148‑1 (MQ=255)
                                                          aCGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATaaa                                                                                     <  1:2684736/148‑1 (MQ=255)
                                                           cGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCTTCTGCAGGGCGTTAAGATTAACGAGAAac                                                                                    >  1:2382026/1‑148 (MQ=255)
                                                           cGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAac                                                                                    >  1:1994162/1‑148 (MQ=255)
                                                           cGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAac                                                                                    >  1:1169883/1‑148 (MQ=255)
                                                               tAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACAt                                                                                >  1:2574528/1‑148 (MQ=255)
                                                                  aaCGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGa                                                                             <  1:3472502/148‑1 (MQ=255)
                                                                  aaCGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGa                                                                             <  1:799658/148‑1 (MQ=255)
                                                                     gtggtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGt                                                                          <  2:1798791/148‑1 (MQ=255)
                                                                        gtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAt                                                                       <  1:2778626/148‑1 (MQ=255)
                                                                         tgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtc                                                                      <  1:3016474/148‑1 (MQ=255)
                                                                           aCGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtcgt                                                                    >  2:2541821/1‑148 (MQ=255)
                                                                                aTTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTc                                                               <  1:662069/148‑1 (MQ=255)
                                                                                   tCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGa                                                            >  2:2704998/1‑148 (MQ=255)
                                                                                     cGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGAtg                                                          <  2:3426327/148‑1 (MQ=255)
                                                                                        cGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGAtgctg                                                       <  1:7131/148‑1 (MQ=255)
                                                                                          gTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCg                                                     <  1:3288185/148‑1 (MQ=255)
                                                                                             cGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTaa                                                  <  2:1571465/148‑1 (MQ=255)
                                                                                               gAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAg                                                <  2:3082624/148‑1 (MQ=255)
                                                                                                aaGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGc                                               <  1:1131248/148‑1 (MQ=255)
                                                                                                aaGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGc                                               >  1:55138/1‑148 (MQ=255)
                                                                                                   cgccgcACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTAc                                            >  1:943494/1‑148 (MQ=255)
                                                                                                        cACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCTTTCGTCAGATGCTGCGTAAAGCTACCATCg                                       >  2:2059890/1‑148 (MQ=255)
                                                                                                          cGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGtt                                     <  1:1030615/148‑1 (MQ=255)
                                                                                                           gACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTa                                    <  2:3228091/148‑1 (MQ=255)
                                                                                                           gACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTa                                    <  1:3202048/148‑1 (MQ=255)
                                                                                                           gACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTa                                    <  2:1233143/148‑1 (MQ=255)
                                                                                                            aCATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTaa                                   <  1:2947999/148‑1 (MQ=255)
                                                                                                            aCATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCAGCGTTaa                                   <  1:2220217/148‑1 (MQ=255)
                                                                                                             cATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAAc                                  >  2:408537/1‑148 (MQ=255)
                                                                                                              aTTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAAcg                                 >  2:2942941/1‑148 (MQ=255)
                                                                                                                 ctgcgtctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCgg                              >  2:1346922/1‑148 (MQ=255)
                                                                                                                 ctgcgtctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCgg                              >  2:1347097/1‑148 (MQ=255)
                                                                                                                  tgcgtctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCggg                             >  1:178797/1‑148 (MQ=255)
                                                                                                                   gcgtctgcgtGGTGTTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGt                            <  1:2273709/148‑1 (MQ=255)
                                                                                                                              gtgtTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACtt                 <  2:2341598/148‑1 (MQ=255)
                                                                                                                               tgtTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTc                <  2:2048074/148‑1 (MQ=255)
                                                                                                                               tgtTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTc                <  2:1642564/148‑1 (MQ=255)
                                                                                                                               tgtTCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTc                <  1:2278796/148‑1 (MQ=255)
                                                                                                                                  tCATGCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTg             >  1:2866163/1‑148 (MQ=255)
                                                                                                                                      gCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAg         >  2:315301/1‑148 (MQ=255)
                                                                                                                                      gCTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAg         <  1:1701552/148‑1 (MQ=255)
                                                                                                                                       cTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAgg        >  1:1947363/1‑148 (MQ=255)
                                                                                                                                       cTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAgg        <  2:809064/148‑1 (MQ=255)
                                                                                                                                       cTGTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACACCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAgg        <  2:551688/148‑1 (MQ=255)
                                                                                                                                         gTTACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCg      <  2:733373/148‑1 (MQ=255)
                                                                                                                                          ttACTCTTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGa     >  1:2065102/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATCGTTCGTCAGATGCTGCGTAAAGCTACCATCGTTAACGCGGGTAGCTCCGACTTCCTGGAAGGCGAACA  >  CP009273/4178825‑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: ‑