Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,361,433 G→T 11.8% P246Q (CCG→CAG)  pspF ← psp operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,361,4330GT11.8% 765.9 / 65.0 238P246Q (CCG→CAG) pspFpsp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (117/93);  new base T (15/13);  total (132/106)
Fisher's exact test for biased strand distribution p-value = 8.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 242 total aligned reads displayed.

CAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATC  >  CP009273/1361285‑1361582
                                                                                                                                                    |                                                                                                                                                     
cAAACTGAGGTGCAGCAACTCTTGTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGTGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCgg                                                                                                                                                      <  2:2653848/150‑1 (MQ=255)
 aaaCTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGa                                                                                                                                                     >  1:1080547/1‑150 (MQ=255)
 aaaCTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGa                                                                                                                                                     >  2:2861965/1‑150 (MQ=255)
  aaCTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGAt                                                                                                                                                    <  2:3969312/150‑1 (MQ=255)
  aaCTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGAt                                                                                                                                                    <  1:2558795/150‑1 (MQ=255)
    cTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATaa                                                                                                                                                  <  2:2377381/150‑1 (MQ=255)
    cTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATaa                                                                                                                                                  >  1:2950929/1‑150 (MQ=255)
     tGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAAt                                                                                                                                                 >  1:888109/1‑150 (MQ=255)
     tGAGTCGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAAt                                                                                                                                                 >  1:3513963/1‑150 (MQ=255)
       aGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCg                                                                                                                                               >  2:1759374/1‑150 (MQ=255)
         ttGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCt                                                                                                                                             <  2:1239818/150‑1 (MQ=255)
          tGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTg                                                                                                                                            >  1:763467/1‑150 (MQ=255)
            cagcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGt                                                                                                                                          <  1:1054277/150‑1 (MQ=255)
             agcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTg                                                                                                                                         >  2:3011705/1‑150 (MQ=255)
                aaCTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGAGGATAATCGCTGGTGCCg                                                                                                                                      >  1:1844022/1‑150 (MQ=255)
                aaCTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGGGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAATGGGATCAATAATGATGTCATCAGGCGGATAATCGCTGGTGCCg                                                                                                                                      >  2:2719860/1‑150 (MQ=255)
                   tctTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGGTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTgg                                                                                                                                   >  1:1152538/1‑150 (MQ=255)
                   tctTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGGTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTgg                                                                                                                                   >  2:3700695/1‑150 (MQ=255)
                     tttttCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCg                                                                                                                                 >  1:1944223/1‑150 (MQ=255)
                      ttttCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGa                                                                                                                                >  2:4631768/1‑150 (MQ=255)
                       tttCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGAt                                                                                                                               >  2:5097772/1‑150 (MQ=255)
                        ttCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGTCGGACGCCGTGTAAAGGGATCAATAATTATGTCATCAAGCGGATAATCGATGGTGACGTGGCGATa                                                                                                                              >  2:5062906/1‑150 (MQ=255)
                        ttCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATa                                                                                                                              <  1:4694759/150‑1 (MQ=255)
                            tgctgcATCTGAAACTCACGTAAATCCCGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGAGAGCGACTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGAGGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACt                                                                                                                          <  2:4513794/150‑1 (MQ=255)
                             gctgcATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGCTAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGGGGCGAGACACTg                                                                                                                         <  2:569335/150‑1 (MQ=255)
                               tgcATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGaa                                                                                                                       >  1:1698016/1‑150 (MQ=255)
                               tgcATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAAGCGCTGGTGCCGAGGCGATACACTGaa                                                                                                                       >  2:3306026/1‑150 (MQ=255)
                                 cATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACg                                                                                                                     >  1:2871175/1‑150 (MQ=255)
                                  aTCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGt                                                                                                                    >  2:21655/1‑150 (MQ=255)
                                   tCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAAGGtt                                                                                                                   <  2:1967756/150‑1 (MQ=255)
                                      gAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTcca                                                                                                                >  1:2736938/1‑150 (MQ=255)
                                      gAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTcca                                                                                                                <  1:5129524/150‑1 (MQ=255)
                                          cTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCcaccac                                                                                                            >  1:4265680/1‑150 (MQ=255)
                                           tCACGTAAATCCAGCGGCAGAGTTGGAAGCGATGTTGTTTCTGAAACTGCGATAGCGTCTTCAGGCCGACGCCGTAAAAAGGGAGCAATAATGAAGTCATCAAGCGCATAATCGCAGTTGCCGAGCAGATACACTGAACGTACCACCACg                                                                                                           >  2:2712257/1‑150 (MQ=255)
                                            cACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGt                                                                                                          >  2:3249709/1‑150 (MQ=255)
                                             aCGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGtt                                                                                                         <  2:3553537/150‑1 (MQ=255)
                                               gTAAATCCAGCGGCAGTTTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGtttt                                                                                                       <  2:4576852/150‑1 (MQ=255)
                                               gTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGtttt                                                                                                       <  2:2945728/150‑1 (MQ=255)
                                                tAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCTGTGGCGATACACTGAACGTTCCACCACGttttt                                                                                                      <  1:3927697/150‑1 (MQ=255)
                                                 aaaTCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTc                                                                                                     <  2:4503520/150‑1 (MQ=255)
                                                      cAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTc                                                                                                <  1:2453287/150‑1 (MQ=255)
                                                       aGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCa                                                                                               >  1:1771988/1‑150 (MQ=255)
                                                         cGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACg                                                                                             <  2:73230/150‑1 (MQ=255)
                                                         cGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACg                                                                                             <  2:316600/150‑1 (MQ=255)
                                                         cGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACg                                                                                             <  2:3714762/150‑1 (MQ=255)
                                                         cGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACg                                                                                             >  2:1404058/1‑150 (MQ=255)
                                                          ggCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGa                                                                                            <  1:1317782/150‑1 (MQ=255)
                                                           gCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGaa                                                                                           <  2:2473341/150‑1 (MQ=255)
                                                           gCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGaa                                                                                           <  1:5171646/150‑1 (MQ=255)
                                                                gtTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATAttt                                                                                      <  2:5020269/150‑1 (MQ=255)
                                                                  tGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTcc                                                                                    <  1:14262/150‑1 (MQ=255)
                                                                    gAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCg                                                                                  >  2:108194/1‑150 (MQ=255)
                                                                         gAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCaa                                                                             >  2:800696/1‑150 (MQ=255)
                                                                         gAGGGGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCaa                                                                             >  2:800841/1‑150 (MQ=255)
                                                                          aGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAAc                                                                            >  1:3335563/1‑150 (MQ=255)
                                                                          aGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAAc                                                                            <  2:3087485/150‑1 (MQ=255)
                                                                             tggtTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGAt                                                                         <  1:3170913/150‑1 (MQ=255)
                                                                             tggtTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGAt                                                                         >  2:4291673/1‑150 (MQ=255)
                                                                             tggtTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGAt                                                                         <  1:3771582/150‑1 (MQ=255)
                                                                             tggtTTCTGAAACGGCGATAGCGTATTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCAACACGTTTTTCATTTCACGAATATTTCCCGGCCAACgtt                                                                         <  1:671475/150‑3 (MQ=255)
                                                                                     gAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGc                                                                 <  1:1900393/150‑1 (MQ=255)
                                                                                       aaCGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCaa                                                               >  1:1244303/1‑150 (MQ=255)
                                                                                        aCGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAAt                                                              <  1:2337956/150‑1 (MQ=255)
                                                                                         cGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATg                                                             <  1:4142719/150‑1 (MQ=255)
                                                                                          ggCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGt                                                            >  2:3537431/1‑150 (MQ=255)
                                                                                            cGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGttt                                                          >  2:1483249/1‑150 (MQ=255)
                                                                                            cGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTCCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGttt                                                          <  1:4585723/150‑1 (MQ=255)
                                                                                             gATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTtc                                                         <  2:4591987/150‑1 (MQ=255)
                                                                                             gATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTtc                                                         >  1:4161192/1‑150 (MQ=255)
                                                                                             gATAGCGTCTTCAGGCGGACGCCGATTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTtc                                                         <  1:4356632/150‑1 (MQ=255)
                                                                                               tAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTtctc                                                       <  2:700137/150‑1 (MQ=255)
                                                                                               tAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTtctc                                                       <  2:2889239/150‑1 (MQ=255)
                                                                                               tAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTtctc                                                       >  1:2638704/1‑150 (MQ=255)
                                                                                                   gTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATCATTCAGCAATGTTTCTCTGgc                                                   >  2:2556809/1‑150 (MQ=255)
                                                                                                    tCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGgcg                                                  >  2:2629163/1‑150 (MQ=255)
                                                                                                      ttCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGgcgcg                                                <  2:62313/150‑1 (MQ=255)
                                                                                                        cAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCt                                              >  2:2488959/1‑150 (MQ=255)
                                                                                                         aGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTc                                             >  2:145573/1‑150 (MQ=255)
                                                                                                           gCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCg                                           <  2:5166982/150‑1 (MQ=255)
                                                                                                             ggACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTa                                         <  2:1266995/150‑1 (MQ=255)
                                                                                                             ggACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTa                                         >  2:168401/1‑150 (MQ=255)
                                                                                                             ggACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTATTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTa                                         >  2:1381526/1‑150 (MQ=255)
                                                                                                                  ccGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAAccc                                    >  1:3685530/1‑150 (MQ=255)
                                                                                                                  ccGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAAccc                                    <  2:4792168/150‑1 (MQ=255)
                                                                                                                   cGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAAcccc                                   <  1:4459929/150‑1 (MQ=255)
                                                                                                                    gTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCg                                  >  2:4096465/1‑150 (MQ=255)
                                                                                                                    gTTAAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCg                                  >  1:1162354/1‑150 (MQ=255)
                                                                                                                     tttAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCgg                                 >  2:3568599/1‑150 (MQ=255)
                                                                                                                      ttAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCggg                                >  1:4827744/1‑150 (MQ=255)
                                                                                                                      ttAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCggg                                >  1:1095926/1‑150 (MQ=255)
                                                                                                                           gggATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACag                           <  1:786224/150‑1 (MQ=255)
                                                                                                                             gATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACagag                         >  1:353017/1‑150 (MQ=255)
                                                                                                                                cAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCa                      >  1:130086/1‑150 (MQ=255)
                                                                                                                                cAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCa                      >  1:3489039/1‑150 (MQ=255)
                                                                                                                                  ataatGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGc                    >  1:5041708/1‑150 (MQ=255)
                                                                                                                                          gTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTc            >  2:3559773/1‑150 (MQ=255)
                                                                                                                                               caAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGac       <  1:3668221/150‑1 (MQ=255)
                                                                                                                                                aaGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGaca      <  2:2016004/150‑1 (MQ=255)
                                                                                                                                                aaGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTGTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGaca      >  1:144741/1‑150 (MQ=255)
                                                                                                                                                   cGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACAt   >  1:4541193/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATC  >  CP009273/1361285‑1361582

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