Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,596,266 C→A 86.7% T446T (ACG→ACT)  ftsY ← Signal Recognition Particle (SRP) receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,596,2660CA86.7% 418.6 / 67.3 180T446T (ACG→ACT) ftsYSignal Recognition Particle (SRP) receptor
Reads supporting (aligned to +/- strand):  ref base C (10/14);  new base A (76/80);  total (86/94)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 180 total aligned reads displayed.

GGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTCGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTTTCTTCATCACGCGGACGATTTTC  >  CP009273/3596121‑3596410
                                                                                                                                                 |                                                                                                                                                
gggCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCACTACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGt                                                                                                                                                <  2:2709714/148‑1 (MQ=255)
gggCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGt                                                                                                                                                <  2:913372/148‑1 (MQ=255)
gggCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGt                                                                                                                                                <  2:555975/148‑1 (MQ=255)
gggCAAAAAGTGCCTCTAGAACGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGGTCGCCGACACCAATGTAGCGGATAGGGATACCACACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGt                                                                                                                                                <  2:1725370/148‑1 (MQ=255)
 ggCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGtt                                                                                                                                               >  2:3173898/1‑148 (MQ=255)
 ggCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGtt                                                                                                                                               >  1:3149902/1‑148 (MQ=255)
 ggCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGtt                                                                                                                                               <  1:1183168/148‑1 (MQ=255)
 ggCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGtt                                                                                                                                               >  2:1535235/1‑148 (MQ=255)
    aaaaaGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGc                                                                                                                                            >  1:3103379/1‑148 (MQ=255)
      aaaGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGt                                                                                                                                          >  2:44371/1‑148 (MQ=255)
      aaaGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGt                                                                                                                                          >  1:1650584/1‑148 (MQ=255)
        aGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGa                                                                                                                                        <  1:1247681/148‑1 (MQ=255)
         gTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGAt                                                                                                                                       >  2:3143391/1‑148 (MQ=255)
            ccTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTCGTTAGCGTGATGcc                                                                                                                                    >  1:1536983/1‑148 (MQ=255)
              tctATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCgg                                                                                                                                  <  2:1927469/148‑1 (MQ=255)
              tctATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCgg                                                                                                                                  >  1:18877/1‑148 (MQ=255)
              tctATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCgg                                                                                                                                  <  2:259393/148‑1 (MQ=255)
                tataAAGTCGTCCGCCTTCAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGtt                                                                                                                                >  2:559251/1‑148 (MQ=255)
                tataAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGtt                                                                                                                                >  2:525728/1‑148 (MQ=255)
                 ataAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTa                                                                                                                               <  1:2742031/148‑1 (MQ=255)
                 ataAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTa                                                                                                                               <  2:2586812/148‑1 (MQ=255)
                   aaaGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAg                                                                                                                             >  1:539964/1‑148 (MQ=255)
                      gtcgtcCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCa                                                                                                                          <  1:2054331/148‑1 (MQ=255)
                      gtcgtcCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCa                                                                                                                          <  1:2054909/148‑1 (MQ=255)
                       tcgtcCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGGTTAGTTAGCGTGATGCCGGTTAAGCCaa                                                                                                                         >  2:1036487/1‑148 (MQ=255)
                           ccGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGc                                                                                                                     >  2:1458909/1‑148 (MQ=255)
                              ccTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTCGTTAGCGTGATGCCGGTTAAGCCAACGGCTTc                                                                                                                  >  1:1531490/1‑148 (MQ=255)
                               cTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCa                                                                                                                 >  1:1375343/1‑148 (MQ=255)
                                    aCGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGat                                                                                                            <  1:3173898/148‑2 (MQ=255)
                                    aCGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGaa                                                                                                            <  1:1112834/148‑1 (MQ=255)
                                      ggACGCAAATCCTCAATACGTTCGCCGACACCACTGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTCGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACa                                                                                                          <  1:2148675/148‑1 (MQ=255)
                                      ggACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTCGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACa                                                                                                          <  1:2149002/148‑1 (MQ=255)
                                       gACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAg                                                                                                         >  2:608777/1‑148 (MQ=255)
                                        aCGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGt                                                                                                        <  1:2927398/148‑1 (MQ=255)
                                        aCGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGt                                                                                                        <  1:3029741/148‑1 (MQ=255)
                                             aaTCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTCGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGc                                                                                                   <  1:155270/148‑1 (MQ=255)
                                             aaTCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGc                                                                                                   >  2:2715379/1‑148 (MQ=255)
                                             aaTCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGc                                                                                                   >  2:723471/1‑148 (MQ=255)
                                             aaTCCTCAATACGTTCGCCGACACCAATGGAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGc                                                                                                   <  2:1551229/148‑1 (MQ=255)
                                                  tCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGc                                                                                              <  2:135994/148‑1 (MQ=255)
                                                  tCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGc                                                                                              <  2:18877/148‑1 (MQ=255)
                                                       aCGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTAcc                                                                                         >  2:486578/1‑148 (MQ=255)
                                                       aCGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTAcc                                                                                         >  1:1648508/1‑148 (MQ=255)
                                                           tCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGt                                                                                     >  2:3094607/1‑148 (MQ=255)
                                                                gACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGc                                                                                <  1:2669883/148‑1 (MQ=255)
                                                                 acacCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGcc                                                                               >  1:1673053/1‑148 (MQ=255)
                                                                 acacCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGcc                                                                               >  1:1404989/1‑148 (MQ=255)
                                                                  cacCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGccc                                                                              >  2:2552645/1‑148 (MQ=255)
                                                                   acCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCg                                                                             >  1:1042272/1‑148 (MQ=255)
                                                                      aaTGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTg                                                                          <  1:384248/148‑1 (MQ=255)
                                                                           aGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTCGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGc                                                                     >  2:2015198/1‑148 (MQ=255)
                                                                            gCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCa                                                                    >  2:1801244/1‑148 (MQ=255)
                                                                            gCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCa                                                                    >  2:1801208/1‑148 (MQ=255)
                                                                               gATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACACCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCa                                                                 >  2:1684230/1‑148 (MQ=255)
                                                                                  aGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGTCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATa                                                              <  1:72813/148‑1 (MQ=255)
                                                                                  aGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGTCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATa                                                              <  1:72832/148‑1 (MQ=255)
                                                                                  aGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATa                                                              >  2:594357/1‑148 (MQ=255)
                                                                                  aGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATa                                                              <  1:2905605/148‑1 (MQ=255)
                                                                                  aGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATa                                                              <  1:2905622/148‑1 (MQ=255)
                                                                                    ggATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGt                                                            <  2:1635132/148‑1 (MQ=255)
                                                                                     gATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTc                                                           >  2:989396/1‑148 (MQ=255)
                                                                                      aTACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCa                                                          <  2:767843/148‑1 (MQ=255)
                                                                                         ccAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCa                                                       <  2:807819/148‑1 (MQ=255)
                                                                                          cAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTCGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCAt                                                      >  2:2929011/1‑148 (MQ=255)
                                                                                          cAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCAt                                                      <  2:1093759/148‑1 (MQ=255)
                                                                                             aCTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAAc                                                   >  2:577646/1‑148 (MQ=255)
                                                                                             aCTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAAc                                                   <  1:3001538/148‑1 (MQ=255)
                                                                                               tGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTCGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACtt                                                 <  1:2878638/148‑1 (MQ=255)
                                                                                               tGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACtt                                                 >  2:1291854/1‑148 (MQ=255)
                                                                                                 gTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCa                                               >  1:3071615/1‑148 (MQ=255)
                                                                                                  tCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCAt                                              >  2:2337383/1‑148 (MQ=255)
                                                                                                    aGCCACCGAGAAAATTACCCCGCCTGTCGCCGTGCCGTCCAGTTTCGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATgc                                            <  1:2290942/148‑1 (MQ=255)
                                                                                                      ccACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATgcgg                                          >  2:3295592/1‑148 (MQ=255)
                                                                                                      ccACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATgcgg                                          >  1:1152785/1‑148 (MQ=255)
                                                                                                           gagaAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTCGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGTGGCGCtt                                     <  1:2398245/148‑1 (MQ=255)
                                                                                                             gaAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCa                                   >  1:1071135/1‑148 (MQ=255)
                                                                                                              aaaaTTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCaa                                  <  1:2834768/148‑1 (MQ=255)
                                                                                                              aaaaTTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCaa                                  >  2:764287/1‑148 (MQ=255)
                                                                                                              aaaaTTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCaa                                  >  1:2305250/1‑148 (MQ=255)
                                                                                                                 aTTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGt                               >  1:1960257/1‑148 (MQ=255)
                                                                                                                 aTTACCCCGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGt                               >  1:1960939/1‑148 (MQ=255)
                                                                                                                     ccccGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATGGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAg                           >  1:2222404/1‑148 (MQ=255)
                                                                                                                       ccGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGtt                         >  1:1459217/1‑148 (MQ=255)
                                                                                                                       ccGCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAAATTCATGCGGCGCTTCAACGTCGAGtt                         >  1:1023883/1‑148 (MQ=255)
                                                                                                                         gCCTTTCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTttc                       <  2:479001/148‑1 (MQ=255)
                                                                                                                            tttCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTttcttc                    <  2:1209989/148‑1 (MQ=255)
                                                                                                                            tttCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTttcttc                    >  2:2431111/1‑148 (MQ=255)
                                                                                                                            tttCGCCGTGCCGTAGAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTttcttc                    <  1:3226979/148‑1 (MQ=255)
                                                                                                                             ttCGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTTTCTtca                   >  1:1183894/1‑148 (MQ=255)
                                                                                                                               cGCCGTGCCGTCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTTTCTtcatc                 <  2:2247304/148‑1 (MQ=255)
                                                                                                                                    tgccgtCCAGTTTAGTTAGCGTGATGCCGCTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTTTCTTCATCAcgcg            <  2:1666671/148‑1 (MQ=255)
                                                                                                                                     gccgtCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTTTCTTCATCACGCgg           <  2:2752749/148‑1 (MQ=255)
                                                                                                                                     gccgtCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTTTCTTCATCACGCgg           <  1:2992600/148‑1 (MQ=255)
                                                                                                                                      ccgtCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTTTCTTCATCACGCGGa          <  1:976355/148‑1 (MQ=255)
                                                                                                                                        gtCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTTTCTTCATCACGCGGACg        <  2:580697/148‑1 (MQ=255)
                                                                                                                                         tCCAGTTTAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTTTCTTCATCACGCGGACGa       >  1:695221/1‑148 (MQ=255)
                                                                                                                                          ccAGTTTCGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTTTCTTCATCACGCGGACGAt      >  1:3290509/1‑148 (MQ=255)
                                                                                                                                              tttAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTTTCTTCATCACGCGGACGATTttc  <  1:562232/148‑1 (MQ=255)
                                                                                                                                              tttAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTTTCTTCATCACGCGGACGATTttc  <  1:180555/148‑1 (MQ=255)
                                                                                                                                              tttAGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTTTCTTCATCACGCGGACGATTttc  <  1:1091596/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTAGCGGATAGGGATACCAAACTGGTCAGCCACCGAGAAAATTACCCCGCCTTTCGCCGTGCCGTCCAGTTTCGTTAGCGTGATGCCGGTTAAGCCAACGGCTTCATGGAACAGTTTGGCCTGGCTTACCGCGTTCTGCCCGGTGCTGGCATCAATAGTCAGCATAACTTCATGCGGCGCTTCAACGTCGAGTTTCTTCATCACGCGGACGATTTTC  >  CP009273/3596121‑3596410

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: ‑