Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,361,433 G→T 12.7% 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,4330GT12.7% 920.0 / 84.5 266P246Q (CCG→CAG) pspFpsp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (115/117);  new base T (16/18);  total (131/135)
Fisher's exact test for biased strand distribution p-value = 8.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 272 total aligned reads displayed.

CAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACAT  >  CP009273/1361285‑1361581
                                                                                                                                                    |                                                                                                                                                    
cAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCg                                                                                                                                                      >  1:4248306/1‑149 (MQ=255)
cAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCg                                                                                                                                                      >  1:3984323/1‑149 (MQ=255)
cAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCg                                                                                                                                                      >  1:3983037/1‑149 (MQ=255)
  aaCTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGa                                                                                                                                                    <  2:2977030/149‑1 (MQ=255)
  aaCTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGa                                                                                                                                                    >  1:989506/1‑149 (MQ=255)
  aaCTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGa                                                                                                                                                    >  1:1702073/1‑149 (MQ=255)
  aaCTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGa                                                                                                                                                    >  1:1319150/1‑149 (MQ=255)
      gAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAAt                                                                                                                                                <  2:4276223/149‑1 (MQ=255)
        gTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTCAAGGGATCAATAATGATGTCATCAAGCTGATACTCg                                                                                                                                              >  2:3525845/1‑149 (MQ=255)
        gTTGCAGCAACTCTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGc                                                                                                                                             >  1:562667/1‑149 (MQ=255)
           gcagcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTg                                                                                                                                           >  2:2868505/1‑149 (MQ=255)
           gcagcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTg                                                                                                                                           <  1:3236282/149‑1 (MQ=255)
           gcagcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTg                                                                                                                                           <  1:2868505/149‑1 (MQ=255)
             agcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGt                                                                                                                                         >  2:220349/1‑149 (MQ=255)
             agcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGt                                                                                                                                         >  2:4627424/1‑149 (MQ=255)
              gcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTg                                                                                                                                        >  1:4832570/1‑149 (MQ=255)
              gcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTg                                                                                                                                        >  1:3849256/1‑149 (MQ=255)
              gcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTg                                                                                                                                        >  1:2236525/1‑149 (MQ=255)
              gcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTg                                                                                                                                        >  1:1792581/1‑149 (MQ=255)
               caACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGc                                                                                                                                       <  2:1010480/149‑1 (MQ=255)
                  ctctTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGt                                                                                                                                    <  2:1637597/149‑1 (MQ=255)
                    ctTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTgg                                                                                                                                  <  2:1877864/149‑1 (MQ=255)
                         tCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATa                                                                                                                             >  1:4729478/1‑149 (MQ=255)
                         tCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATa                                                                                                                             <  2:4387405/149‑1 (MQ=255)
                          ccTGCTGCATCTGAATCTCACGTAAATCCAGCGGCAGTGTTGGTAGCTAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATac                                                                                                                            <  2:3774839/149‑1 (MQ=255)
                          ccTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATac                                                                                                                            <  2:4126542/149‑1 (MQ=255)
                             gctgcATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACt                                                                                                                         <  2:2201425/149‑1 (MQ=255)
                                gcATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGaa                                                                                                                      <  1:859083/149‑1 (MQ=255)
                                 cATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAAc                                                                                                                     <  2:1363928/149‑1 (MQ=255)
                                   tCTGAAACTCACTTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGt                                                                                                                   <  1:2657464/149‑1 (MQ=255)
                                     tGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTc                                                                                                                 <  1:735898/149‑1 (MQ=255)
                                     tGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTc                                                                                                                 <  1:588041/149‑1 (MQ=255)
                                       aaaCTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTcca                                                                                                               >  1:1459991/1‑149 (MQ=255)
                                       aaaCTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTcca                                                                                                               >  2:4030850/1‑149 (MQ=255)
                                       aaaCTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTcca                                                                                                               >  2:1434333/1‑149 (MQ=255)
                                         aCTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTccacc                                                                                                             <  1:3822892/149‑1 (MQ=255)
                                          cTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTccacca                                                                                                            >  2:1116855/1‑149 (MQ=255)
                                                  aaTCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTc                                                                                                    <  1:2445165/149‑1 (MQ=255)
                                                     ccAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAAt                                                                                                 <  2:2280823/149‑1 (MQ=255)
                                                      cAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAAtt                                                                                                <  1:1862582/149‑1 (MQ=255)
                                                           gCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGa                                                                                           >  1:4683125/1‑149 (MQ=255)
                                                             aGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAat                                                                                         >  1:1095431/1‑149 (MQ=255)
                                                               tgtTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAatat                                                                                       <  1:3595517/149‑1 (MQ=255)
                                                               tgtTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAatat                                                                                       >  1:3364250/1‑149 (MQ=255)
                                                                gtTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATAtt                                                                                      >  1:4323941/1‑149 (MQ=255)
                                                                  tGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTc                                                                                    <  2:1958420/149‑1 (MQ=255)
                                                                   ggAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTcc                                                                                   >  2:3161012/1‑149 (MQ=255)
                                                                    gAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTccc                                                                                  <  1:2482289/149‑1 (MQ=255)
                                                                        cGAGTTGTTTTCTGAAACGTCGATAGCGTCTTCAGGGGGACGCAGTTTATAGGGATCAATAATGATGTCATCAAGCGTATATTCGCTGGTGCCGTGGAGATACACTTAACGTTCCTCCACGTTTTTCAATTCACGAATATTTCCCGGcc                                                                              <  2:4775172/149‑1 (MQ=255)
                                                                         gAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCa                                                                             <  1:3086403/149‑1 (MQ=255)
                                                                           ggtggtTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAAc                                                                           <  2:1267453/149‑1 (MQ=255)
                                                                             tggtTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGa                                                                         <  2:4116533/149‑1 (MQ=255)
                                                                             tggtTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGa                                                                         <  2:660814/149‑1 (MQ=255)
                                                                             tggtTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGa                                                                         <  1:3379602/149‑1 (MQ=255)
                                                                                 ttCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAAt                                                                     <  1:2630741/149‑1 (MQ=255)
                                                                                  tCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAAtt                                                                    <  1:1140975/149‑1 (MQ=255)
                                                                                  tCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAAtt                                                                    <  2:5139253/149‑1 (MQ=255)
                                                                                      aaaCGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGc                                                                <  2:1889726/149‑1 (MQ=255)
                                                                                      aaaCGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGc                                                                <  2:1775124/149‑1 (MQ=255)
                                                                                        aCGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCaa                                                              >  2:5177075/1‑149 (MQ=255)
                                                                                        aCGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCaa                                                              <  1:2067205/149‑1 (MQ=255)
                                                                                         cGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAAt                                                             <  2:689022/149‑1 (MQ=255)
                                                                                          ggCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATg                                                            <  1:3478227/149‑1 (MQ=255)
                                                                                           gCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGt                                                           <  1:1811343/149‑1 (MQ=255)
                                                                                            cGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGtt                                                          <  2:2472435/149‑1 (MQ=255)
                                                                                            cGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATGTCCCGGCCAACGATAATTCAGCAATGtt                                                          <  1:5125202/149‑1 (MQ=255)
                                                                                             gATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGttt                                                         >  2:448289/1‑149 (MQ=255)
                                                                                                aGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTtctc                                                      <  2:798688/149‑1 (MQ=255)
                                                                                                 gCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTctct                                                     <  1:375629/149‑1 (MQ=255)
                                                                                                    tCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGgc                                                  >  2:2104327/1‑149 (MQ=255)
                                                                                                    tCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGgc                                                  >  1:1073317/1‑149 (MQ=255)
                                                                                                    tCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAAGCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGgc                                                  >  1:691213/1‑149 (MQ=255)
                                                                                                      ttCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGgcgc                                                >  1:1568040/1‑149 (MQ=255)
                                                                                                      ttCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGgcgc                                                <  2:3645631/149‑1 (MQ=255)
                                                                                                        cAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGgcgcgc                                              <  2:866737/149‑1 (MQ=255)
                                                                                                          ggCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTc                                            >  1:1397490/1‑149 (MQ=255)
                                                                                                          ggCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTc                                            >  2:3313774/1‑149 (MQ=255)
                                                                                                          ggCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTc                                            >  1:5050295/1‑149 (MQ=255)
                                                                                                           gCGGACGCCGTTTAAAGGGATCAGTAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCGGCAATGTTTCTCTGGCGCGCTcc                                           >  1:3575082/1‑149 (MQ=255)
                                                                                                           gCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTcc                                           >  1:3317615/1‑149 (MQ=255)
                                                                                                                 gCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAAc                                     >  1:5182679/1‑149 (MQ=255)
                                                                                                                  ccGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAAcc                                    >  2:4965678/1‑149 (MQ=255)
                                                                                                                    gTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAAcccc                                  >  1:3192778/1‑149 (MQ=255)
                                                                                                                    gTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAAcccc                                  >  2:2667848/1‑149 (MQ=255)
                                                                                                                    gTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAAcccc                                  >  2:2580076/1‑149 (MQ=255)
                                                                                                                        aaaGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGa                              >  1:1024851/1‑149 (MQ=255)
                                                                                                                         aaGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGaa                             <  1:3243001/149‑1 (MQ=255)
                                                                                                                         aaGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGaa                             <  1:3018916/149‑1 (MQ=255)
                                                                                                                         aaGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGaa                             >  2:4916606/1‑149 (MQ=255)
                                                                                                                           gggATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACa                           >  1:1014929/1‑149 (MQ=255)
                                                                                                                              aTCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACagag                        <  2:2931869/149‑1 (MQ=255)
                                                                                                                              aTCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACagag                        <  2:4717170/149‑1 (MQ=255)
                                                                                                                               tCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAgg                       <  1:2128150/149‑1 (MQ=255)
                                                                                                                                cAATCATGATGTCATCTAGCGGATAATCGCTGGTGCCTTGGCGATACACTGAACTTTCCACCACGTTTTTCTATTCTCGAATATGGCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGATGc                      <  2:199757/149‑1 (MQ=255)
                                                                                                                                    aatGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCt                  <  1:196211/149‑1 (MQ=255)
                                                                                                                                           tcatcaAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTc           <  2:3965327/149‑1 (MQ=255)
                                                                                                                                           tcatcaAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTc           <  1:2802916/149‑1 (MQ=255)
                                                                                                                                            catcaAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTcc          <  1:4462883/149‑1 (MQ=255)
                                                                                                                                                aaGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGac      >  2:3877751/1‑149 (MQ=255)
                                                                                                                                                 aGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGaca     <  2:2053442/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACAT  >  CP009273/1361285‑1361581

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: