Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,361,433 G→T 5.6% 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,4330GT 5.6% 1163.6 / 38.3 304P246Q (CCG→CAG) pspFpsp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (141/146);  new base T (6/11);  total (147/157)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 311 total aligned reads displayed.

CAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATC  >  CP009273/1361285‑1361582
                                                                                                                                                    |                                                                                                                                                     
cAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCgg                                                                                                                                                      >  2:2916256/1‑150 (MQ=255)
cAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCgg                                                                                                                                                      <  1:24188/150‑1 (MQ=255)
cAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCgg                                                                                                                                                      >  1:2146650/1‑150 (MQ=255)
cAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGGTTAAAGGGATCAATAATGATGTCATCAAGCgg                                                                                                                                                      >  2:148111/1‑150 (MQ=255)
 aaaCTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGa                                                                                                                                                     >  2:534063/1‑150 (MQ=255)
     tGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAAt                                                                                                                                                 >  1:3940311/1‑150 (MQ=255)
     tGAGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAAt                                                                                                                                                 <  2:3839622/150‑1 (MQ=255)
     agAGTTGCAGCAACTCTTTTGCCTGCTGCAGATGAAACTCACGTAAATCCCGCGGCAGAGTTGGACGCGAGGTGGTATCTGAAACGGCGAGAGCGTCTTCAGGCGGACGCCGTTTACAGGGAACACTAATGCGGTCCTCAAGCGGAtgag                                                                                                                                                 >  2:2919853/2‑147 (MQ=255)
       aGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCg                                                                                                                                               <  1:3551827/150‑1 (MQ=255)
       aGTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCg                                                                                                                                               <  1:4768907/150‑1 (MQ=255)
        gTTGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGc                                                                                                                                              >  1:3784449/1‑150 (MQ=255)
        gTTGCAGCAACTCTTTTTCCAGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGACGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGGTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGc                                                                                                                                              >  2:11303/1‑150 (MQ=255)
         ttGCAGCAACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCt                                                                                                                                             <  1:2665992/150‑1 (MQ=255)
            cagcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTTTTGGAAGCGATGTGGTTTCTGAAAAGGCGATAGCGTCTTCAGGCGGAAGCCGTTTAAAGGGATAAATAATGATGTAATCAATCGTATAATCGCTgga                                                                                                                                          >  2:817018/1‑149 (MQ=255)
             agcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTg                                                                                                                                         >  1:3032062/1‑150 (MQ=255)
              gcaACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGc                                                                                                                                        >  2:1975089/1‑150 (MQ=255)
               caACTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGcc                                                                                                                                       >  1:724434/1‑150 (MQ=255)
                aaCTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCg                                                                                                                                      <  1:1904185/150‑1 (MQ=255)
                aaCTCTTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCg                                                                                                                                      >  1:1382828/1‑150 (MQ=255)
                  ctctTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTg                                                                                                                                    >  1:1541763/1‑150 (MQ=255)
                   tctTTTTCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTgg                                                                                                                                   <  1:650886/150‑1 (MQ=255)
                       tttCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGAt                                                                                                                               >  2:3263190/1‑150 (MQ=255)
                        ttCCTGCTGCATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATa                                                                                                                              <  1:4118004/150‑1 (MQ=255)
                            tgctgcATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACt                                                                                                                          <  2:721200/150‑1 (MQ=255)
                            tgctgcATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACt                                                                                                                          <  1:871503/150‑1 (MQ=255)
                             gctgcATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTg                                                                                                                         <  2:4633434/150‑1 (MQ=255)
                              ctgcATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGa                                                                                                                        >  2:2267424/1‑150 (MQ=255)
                                tcATCTTAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGTGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAAc                                                                                                                      <  1:5357366/149‑1 (MQ=255)
                                 cATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACg                                                                                                                     >  1:3461231/1‑150 (MQ=255)
                                 cATCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACg                                                                                                                     <  1:514077/150‑1 (MQ=255)
                                  aTCTGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGt                                                                                                                    >  1:1020575/1‑150 (MQ=255)
                                     tGAAACTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGTGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTcc                                                                                                                 >  2:4151498/1‑150 (MQ=255)
                                         aCTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTccacca                                                                                                             >  2:4237012/1‑150 (MQ=255)
                                         aCTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTccacca                                                                                                             <  1:100075/150‑1 (MQ=255)
                                          cTCACGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCcaccac                                                                                                            <  2:1620516/150‑1 (MQ=255)
                                             aCGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGtt                                                                                                         >  2:2342854/1‑150 (MQ=255)
                                              cGTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGttt                                                                                                        <  2:3710694/150‑1 (MQ=255)
                                               gTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGtttt                                                                                                       <  2:1327904/150‑1 (MQ=255)
                                               gTAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGtttt                                                                                                       >  1:1600993/1‑150 (MQ=255)
                                                tAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGttttt                                                                                                      <  2:1569309/150‑1 (MQ=255)
                                                tAAATCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGttttt                                                                                                      <  2:4297751/150‑1 (MQ=255)
                                                 aaaTCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTc                                                                                                     >  2:3450520/1‑150 (MQ=255)
                                                 aaaTCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTc                                                                                                     <  1:4850139/150‑1 (MQ=255)
                                                   aTCCAGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCaa                                                                                                   >  2:4207595/1‑150 (MQ=255)
                                                       aGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCa                                                                                               >  1:3302444/1‑150 (MQ=255)
                                                       aGCGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCa                                                                                               <  1:3185617/150‑1 (MQ=255)
                                                         cGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGCTCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACg                                                                                             >  2:88729/1‑150 (MQ=255)
                                                         cGGCAGTGTTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACg                                                                                             <  2:5067762/150‑1 (MQ=255)
                                                              gtgtTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAatat                                                                                        <  2:4590184/150‑1 (MQ=255)
                                                              gtgtTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAatat                                                                                        >  1:3191493/1‑150 (MQ=255)
                                                                gtTGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATAttt                                                                                      <  2:4092528/150‑1 (MQ=255)
                                                                  tGGAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTcc                                                                                    <  1:108854/150‑1 (MQ=255)
                                                                    gAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCg                                                                                  >  1:1187743/1‑150 (MQ=255)
                                                                    gAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCg                                                                                  <  1:1742111/150‑1 (MQ=255)
                                                                    gAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCg                                                                                  >  1:5009602/1‑150 (MQ=255)
                                                                    gAAGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCg                                                                                  <  1:3366816/150‑1 (MQ=255)
                                                                      aGCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGACGTATAAAGGGATCAATAATGATGTCATCAAGCGGATAAACGCTGGTGCAGTCGCCATACACTGAACGATCCACCACGTTTTACAATTCACGAATATTTccctga                                                                                >  2:564719/1‑147 (MQ=255)
                                                                       gCGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGcc                                                                               <  1:2716706/150‑1 (MQ=255)
                                                                        cGAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACACCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCa                                                                              >  1:3103507/1‑150 (MQ=255)
                                                                         gAGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCaa                                                                             >  1:4172114/1‑150 (MQ=255)
                                                                          aGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAAc                                                                            >  2:2314699/1‑150 (MQ=255)
                                                                          aGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAAc                                                                            <  2:244197/150‑1 (MQ=255)
                                                                          aGGTGGTTTCTGAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAAc                                                                            >  1:2951763/1‑150 (MQ=255)
                                                                                     gAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGc                                                                 <  1:5195638/150‑1 (MQ=255)
                                                                                     gAAACGGCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGc                                                                 <  1:4923939/150‑1 (MQ=255)
                                                                                          ggCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGt                                                            <  2:2222080/150‑1 (MQ=255)
                                                                                          ggCGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGt                                                            >  2:1084058/1‑150 (MQ=255)
                                                                                            cGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGttt                                                          <  2:4707021/150‑1 (MQ=255)
                                                                                            cGATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGttt                                                          <  1:2805188/150‑1 (MQ=255)
                                                                                             gATAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTtc                                                         >  1:454085/1‑150 (MQ=255)
                                                                                              aTAGCGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTtct                                                        >  2:4122736/1‑150 (MQ=255)
                                                                                                aGCGTCTTCAGGCGGACGCAGTTTAAAGGGATCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTctct                                                      <  1:3518963/150‑1 (MQ=255)
                                                                                                  cGTCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTgg                                                    <  1:5462459/150‑1 (MQ=255)
                                                                                                    tCTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGgcg                                                  <  1:5368901/150‑1 (MQ=255)
                                                                                                     cTTCAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGgcgc                                                 >  1:2326245/1‑150 (MQ=255)
                                                                                                        cAGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCt                                              <  2:2511549/150‑1 (MQ=255)
                                                                                                         aGGCGGACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTc                                             >  2:1923431/1‑150 (MQ=255)
                                                                                                             ggACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTa                                         >  1:1552552/1‑150 (MQ=255)
                                                                                                              gACGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTaa                                        <  1:968069/150‑1 (MQ=255)
                                                                                                               aCGCCGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTaaa                                       >  2:4766457/1‑150 (MQ=255)
                                                                                                                   cGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAAcccc                                   >  2:872972/1‑150 (MQ=255)
                                                                                                                   cGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAAcccc                                   >  2:75809/1‑150 (MQ=255)
                                                                                                                   cGTTTAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAAcccc                                   <  1:1208466/150‑1 (MQ=255)
                                                                                                                   cGTTTAAAGGGAGCAATAATGATGTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAAcccc                                   <  2:2939345/150‑1 (MQ=255)
                                                                                                                      ttAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCggg                                >  2:3398361/1‑150 (MQ=255)
                                                                                                                       tAAAGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGa                               <  2:2078635/150‑1 (MQ=255)
                                                                                                                          aGGGATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACa                            <  2:2178139/150‑1 (MQ=255)
                                                                                                                           gggATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGACCag                           >  2:2629966/1‑150 (MQ=255)
                                                                                                                            ggATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACaga                          >  1:506782/1‑150 (MQ=255)
                                                                                                                            ggATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACaga                          >  1:2545004/1‑150 (MQ=255)
                                                                                                                             gATCAATAATGATGTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACagag                         <  2:503707/150‑1 (MQ=255)
                                                                                                                                     atgatgTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTg                 <  1:2619808/150‑1 (MQ=255)
                                                                                                                                       gatgTCATCAAGCTGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCTGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGAt               >  2:5206602/1‑150 (MQ=255)
                                                                                                                                         tgTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGAttt             >  2:2389902/1‑150 (MQ=255)
                                                                                                                                         tgTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGAttt             <  1:4039153/150‑1 (MQ=255)
                                                                                                                                         tgTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGAttt             <  2:3590577/150‑1 (MQ=255)
                                                                                                                                          gTCATCAAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTc            >  2:4746879/1‑150 (MQ=255)
                                                                                                                                           tcatcaAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTcc           >  1:4474954/1‑150 (MQ=255)
                                                                                                                                               caAGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGac       >  1:920239/1‑150 (MQ=255)
                                                                                                                                                aaGCGGATAATCGCTGGTGCCGTGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTCGATTTCCCGaca      <  1:1491309/150‑1 (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: