Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,446,567 T→A 6.5% A21A (GCA→GCT rpsJ ← 30S ribosomal subunit protein S10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,446,5670TA 6.5% 868.0 / 37.7 276A21A (GCA→GCTrpsJ30S ribosomal subunit protein S10
Reads supporting (aligned to +/- strand):  ref base T (95/163);  new base A (9/9);  total (104/172)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 279 total aligned reads displayed.

GTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCT  >  CP009273/3446420‑3446712
                                                                                                                                                   |                                                                                                                                                 
gTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGtt                                                                                                                                                   >  1:1599849/1‑148 (MQ=255)
 tGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTg                                                                                                                                                  >  2:632945/1‑148 (MQ=255)
 tGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTg                                                                                                                                                  >  1:2159698/1‑148 (MQ=255)
  gAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGc                                                                                                                                                 <  1:2803148/148‑1 (MQ=255)
     tACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTg                                                                                                                                              >  2:3489742/1‑148 (MQ=255)
      aCGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTga                                                                                                                                             <  2:425201/148‑1 (MQ=255)
      aCGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTga                                                                                                                                             <  2:2987782/148‑1 (MQ=255)
      aCGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTga                                                                                                                                             <  2:353725/148‑1 (MQ=255)
       cGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTgat                                                                                                                                            >  1:1799991/1‑148 (MQ=255)
       cGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGAACACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTgat                                                                                                                                            >  2:59082/1‑148 (MQ=255)
       cGGAGTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTgat                                                                                                                                            <  2:1016702/148‑1 (MQ=255)
          aTTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTgatcga                                                                                                                                         >  1:1228422/1‑148 (MQ=255)
            ttCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTgatcgatc                                                                                                                                       >  2:1043574/1‑148 (MQ=255)
              cGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAg                                                                                                                                     >  1:1437901/1‑148 (MQ=255)
                 aCTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACg                                                                                                                                  <  1:1635452/148‑1 (MQ=255)
                 aCTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACg                                                                                                                                  <  1:650630/148‑1 (MQ=255)
                   tGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACgat                                                                                                                                >  1:731838/1‑148 (MQ=255)
                   tGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACgat                                                                                                                                >  1:204130/1‑148 (MQ=255)
                      tCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACgatgat                                                                                                                             <  1:911434/148‑1 (MQ=255)
                      tCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACgatgat                                                                                                                             >  1:3062249/1‑148 (MQ=255)
                         cgcgcgTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCaa                                                                                                                          <  2:1622258/148‑1 (MQ=255)
                              gTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCt                                                                                                                     <  1:2376620/148‑1 (MQ=255)
                               tCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCtt                                                                                                                    >  2:2396875/1‑148 (MQ=255)
                                  ttgttgACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCa                                                                                                                 >  1:248927/1‑148 (MQ=255)
                                    gttgACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAgg                                                                                                               <  2:2549094/148‑1 (MQ=255)
                                      tgACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCg                                                                                                             >  1:2959279/1‑148 (MQ=255)
                                       gACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCgg                                                                                                            <  2:708361/148‑1 (MQ=255)
                                        aCGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGa                                                                                                           <  2:2930353/148‑1 (MQ=255)
                                        aCGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGAAGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGa                                                                                                           <  2:2930224/148‑1 (MQ=255)
                                         cGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGAt                                                                                                          >  1:1003092/1‑148 (MQ=255)
                                             cGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACgg                                                                                                      <  1:3779080/148‑1 (MQ=255)
                                               gggAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGAt                                                                                                    >  2:2681432/1‑148 (MQ=255)
                                               gggAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGAt                                                                                                    <  2:555330/148‑1 (MQ=255)
                                                  agaTCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTAGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCt                                                                                                 >  2:3684749/1‑148 (MQ=255)
                                                    aTCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCttt                                                                                               >  2:1382392/1‑148 (MQ=255)
                                                     tCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTg                                                                                              <  1:2587010/148‑1 (MQ=255)
                                                     tCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTg                                                                                              >  2:1557583/1‑148 (MQ=255)
                                                      cAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTgg                                                                                             >  1:2132674/1‑148 (MQ=255)
                                                       aGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGt                                                                                            <  1:2190621/148‑1 (MQ=255)
                                                       aGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGt                                                                                            <  1:3574518/148‑1 (MQ=255)
                                                        gAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGtt                                                                                           >  1:423535/1‑148 (MQ=255)
                                                        gAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGtt                                                                                           >  1:3340946/1‑148 (MQ=255)
                                                             gagAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCa                                                                                      <  2:3655693/146‑1 (MQ=255)
                                                                  gcgcTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATgaga                                                                                 >  1:670783/1‑148 (MQ=255)
                                                                   cgcTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGAc                                                                                <  2:1247871/148‑1 (MQ=255)
                                                                    gcTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGAcc                                                                               <  1:2596916/148‑1 (MQ=255)
                                                                      tctTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCag                                                                             <  1:3011035/148‑1 (MQ=255)
                                                                      tctTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTAGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCag                                                                             <  2:3059670/148‑1 (MQ=255)
                                                                       ctTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCaga                                                                            >  1:1983904/1‑148 (MQ=255)
                                                                        tttGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCagag                                                                           <  1:32868/148‑1 (MQ=255)
                                                                         ttGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGc                                                                          <  1:1308355/148‑1 (MQ=255)
                                                                         ttGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGc                                                                          >  2:1471062/1‑148 (MQ=255)
                                                                          tGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCt                                                                         <  1:2688802/148‑1 (MQ=255)
                                                                           gCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTc                                                                        <  1:1868464/148‑1 (MQ=255)
                                                                               gtCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAat                                                                    <  1:2387321/148‑1 (MQ=255)
                                                                               gtCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAat                                                                    <  1:3126534/148‑1 (MQ=255)
                                                                                tCGGCAGCGGGATCGGACCACGGAACTGCGCACCAGTGCGCTTTGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGAGCAGAGCTCCAatt                                                                   <  1:1363029/148‑1 (MQ=255)
                                                                                 aggaatcGGGATCGGACCACGGACCTGCGCAACAGTGAGCTTGGCAGTCTCGACTATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACTCTTTCAGGAGGATACTGATTCTTTGGTTCTGCATGATACCAGAGCTCCAatta                                                                  >  2:2452314/7‑148 (MQ=255)
                                                                                  ggCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAattat                                                                 <  1:2401171/148‑1 (MQ=255)
                                                                                   gCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAattatt                                                                <  2:2545169/148‑1 (MQ=255)
                                                                                     aGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTAtttt                                                              <  1:529738/148‑1 (MQ=255)
                                                                                      gCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTta                                                             <  2:3703228/148‑1 (MQ=255)
                                                                                      gCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTAGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTta                                                             <  1:1341686/148‑1 (MQ=255)
                                                                                         ggATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATaa                                                          <  1:1865193/148‑1 (MQ=255)
                                                                                          gATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATaaa                                                         <  2:2352643/148‑1 (MQ=255)
                                                                                                aCCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGaaaa                                                   >  1:1231748/1‑148 (MQ=255)
                                                                                                aCCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGaaaa                                                   <  1:156774/148‑1 (MQ=255)
                                                                                                  cACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATg                                                 >  2:935214/1‑148 (MQ=255)
                                                                                                  cACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATg                                                 <  1:3489742/148‑1 (MQ=255)
                                                                                                  cACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATg                                                 >  1:371606/1‑148 (MQ=255)
                                                                                                   aCGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGa                                                <  2:1799991/148‑1 (MQ=255)
                                                                                                   aCGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGa                                                <  2:2794500/148‑1 (MQ=255)
                                                                                                      gACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTa                                             <  2:3220569/148‑1 (MQ=255)
                                                                                                      gACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTa                                             <  1:1765048/148‑1 (MQ=255)
                                                                                                         cTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTActc                                          >  2:691984/1‑148 (MQ=255)
                                                                                                          tGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTActcc                                         >  1:1137574/1‑148 (MQ=255)
                                                                                                           gcgcACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTActcct                                        >  1:2295959/1‑148 (MQ=255)
                                                                                                              cACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAg                                     >  1:3315838/1‑148 (MQ=255)
                                                                                                               aCCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGa                                    >  2:1085347/1‑148 (MQ=255)
                                                                                                                   gTGCTCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCa                                <  2:1695240/148‑1 (MQ=255)
                                                                                                                    tGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCAt                               <  2:1993151/148‑1 (MQ=255)
                                                                                                                     gcgcTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCCGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCAtt                              >  1:1534955/1‑148 (MQ=255)
                                                                                                                        cTTGGCTGTATCGACTTTTTCCTCGGTTTCTTGATCTATCCGCCGATGATCTAACTCTTTCAGGCGGATACGGTTTCTTTGGTTCTGCATGAGTCCAGAGCTCCAATTATGTTATTAACGCAAATGATTACTCCTCAGACCCATTACg                           <  1:2901474/148‑1 (MQ=255)
                                                                                                                         ttGGCAGTCTCGACGATTTCCGCGGTAGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGa                          <  2:1779694/148‑1 (MQ=255)
                                                                                                                           ggCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGAtt                        <  2:774878/148‑1 (MQ=255)
                                                                                                                            gCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTg                       <  2:3739292/148‑1 (MQ=255)
                                                                                                                             cAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGa                      <  2:2082031/148‑1 (MQ=255)
                                                                                                                             cAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGa                      <  1:318840/148‑1 (MQ=255)
                                                                                                                             cAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGa                      >  1:829651/1‑148 (MQ=255)
                                                                                                                                 ctcGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATggg                  >  1:3211812/1‑148 (MQ=255)
                                                                                                                                 ctcGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATggg                  <  1:599119/148‑1 (MQ=255)
                                                                                                                                   cgacgaTTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGgag                >  2:3374931/1‑148 (MQ=255)
                                                                                                                                    gacgaTTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGgaga               <  1:1944684/148‑1 (MQ=255)
                                                                                                                                     acgaTTTCCGCGGTAGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGagag              >  2:755370/1‑148 (MQ=255)
                                                                                                                                       gaTTTCCGCGGTAGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAgtg            <  1:3236255/148‑1 (MQ=255)
                                                                                                                                         tttCCGCGGTAGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTa          >  1:2511305/1‑148 (MQ=255)
                                                                                                                                              gcgGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGt     <  1:256492/148‑1 (MQ=255)
                                                                                                                                                ggTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTc   <  1:632882/148‑1 (MQ=255)
                                                                                                                                                 gTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCt  >  2:400684/1‑148 (MQ=255)
                                                                                                                                                 gTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCt  <  2:2646860/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTGAGTACGGATTTCGTACTGATCGCGCGCGTCTTTGTTGACGTGCGGGGAGATCAGAACAGTGAAGCGCTCTTTGCGTGTCGGCAGCGGGATCGGACCACGGACCTGCGCACCAGTGCGCTTGGCAGTCTCGACGATTTCCGCGGTTGCTTGATCGATCAGACGATGATCAAACGCTTTCAGGCGGATACGGATTCTTTGGTTCTGCATGAGACCAGAGCTCCAATTATTTTATAAACGAAAATGATTACTCCTCAGACCCATTACGATTGATGGGAGAGTGTAACCGTTCT  >  CP009273/3446420‑3446712

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: