Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,392,594 A→C 10.4% A248A (GCT→GCG mreC ← cell wall structural complex MreBCD transmembrane component MreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,392,5940AC10.4% 909.5 / 11.3 368A248A (GCT→GCGmreCcell wall structural complex MreBCD transmembrane component MreC
Reads supporting (aligned to +/- strand):  ref base A (187/128);  new base C (3/47);  total (190/178)
Fisher's exact test for biased strand distribution p-value = 1.16e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 376 total aligned reads displayed.

GCATCATCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATCCGTACA  >  CP009273/3392446‑3392743
                                                                                                                                                    |                                                                                                                                                     
gCATCATCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAg                                                                                                                                                      <  2:1105285/150‑1 (MQ=255)
    catcTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAgcgcg                                                                                                                                                  <  1:3504395/150‑1 (MQ=255)
      tcTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCt                                                                                                                                                <  1:2861762/150‑1 (MQ=255)
      tcTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCt                                                                                                                                                <  2:4011928/150‑1 (MQ=255)
         gCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTggg                                                                                                                                             <  2:1070465/150‑1 (MQ=255)
          cATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGt                                                                                                                                            >  2:918030/1‑150 (MQ=255)
          cATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGt                                                                                                                                            <  2:2838802/150‑1 (MQ=255)
            tCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTAt                                                                                                                                          >  2:355350/1‑150 (MQ=255)
            tCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTAt                                                                                                                                          >  1:4081120/1‑150 (MQ=255)
            tCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTAt                                                                                                                                          <  1:3783289/150‑1 (MQ=255)
             cAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATc                                                                                                                                         <  1:3870899/150‑1 (MQ=255)
                aCGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAg                                                                                                                                      <  2:4466208/150‑1 (MQ=255)
                aCGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAg                                                                                                                                      <  1:3129553/150‑1 (MQ=255)
                 cGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGt                                                                                                                                     >  1:2653476/1‑150 (MQ=255)
                  gTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATCAGCGCGCTGGGTATCGAGtt                                                                                                                                    >  2:4148072/1‑150 (MQ=255)
                    tCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGtttt                                                                                                                                  >  1:1838313/1‑150 (MQ=255)
                     cATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGATTCACAGTATAAGCGCGCTGGGTATCGAGTTTTa                                                                                                                                 >  2:2344458/1‑150 (MQ=255)
                     cATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTa                                                                                                                                 <  1:1200277/150‑1 (MQ=255)
                     cATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTa                                                                                                                                 >  1:2575724/1‑150 (MQ=255)
                       ttAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACg                                                                                                                               <  1:1429784/150‑1 (MQ=255)
                                  aTGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAgaga                                                                                                                        <  2:3128182/146‑1 (MQ=255)
                                  aTGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAAc                                                                                                                    >  1:254142/1‑150 (MQ=255)
                                   tGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAAcc                                                                                                                   <  1:198305/150‑1 (MQ=255)
                                      aCCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAACAGACAACcgcg                                                                                                                >  1:3664668/1‑150 (MQ=255)
                                        ctctTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCTGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGAc                                                                                                              >  1:1257940/1‑150 (MQ=255)
                                         tctTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGAcc                                                                                                             >  2:2551842/1‑150 (MQ=255)
                                         tctTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGAcc                                                                                                             <  1:3108258/150‑1 (MQ=255)
                                               ggCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAg                                                                                                       >  2:3237247/1‑150 (MQ=255)
                                                  gTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTCTAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCt                                                                                                    >  2:708183/1‑150 (MQ=255)
                                                   tCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCtt                                                                                                   <  2:2787480/150‑1 (MQ=255)
                                                    cATGGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTc                                                                                                  <  1:214557/150‑1 (MQ=255)
                                                      tCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCg                                                                                                >  2:3735883/1‑150 (MQ=255)
                                                      tCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCg                                                                                                <  1:593400/150‑1 (MQ=255)
                                                       cGGGTTAGCGCCGTTACGATCTGCCCCCCACCGCAGCAGCAGATAACGCAACCGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTCTCAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCgg                                                                                               <  2:1674885/150‑1 (MQ=255)
                                                       cGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCgg                                                                                               >  2:4298771/1‑150 (MQ=255)
                                                         ggTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGa                                                                                             <  2:633578/150‑1 (MQ=255)
                                                            tAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATACGCGCGCTGGGTAGCGATTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAAc                                                                                          <  1:2963238/150‑1 (MQ=255)
                                                            tAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATACGCGCGCGGGGTATCGAGTTTTACGGAATAGACAACCGCGACCGGATAGCCTTCCGGGAAAc                                                                                          <  1:552683/150‑1 (MQ=255)
                                                             aGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACg                                                                                         >  1:3993587/1‑150 (MQ=255)
                                                              gcgcCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGa                                                                                        <  1:1967619/150‑1 (MQ=255)
                                                               cgcCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGAc                                                                                       >  2:201703/1‑150 (MQ=255)
                                                                gcCGTTACGATCGGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGATTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGAcc                                                                                      <  2:243869/150‑1 (MQ=255)
                                                                   gTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGcc                                                                                   >  1:4756257/1‑150 (MQ=255)
                                                                    ttACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGccc                                                                                  >  2:2488836/1‑150 (MQ=255)
                                                                        gATCGGCCCCCCACCGCCGCAGCAGATAACGCAAAGGTTGCAGCCCTGCTGTCGGACGCGCCTGAATCACAGTATACGCGCGGTGGTTTTCGTGTTTTCCGTAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGAc                                                                              <  1:2788894/150‑1 (MQ=255)
                                                                            tGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGa                                                                          >  1:151714/1‑150 (MQ=255)
                                                                             gCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGGTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGaa                                                                         <  2:1418317/150‑1 (MQ=255)
                                                                               ccccccccGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTTAATCACTTTATAAGCGCGCTGGGTATCGTGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGt                                                                       <  2:2970523/150‑1 (MQ=255)
                                                                                    ccagcagcagcatATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTTAACCACCGTATAAGCGCGCGGGTTATCGAGTTTTACGTAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCa                                                                  <  2:488112/149‑1 (MQ=255)
                                                                                     ccgcagcagcatATAACGCAAAGGTTGCAGCCCTGCTGTCGGACGCGCGTTAATCACAGTATAAGCGCGGTGGGTATCGATTTTTAGGTAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAg                                                                 <  1:1924829/148‑1 (MQ=255)
                                                                                      agcagcagcagATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCGGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGc                                                                <  2:2646235/150‑1 (MQ=255)
                                                                                          gcagcagATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACAt                                                            >  1:4594654/1‑150 (MQ=255)
                                                                                           cagcagATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATc                                                           >  2:4539541/1‑150 (MQ=255)
                                                                                           cagcagATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATc                                                           >  2:3252708/1‑150 (MQ=255)
                                                                                            ctcataTACCGAAAAGGTGGCACCCGTGCAGTCGGACGCGCCTTAATAACTTTATACGCGCGGTGGTTATCGTGTTTTAGGGAAGACACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCa                                                          <  1:2597933/145‑1 (MQ=255)
                                                                                            ccgggATAACGCAAACGTGGCAGCCCTGCTGGCGGACGCGCCTGAATCACAGTATACGCGCGGGGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCa                                                          <  2:536194/146‑1 (MQ=255)
                                                                                            agcagATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCa                                                          >  1:4773998/1‑150 (MQ=255)
                                                                                             gcagATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCAc                                                         <  2:2916317/150‑1 (MQ=255)
                                                                                               agATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGACACGACCGCCCAGACCGGAAGTCACCAGCACATCACCa                                                       >  1:1298816/1‑150 (MQ=255)
                                                                                               agATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCa                                                       >  2:1434692/1‑150 (MQ=255)
                                                                                                gATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCaa                                                      >  1:415841/1‑150 (MQ=255)
                                                                                                  tAACGAAAAGTTTGCTGCCGTGCTGTCGGACGCGCTTTAATCAAATTATACGCGGGCTGGTTATCGTTTTTTAGGTAAGATAAAACCGCGACCGTATAGCCTTCCGGTAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAaca                                                    <  1:4680501/150‑1 (MQ=255)
                                                                                                   aaCGCAAAGGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAacac                                                   <  2:1872265/150‑1 (MQ=255)
                                                                                                     cGCAACGGTTGCAGCCCTGCAGTGGGACGCGCCTTAAACACTGTATACGCGCGGGGGGTTTGGTTTTTTCGGGAAGAGAAAACCGCGACCGGATAGCCTTCCGGTAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGa                                                 <  2:1216345/150‑1 (MQ=255)
                                                                                                     cGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGa                                                 >  1:4499824/1‑150 (MQ=255)
                                                                                                     cGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGa                                                 >  1:2999061/1‑150 (MQ=255)
                                                                                                        aaaCGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAata                                              >  1:452466/1‑150 (MQ=255)
                                                                                                         aaCGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAatat                                             >  2:672718/1‑150 (MQ=255)
                                                                                                          aCGTTGCAGCCCTGCAGTGGGACGCGCCTGAATCACTTTTTACGCGCGCGGGTTATCGTTTTTTCCGTAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATc                                            <  1:3915082/150‑1 (MQ=255)
                                                                                                           cGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATcc                                           >  2:2704706/1‑150 (MQ=255)
                                                                                                             ttGCCGCCGTGATGTGGGCCGCGCCTTAAACACATTATACGCGCGTGGGTTTTCGTTTTTTAGGGAAGAAAAAACCGCGACCGGATCGCCTTCGGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGt                                         <  1:546561/150‑1 (MQ=255)
                                                                                                                 aGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTc                                     >  1:1092210/1‑150 (MQ=255)
                                                                                                                           tCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGAt                           >  1:4239838/1‑150 (MQ=255)
                                                                                                                            cGGACGCGCCTTAATCACTTTATACGCGCGGGGGTTATCTTTTTTTACGTAATATACAACCGCGACCGTATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATg                          <  1:3348905/150‑1 (MQ=255)
                                                                                                                            cGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATg                          >  1:693346/1‑150 (MQ=255)
                                                                                                                            cGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATg                          >  2:132287/1‑150 (MQ=255)
                                                                                                                              gACGCGCCTTAATCACAGTATAAGCGCGCTGGGTATCGATTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAGACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCt                        <  1:580017/150‑1 (MQ=255)
                                                                                                                                cgcgCCTTAATCACAGTATACGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCa                      <  2:104168/150‑1 (MQ=255)
                                                                                                                                cgcgCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCa                      >  1:3487260/1‑150 (MQ=255)
                                                                                                                                  cgcCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAg                    >  2:4036626/1‑150 (MQ=255)
                                                                                                                                    ccTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCt                  >  2:4674336/1‑150 (MQ=255)
                                                                                                                                     cTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTg                 >  1:1596412/1‑150 (MQ=255)
                                                                                                                                      tGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGc                >  2:1584317/1‑150 (MQ=255)
                                                                                                                                      tGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCACCGGCAGATGCTCAAGCTGc                >  1:2278571/1‑150 (MQ=255)
                                                                                                                                       gAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCa               >  1:388592/1‑150 (MQ=255)
                                                                                                                                        aaTCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCaa              >  2:2094454/1‑150 (MQ=255)
                                                                                                                                        aaTCAACGTTTAAGCGGGGTGGGTTTGGTGTTTTACTTAAGAGACACCCGCGACCGGATAGCCTTCCGGTAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCaa              <  2:4095751/150‑1 (MQ=255)
                                                                                                                                         aTCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCaaa             >  1:3256349/1‑150 (MQ=255)
                                                                                                                                           cacaGTATATGCGCGCTGGGTATCGAGTTTTACGGATGAGACAACCGCGACCGGATAGCCTTCGGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAatc           <  2:2245090/150‑1 (MQ=255)
                                                                                                                                             caGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAatcat         >  2:50233/1‑150 (MQ=255)
                                                                                                                                              aGTATACGCGCGGTGGGTATCTTGTTTTACGTATGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAatcatc        <  2:152356/150‑1 (MQ=255)
                                                                                                                                              aGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAatcatc        >  1:1761951/1‑150 (MQ=255)
                                                                                                                                               ttataCGCGCGTTGGTTATGGTTTTTTAGGTATGATAAAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATcc       <  1:3813854/149‑1 (MQ=255)
                                                                                                                                               gTATACGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATcc       <  2:3769685/150‑1 (MQ=255)
                                                                                                                                               gTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATcc       >  2:613365/1‑150 (MQ=255)
                                                                                                                                                tataAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATCCg      <  2:3836413/150‑1 (MQ=255)
                                                                                                                                                 ataAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATCCGt     >  1:2483526/1‑150 (MQ=255)
                                                                                                                                                  taAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATCCGTa    >  2:3324455/1‑150 (MQ=255)
                                                                                                                                                  taAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATCCGTa    <  2:3864793/150‑1 (MQ=255)
                                                                                                                                                   aaGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATCCGTAc   >  2:4688712/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCATCATCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATCCGTACA  >  CP009273/3392446‑3392743

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: