Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,617,102 G→T 20.2% D234Y (GAT→TAT)  ppk → polyphosphate kinase, component of RNA degradosome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,617,1020GT20.2% 633.5 / 122.9 231D234Y (GAT→TAT) ppkpolyphosphate kinase, component of RNA degradosome
Reads supporting (aligned to +/- strand):  ref base G (98/86);  new base T (29/18);  total (127/104)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

ATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCC  >  CP009273/2616953‑2617251
                                                                                                                                                     |                                                                                                                                                     
aTTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACg                                                                                                                                                       <  2:2665839/150‑1 (MQ=255)
  ttACCGCCAGTAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGAt                                                                                                                                                     <  1:4236733/150‑1 (MQ=255)
   tACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGAtt                                                                                                                                                    <  1:3147478/150‑1 (MQ=255)
   tACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGAtt                                                                                                                                                    <  2:615381/150‑1 (MQ=255)
   tACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGAtt                                                                                                                                                    <  1:2122121/150‑1 (MQ=255)
       gCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGt                                                                                                                                                >  1:4719866/1‑150 (MQ=255)
       gCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGt                                                                                                                                                >  1:1942695/1‑150 (MQ=255)
        ccAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTg                                                                                                                                               >  1:2959950/1‑150 (MQ=255)
          aGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCa                                                                                                                                             >  2:553001/1‑150 (MQ=255)
               cgccgcGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATgaga                                                                                                                                        <  2:849940/150‑1 (MQ=255)
                 ccgcGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGTGCATGAGATg                                                                                                                                      <  1:4077414/150‑1 (MQ=255)
                 ccgcGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGTGCATGAGATg                                                                                                                                      >  1:3710795/1‑150 (MQ=255)
                    cgTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGaa                                                                                                                                   >  1:2763113/1‑150 (MQ=255)
                      tCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTACCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAagc                                                                                                                                 <  1:766832/150‑1 (MQ=255)
                           gCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAagccagcc                                                                                                                            >  2:4835644/1‑150 (MQ=255)
                             aaGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTg                                                                                                                          >  1:2526889/1‑150 (MQ=255)
                             aaGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTg                                                                                                                          <  1:1164788/150‑1 (MQ=255)
                               gCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGAt                                                                                                                        >  1:4776289/1‑150 (MQ=255)
                                ccGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATg                                                                                                                       <  1:574548/150‑1 (MQ=255)
                                ccGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATg                                                                                                                       <  1:3006935/150‑1 (MQ=255)
                                 cGATGATTCTTCTGGATAACATTCTGCGTTAATGCCTTGATGATTTTTTCATAGGCTACTTTTATTATGACGCGCTGAATGCCTATTGAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGAATGAGATGGAAGCCAGCCTGTTgg                                                                                                                      <  2:2800168/150‑1 (MQ=255)
                                        tcttctGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGAt                                                                                                               <  1:2167822/150‑1 (MQ=255)
                                            ctGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCt                                                                                                           <  1:800186/150‑1 (MQ=255)
                                            ctGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCt                                                                                                           <  2:148198/150‑1 (MQ=255)
                                              ggATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTc                                                                                                         >  2:1750639/1‑150 (MQ=255)
                                                aTAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCa                                                                                                       >  2:875248/1‑150 (MQ=255)
                                                aTAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCa                                                                                                       <  1:3309188/150‑1 (MQ=255)
                                                 tAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAg                                                                                                      >  2:4096410/1‑150 (MQ=255)
                                                     aTTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGtctc                                                                                                  >  1:2092494/1‑150 (MQ=255)
                                                       tCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCaa                                                                                                <  2:3989316/150‑1 (MQ=255)
                                                        cTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAag                                                                                               <  1:340823/150‑1 (MQ=255)
                                                        cTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAag                                                                                               <  2:2134910/150‑1 (MQ=255)
                                                        cTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAag                                                                                               <  2:3486185/150‑1 (MQ=255)
                                                         tGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAagc                                                                                              >  1:4069153/1‑150 (MQ=255)
                                                           cGTTACTGCCTTGATGATATTTTCAAAGTCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGACGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGATGATGTCTTCCAGACTCAagcag                                                                                            >  1:3871750/1‑150 (MQ=255)
                                                           cGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAagcag                                                                                            <  1:3389239/150‑1 (MQ=255)
                                                               aCTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGtt                                                                                        <  2:2860488/150‑1 (MQ=255)
                                                               aCTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGtt                                                                                        >  2:2007702/1‑150 (MQ=255)
                                                                cTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGttt                                                                                       <  1:4854194/150‑1 (MQ=255)
                                                                cTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGttt                                                                                       >  2:1860744/1‑150 (MQ=255)
                                                                   ccTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAAc                                                                                    >  2:784714/1‑150 (MQ=255)
                                                                     ttGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAActg                                                                                  <  2:134154/150‑1 (MQ=255)
                                                                     ttGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAActg                                                                                  <  1:2238731/150‑1 (MQ=255)
                                                                      tgatgaTATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAActgc                                                                                 >  1:1288900/1‑150 (MQ=255)
                                                                      tgatgaTATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAActgc                                                                                 <  1:1035498/150‑1 (MQ=255)
                                                                       gatgatATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAActgct                                                                                >  1:4164520/1‑150 (MQ=255)
                                                                       gatgatATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAActgct                                                                                <  1:1809724/150‑1 (MQ=255)
                                                                          gatATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAg                                                                             >  1:2692058/1‑150 (MQ=255)
                                                                          gatATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAg                                                                             >  2:3470161/1‑150 (MQ=255)
                                                                          gatATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAg                                                                             >  2:1017275/1‑150 (MQ=255)
                                                                            tatTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGcc                                                                           >  2:31304/1‑150 (MQ=255)
                                                                            tatTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGcc                                                                           >  2:32990/1‑150 (MQ=255)
                                                                             atTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCg                                                                          >  1:3869258/1‑150 (MQ=255)
                                                                              ttttCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCgg                                                                         >  1:4406704/1‑150 (MQ=255)
                                                                               tttCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGt                                                                        >  2:3581658/1‑150 (MQ=255)
                                                                                  cAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCg                                                                     >  2:3361810/1‑150 (MQ=255)
                                                                                   aaaGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGt                                                                    <  2:1237435/150‑1 (MQ=255)
                                                                                          tcttTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGttt                                                             <  2:2393184/150‑1 (MQ=255)
                                                                                           cttTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTa                                                            <  2:3741894/150‑1 (MQ=255)
                                                                                           cttTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTa                                                            <  2:579343/150‑1 (MQ=255)
                                                                                           cttTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCACGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTa                                                            >  2:4160310/1‑150 (MQ=255)
                                                                                              tGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCa                                                         >  1:4814266/1‑150 (MQ=255)
                                                                                              tGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCa                                                         <  2:444759/150‑1 (MQ=255)
                                                                                               gATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAg                                                        <  1:1070152/150‑1 (MQ=255)
                                                                                                aTTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAgc                                                       <  1:384718/150‑1 (MQ=255)
                                                                                                aTTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAgc                                                       >  1:2803150/1‑150 (MQ=255)
                                                                                                   aTGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGcgcg                                                    >  2:179958/1‑150 (MQ=255)
                                                                                                    tGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGa                                                   <  2:1159828/150‑1 (MQ=255)
                                                                                                     gACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGat                                                  <  2:4357884/150‑1 (MQ=255)
                                                                                                     gACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGat                                                  <  2:2692058/150‑1 (MQ=255)
                                                                                                       cgcgCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGatat                                                <  2:4063108/150‑1 (MQ=255)
                                                                                                         cgcTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGc                                              <  1:3351995/150‑1 (MQ=255)
                                                                                                          gcTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGcc                                             >  1:232777/1‑150 (MQ=255)
                                                                                                           cTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGccc                                            >  1:1744135/1‑150 (MQ=255)
                                                                                                           cTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGccc                                            >  1:3605418/1‑150 (MQ=255)
                                                                                                            tGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCa                                           >  1:1541272/1‑150 (MQ=255)
                                                                                                              aaTGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGACCAAt                                         >  2:4184099/1‑150 (MQ=255)
                                                                                                               aTGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATg                                        >  1:4466532/1‑150 (MQ=255)
                                                                                                                tGCCTATTCAATGAAGATGACCCGCGATGCCGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATgc                                       <  2:71607/150‑1 (MQ=255)
                                                                                                                            gAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAgt                           >  2:4046997/1‑150 (MQ=255)
                                                                                                                             aaGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAgtg                          >  1:4564178/1‑150 (MQ=255)
                                                                                                                                 tGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTAc                      <  1:875249/150‑1 (MQ=255)
                                                                                                                                  gACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTAcg                     <  1:4802184/150‑1 (MQ=255)
                                                                                                                                     ccGCGATGCCGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGa                  >  2:4615816/1‑150 (MQ=255)
                                                                                                                                     ccGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGa                  >  1:1314510/1‑150 (MQ=255)
                                                                                                                                        cgATGCCGAATACTATTTAGTGCATGAGATGTAAGCTAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCTATATGCCCAATGCGCTGGTTGAAGTGTTACGCGaaaa               >  1:3450230/1‑150 (MQ=255)
                                                                                                                                            gCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTg           >  2:4691012/1‑150 (MQ=255)
                                                                                                                                             ccGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGa          >  1:616730/1‑150 (MQ=255)
                                                                                                                                              cGAATACTATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGAc         >  1:2041106/1‑150 (MQ=255)
                                                                                                                                              cGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGAc         >  2:1933748/1‑150 (MQ=255)
                                                                                                                                              cGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGAc         >  2:3162326/1‑150 (MQ=255)
                                                                                                                                               gAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACt        <  2:1795917/150‑1 (MQ=255)
                                                                                                                                                  tACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTAtt     >  2:2531329/1‑150 (MQ=255)
                                                                                                                                                  tACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTAtt     <  2:366228/150‑1 (MQ=255)
                                                                                                                                                  tACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTAtt     >  2:103132/1‑150 (MQ=255)
                                                                                                                                                   aCGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTAttt    >  1:4576991/1‑150 (MQ=255)
                                                                                                                                                   aCGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTAttt    >  1:699845/1‑150 (MQ=255)
                                                                                                                                                    cGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTc   <  1:1275875/150‑1 (MQ=255)
                                                                                                                                                     gATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTcc  >  1:1473282/1‑150 (MQ=255)
                                                                                                                                                     gATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTcc  >  2:3218899/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTTACCGCCAGAAGCGCCGCGTCGACGCAAGCCGATGATTCTTCTGGATAACATTCTGCGTTACTGCCTTGATGATATTTTCAAAGGCTTCTTTGATTATGACGCGCTGAATGCCTATTCAATGAAGATGACCCGCGATGCCGAATACGATTTAGTGCATGAGATGGAAGCCAGCCTGATGGAGTTGATGTCTTCCAGTCTCAAGCAGCGTTTAACTGCTGAGCCGGTGCGTTTTGTTTATCAGCGCGATATGCCCAATGCGCTGGTTGAAGTGTTACGCGAAAAACTGACTATTTCC  >  CP009273/2616953‑2617251

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: