Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,392,594 A→C 8.7% 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,5940AC 8.7% 690.5 / 4.1 301A248A (GCT→GCG) mreCcell wall structural complex MreBCD transmembrane component MreC
Reads supporting (aligned to +/- strand):  ref base A (147/113);  new base C (4/32);  total (151/150)
Fisher's exact test for biased strand distribution p-value = 1.58e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 306 total aligned reads displayed.

GGCATCATCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATCCGTACA  >  CP009273/3392445‑3392743
                                                                                                                                                     |                                                                                                                                                     
ggCATCATCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATaa                                                                                                                                                       <  1:5276391/150‑1 (MQ=255)
  catcatCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAgc                                                                                                                                                     <  1:3799312/150‑1 (MQ=255)
  catcatCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAgc                                                                                                                                                     >  1:4017581/1‑150 (MQ=255)
   atcatcTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAgcg                                                                                                                                                    >  1:1928827/1‑150 (MQ=255)
    tcatcTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAgcgc                                                                                                                                                   <  1:2997784/150‑1 (MQ=255)
      atcTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAgcgcgc                                                                                                                                                 >  1:2705635/1‑150 (MQ=255)
      atcTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAgcgcgc                                                                                                                                                 >  2:2665090/1‑150 (MQ=255)
      atcTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAgcgcgc                                                                                                                                                 >  2:3046383/1‑150 (MQ=255)
          gCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATCAGCGCGCTggg                                                                                                                                             >  2:243239/1‑150 (MQ=255)
          gCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTggg                                                                                                                                             >  1:884330/1‑150 (MQ=255)
           cATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGt                                                                                                                                            <  1:337694/150‑1 (MQ=255)
            aTCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTa                                                                                                                                           >  2:381563/1‑150 (MQ=255)
              cAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATc                                                                                                                                         <  1:3292136/150‑1 (MQ=255)
               aGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCg                                                                                                                                        <  2:3642435/150‑1 (MQ=255)
                gACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGa                                                                                                                                       <  2:1535560/150‑1 (MQ=255)
                   gTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGtt                                                                                                                                    >  1:2644494/1‑150 (MQ=255)
                    ttCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGttt                                                                                                                                   <  1:3741643/150‑1 (MQ=255)
                    ttCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGttt                                                                                                                                   <  1:4258646/150‑1 (MQ=255)
                     tCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGtttt                                                                                                                                  <  1:1060042/150‑1 (MQ=255)
                     tCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGtttt                                                                                                                                  <  1:2922104/150‑1 (MQ=255)
                        ttAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACg                                                                                                                               <  2:3491754/150‑1 (MQ=255)
                          aGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGa                                                                                                                             >  2:3177493/1‑150 (MQ=255)
                          aGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGa                                                                                                                             <  1:468890/150‑1 (MQ=255)
                          aGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGa                                                                                                                             >  1:4998722/1‑150 (MQ=255)
                            cAACACGATGCACCTCGTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAag                                                                                                                           <  2:3692564/150‑1 (MQ=255)
                             aaCACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAaga                                                                                                                          <  1:3168108/150‑1 (MQ=255)
                                acGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGAc                                                                                                                       <  2:2798845/150‑1 (MQ=255)
                                   aTGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAAc                                                                                                                    <  2:4387837/150‑1 (MQ=255)
                                   aTGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAAc                                                                                                                    <  2:900971/150‑1 (MQ=255)
                                     gCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACcg                                                                                                                  >  2:2265347/1‑150 (MQ=255)
                                         ctctTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCACACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATACGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGAc                                                                                                              >  2:3370332/1‑150 (MQ=255)
                                              ccGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGGTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCggtt                                                                                                         <  2:3257599/150‑3 (MQ=255)
                                                      aTCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTcc                                                                                                 >  2:4029375/1‑150 (MQ=255)
                                                      aTCGGGTTAGCGCCGTTACAATCTGCCCCCCACCGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTcc                                                                                                 <  1:4166860/150‑1 (MQ=255)
                                                           gTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCGGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGaa                                                                                            <  2:1018871/150‑1 (MQ=255)
                                                             tAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTCTCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAAc                                                                                          >  2:1933546/1‑150 (MQ=255)
                                                               gcgcCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGACGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGa                                                                                        <  1:1687715/150‑1 (MQ=255)
                                                                   cGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGc                                                                                    >  2:1447243/1‑150 (MQ=255)
                                                                     ttACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTTAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGccc                                                                                  <  2:286503/150‑1 (MQ=255)
                                                                         gATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTAGACGCGCGCTGGGTCTCGAGTTTTACGGAAGAGACAACCGCGACCGGAGAGCCTTCCGGGAACCGACCGCCCAGAc                                                                              >  2:4408665/1‑150 (MQ=255)
                                                                          aTCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGAcc                                                                             >  2:4273297/1‑150 (MQ=255)
                                                                           tCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCg                                                                            <  2:3437786/150‑1 (MQ=255)
                                                                            cTGCCCCCCACAGCCGCCGCAGATAACGCAACCGTTGCCGCCCTGCAGTCGGACGCGCCTGAAGCACCGTATACCCGCGCGGGGTAGCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCgg                                                                           <  2:2838069/150‑1 (MQ=255)
                                                                             tGCCCCCCACAGCAGCTGCATATAACGCAAAGGTTGCAGCCCGGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCGGGGTTTCGTGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGa                                                                          <  2:5104419/150‑1 (MQ=255)
                                                                             tGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACTGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGa                                                                          <  2:1266229/150‑1 (MQ=255)
                                                                             gcccccccACCGCCGCAGCAGATACCGAAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAAACACCGTATACGCGCGCGGGGTATCGTGTTTTACGTAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGa                                                                          <  2:2458329/148‑1 (MQ=255)
                                                                               ccccccACCGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATACGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAg                                                                        <  1:4883131/150‑1 (MQ=255)
                                                                                 ccccACCGCAGCAGCAGATAACGCAACCGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTc                                                                      <  2:4637575/150‑1 (MQ=255)
                                                                                 caccACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACAGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTc                                                                      >  2:4964304/3‑150 (MQ=255)
                                                                                       cgcagcagaatATACCGCAACCGTTGCAGCCCTGCAGTCGGACGCCCCTGAATCACAGTATCAGCGCGCTGGTTATCGTGTTTTACGGACTAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGc                                                                <  2:4517461/149‑1 (MQ=255)
                                                                                       agcagcagcagATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGc                                                                >  1:343345/1‑150 (MQ=255)
                                                                                         cagcagcagATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGcac                                                              >  1:3515939/1‑150 (MQ=255)
                                                                                         cagcagcagATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGcac                                                              >  1:4335072/1‑150 (MQ=255)
                                                                                          agcagcagATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGcaca                                                             >  1:1342810/1‑150 (MQ=255)
                                                                                          agcagcagATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGcaca                                                             >  1:694508/1‑150 (MQ=255)
                                                                                              gcagATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCACTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCAc                                                         >  2:778638/1‑150 (MQ=255)
                                                                                               cagATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCAcc                                                        >  2:2070762/1‑150 (MQ=255)
                                                                                               cagATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCAcc                                                        <  2:3912321/150‑1 (MQ=255)
                                                                                                  aTAACGCAACCGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAac                                                     <  2:702856/150‑1 (MQ=255)
                                                                                                  aTAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAATGACCGCCCAGACCGGAAGTCACCAGCACATCACCAac                                                     >  2:4958696/1‑150 (MQ=255)
                                                                                                   tAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAaca                                                    >  1:1746527/1‑150 (MQ=255)
                                                                                                        cAAAGGTTGCCGCCCTGTTGTCGGACGCGCCTTAATAACATTATAAGCGGGTGGGTTTTCGAGTTTTTCGTAAGAGAAAACCGCGACCGGTTAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAat                                               <  1:1255482/150‑1 (MQ=255)
                                                                                                        cAAAGGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAat                                               <  2:1143365/150‑1 (MQ=255)
                                                                                                          ctgggtGCCGCCCTGCTGTCGTACGCGCCTTAATCACTTTTTAGGGGGGGTGGTTATCGTTTTTTAGGTAAGATAAAACCGCGACCGGATAGCCTTCCGGTAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAatat                                             <  1:3218765/145‑1 (MQ=255)
                                                                                                          caggTTGCAGCCCTGATGTCGGACGCGCCTGATTCACAGTATACGCGCGGTGGGTATGGTGTTTTAGGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAatat                                             <  2:2474402/147‑1 (MQ=255)
                                                                                                          aaCGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAatat                                             <  2:1612841/150‑1 (MQ=255)
                                                                                                           ccGTGGCAGCCCTGCTGTGGGACGCGCCTTAATCACCGTATAAGCGCGCGGGGTATCGTGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATc                                            <  2:4573860/149‑1 (MQ=255)
                                                                                                           aCGTTGAAGCCCTGCTGTCGGACGCGCCTGAATAACTGTATAAGCGCGGTGGGTATCGAGTTTTAGGTAAGATACAACCGCGACGGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATc                                            <  2:4099665/150‑1 (MQ=255)
                                                                                                            cGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGTAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATcc                                           <  2:4805793/150‑1 (MQ=255)
                                                                                                            cGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATcc                                           >  2:5314428/1‑150 (MQ=255)
                                                                                                             gTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCg                                          >  1:142595/1‑150 (MQ=255)
                                                                                                              ttGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGATACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGt                                         <  2:5226431/150‑1 (MQ=255)
                                                                                                                gCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTAt                                       >  1:1189230/1‑150 (MQ=255)
                                                                                                                 cAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTAtt                                      >  1:2045686/1‑150 (MQ=255)
                                                                                                                 cAGCCCTGATGTCGGACGCGCCTTAATCACTTTATACGCGGGCGGGGTATCGAGTTTTACGTACGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTAtt                                      <  2:1849589/150‑1 (MQ=255)
                                                                                                                  aGCCCTGCAGTCGGACGCGCCTGACTCACAGTATCAGCGCGCTGGGTATCGAGTTTAACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCCCCAGCACCTCACCAACACGAATATCCGTATTc                                     <  2:1755042/150‑1 (MQ=255)
                                                                                                                   gCCCTGCTGTCGGACGCGCCTGAAACACATTATATGCGCGCTGGGTATCGCGTTTTACGTATGAGACAACCGCGCCCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCg                                    <  2:1985452/150‑1 (MQ=255)
                                                                                                                   gCCCTGCAGTCGGACGCGCCTGAATCACAGTATACGCGCGCTGGGTATCGCGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCg                                    <  2:1670462/150‑1 (MQ=255)
                                                                                                                    cccTGCTGTGGGACGCGCCTTAAACACTGTCTACGCGCGCTGGGTATCTATTTTTACGTAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGc                                   <  2:4926096/150‑1 (MQ=255)
                                                                                                                    cccTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGc                                   >  1:5167853/1‑150 (MQ=255)
                                                                                                                       tGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCgg                                >  2:415030/1‑150 (MQ=255)
                                                                                                                        gCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGACGAGCCGTATTCGCCgtc                               >  2:3305035/1‑148 (MQ=255)
                                                                                                                          aGTAGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAg                             >  1:387246/1‑150 (MQ=255)
                                                                                                                              ggACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGc                         >  1:243469/1‑150 (MQ=255)
                                                                                                                                aCGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTCCGGAAGAGACACCCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTc                       <  1:4392293/150‑1 (MQ=255)
                                                                                                                                   cgcCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAg                    >  2:2280139/1‑150 (MQ=255)
                                                                                                                                     ccTGAATCACAGTATACGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCt                  <  2:2705635/150‑1 (MQ=255)
                                                                                                                                      gttaaTCACATTATCCGCGCGTTGGTTATCTATTTTTACGCAATACAAACCCGCGACCGTATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTg                 <  1:1877095/147‑1 (MQ=255)
                                                                                                                                      cTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTg                 >  1:1658660/1‑150 (MQ=255)
                                                                                                                                        taaTCACTTTATACGCGCGCTGGTTTTCGTGTTTTACGTAAGAGACAACCGCGACCGTATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCa               <  1:2143738/149‑1 (MQ=255)
                                                                                                                                         aaTCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCTGTATTCGCCGGCAGATGCTCAAGCTGCaa              >  1:577376/1‑150 (MQ=255)
                                                                                                                                         aaTCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCac              >  2:1205287/1‑149 (MQ=255)
                                                                                                                                         aaTCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCaa              >  2:4807837/1‑150 (MQ=255)
                                                                                                                                          aTCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCaaa             >  2:1225114/1‑150 (MQ=255)
                                                                                                                                           tCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAat            >  1:3652926/1‑150 (MQ=255)
                                                                                                                                            cacaGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAatc           >  1:354365/1‑150 (MQ=255)
                                                                                                                                                gTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATcc       >  1:4426701/1‑150 (MQ=255)
                                                                                                                                                  ataCGCGCGCTGGTTATCGATTTTTACGTAAGACACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATCCGt     <  1:1643020/150‑1 (MQ=255)
                                                                                                                                                  ataAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATCCGt     <  2:644006/150‑1 (MQ=255)
                                                                                                                                                   taAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATCCGTa    >  1:2374286/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCATCATCTGCATCAGACGTTCATTAGCAACACGATGCACCTCTTCCGGCGTCATCGGGTTAGCGCCGTTACGATCTGCCCCCCACAGCAGCAGCAGATAACGCAAACGTTGCAGCCCTGCAGTCGGACGCGCCTGAATCACAGTATAAGCGCGCTGGGTATCGAGTTTTACGGAAGAGACAACCGCGACCGGATAGCCTTCCGGGAAACGACCGCCCAGACCGGAAGTCACCAGCACATCACCAACACGAATATCCGTATTCGCCGGCAGATGCTCAAGCTGCAAATCATCCGTACA  >  CP009273/3392445‑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: ‑