Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 321,743 G→A 6.2% intergenic (+28/+430) aguA → / ← PP_0267 agmatine deiminase/ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947321,7430GA 6.2% 603.5 / 18.0 192intergenic (+28/+430)aguA/PP_0267agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (88/92);  new base A (6/6);  total (94/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 197 total aligned reads displayed.

GATCCTGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACGGGCCAG  >  NC_002947/321595‑321892
                                                                                                                                                    |                                                                                                                                                     
gctCCCGGCAAGGCTCGTCCCCGACCATGAGGTGCTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCAtgt                                                                                                                                                      <  2:5236620/148‑1 (MQ=255)
 aTCCTGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCAtgtg                                                                                                                                                     >  2:5277514/1‑150 (MQ=255)
    cTGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACg                                                                                                                                                  <  2:1076516/150‑1 (MQ=255)
     tGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACgg                                                                                                                                                 >  2:131624/1‑150 (MQ=255)
     tGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACgg                                                                                                                                                 >  2:4135984/1‑150 (MQ=255)
     tGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACgg                                                                                                                                                 >  2:1512655/1‑150 (MQ=255)
     tGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACgg                                                                                                                                                 >  1:858309/1‑150 (MQ=255)
      ggCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACggg                                                                                                                                                <  2:2650811/150‑1 (MQ=255)
       gCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCACATGACGGGc                                                                                                                                               <  2:2798502/150‑1 (MQ=255)
       gCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCACATGACGGGc                                                                                                                                               <  2:2001046/150‑1 (MQ=255)
         aaGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCct                                                                                                                                             >  1:2652263/1‑150 (MQ=255)
         aaGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCct                                                                                                                                             >  1:4880523/1‑150 (MQ=255)
         aaGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCct                                                                                                                                             <  1:4721633/150‑1 (MQ=255)
         aaGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCct                                                                                                                                             >  2:5212184/1‑150 (MQ=255)
         aaGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCct                                                                                                                                             >  1:5473651/1‑150 (MQ=255)
         aaGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCACATGACGGGCct                                                                                                                                             <  2:3641989/150‑1 (MQ=255)
          aGGATCTTCCCAGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCctc                                                                                                                                            >  1:3156274/1‑150 (MQ=255)
              tCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCtttt                                                                                                                                        >  2:1883460/1‑150 (MQ=255)
                ttCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAGGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCtttttt                                                                                                                                      <  1:5530002/150‑1 (MQ=255)
                ttCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCtttttt                                                                                                                                      <  1:3810540/150‑1 (MQ=255)
                ttCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCtttttt                                                                                                                                      <  2:3279395/150‑1 (MQ=255)
                gtCCCCGCCCATGAGGTGGTGATGATCCCGGCCCGCGAGCTGCACTTGGGGGGCGGCAACATTCACTGCCAGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCtttttt                                                                                                                                      <  2:1771509/149‑1 (MQ=255)
                     cGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGcgcg                                                                                                                                 >  2:1742418/1‑150 (MQ=255)
                                   tgatgaTCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGaa                                                                                                                   <  2:1120676/150‑1 (MQ=255)
                                    gatgatCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAg                                                                                                                  >  2:4611881/1‑150 (MQ=255)
                                    gatgatCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAg                                                                                                                  >  1:257884/1‑150 (MQ=255)
                                     atgatCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGc                                                                                                                 >  1:1560762/1‑150 (MQ=255)
                                      tgatCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGcc                                                                                                                <  2:4558571/150‑1 (MQ=255)
                                        atCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGa                                                                                                              <  2:2724812/150‑1 (MQ=255)
                                         tCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAg                                                                                                             >  1:4788809/1‑150 (MQ=255)
                                          cccGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGc                                                                                                            <  2:1632443/150‑1 (MQ=255)
                                          accGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTATTGCGCGGGCGACGCCCGGAAGCCGAGc                                                                                                            >  2:3840243/2‑150 (MQ=255)
                                           ccGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCt                                                                                                           <  2:2683151/150‑1 (MQ=255)
                                             gggCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTgg                                                                                                         <  2:288227/150‑1 (MQ=255)
                                                  gcgAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATAc                                                                                                    <  2:3438898/150‑1 (MQ=255)
                                                    gAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGt                                                                                                  >  2:650819/1‑150 (MQ=255)
                                                           tctTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTAt                                                                                           >  1:5092687/1‑150 (MQ=255)
                                                            ctTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTAtt                                                                                          <  1:1174958/150‑1 (MQ=255)
                                                             ttGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTg                                                                                         >  1:4316079/1‑150 (MQ=255)
                                                              tGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGt                                                                                        <  1:800361/150‑1 (MQ=255)
                                                              ggggTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGt                                                                                        <  2:3709593/149‑1 (MQ=255)
                                                                 gTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGtttt                                                                                     >  1:5215162/1‑150 (MQ=255)
                                                                 gTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGtttt                                                                                     >  2:3473428/1‑150 (MQ=255)
                                                                   ggcggcAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTc                                                                                   <  1:921814/150‑1 (MQ=255)
                                                                   ggcggcAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTc                                                                                   <  2:1230140/150‑1 (MQ=255)
                                                                   ggcggcAACAGTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTc                                                                                   <  1:4331190/150‑1 (MQ=255)
                                                                     cggcAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCaa                                                                                 <  2:593726/150‑1 (MQ=255)
                                                                      ggcAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAg                                                                                >  1:362885/1‑150 (MQ=255)
                                                                            aTTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATa                                                                          <  1:1110862/150‑1 (MQ=255)
                                                                            aTTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATa                                                                          <  1:3691191/150‑1 (MQ=255)
                                                                              tCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGc                                                                        >  2:5273933/1‑150 (MQ=255)
                                                                                aCTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAg                                                                         <  1:1076697/147‑1 (MQ=255)
                                                                                aCTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCaa                                                                      <  2:4942480/150‑1 (MQ=255)
                                                                                aCTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCaa                                                                      <  2:3591152/150‑1 (MQ=255)
                                                                                aCTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCaa                                                                      >  1:3061580/1‑150 (MQ=255)
                                                                                  tGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGa                                                                    >  2:1566831/1‑150 (MQ=255)
                                                                                   gCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGAt                                                                   <  1:3235158/150‑1 (MQ=255)
                                                                                    ccTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATg                                                                  <  1:2996726/150‑1 (MQ=255)
                                                                                       gACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCa                                                               >  1:4552450/1‑150 (MQ=255)
                                                                                       gACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCa                                                               <  1:2272787/150‑1 (MQ=255)
                                                                                        aCCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCAc                                                              >  1:3881003/1‑150 (MQ=255)
                                                                                          ccAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATCCGGTTTGTATTGTTTTTCATGCAGATAGCAAGAGGGCACgg                                                            >  1:5340513/1‑150 (MQ=255)
                                                                                          ccAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACgg                                                            >  1:836893/1‑150 (MQ=255)
                                                                                           cAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGt                                                           <  1:25115/150‑1 (MQ=255)
                                                                                           cAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGt                                                           >  1:4321368/1‑150 (MQ=255)
                                                                                            aaCAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTc                                                          >  1:1267305/1‑150 (MQ=255)
                                                                                                  agcCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCcg                                                    >  1:2082763/1‑150 (MQ=255)
                                                                                                   gcCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCcgc                                                   <  2:4868188/150‑1 (MQ=255)
                                                                                                       gcgcCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCct                                               >  2:2478412/1‑150 (MQ=255)
                                                                                                          ccGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCAGGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTg                                            <  2:3607560/150‑1 (MQ=255)
                                                                                                           cGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGa                                           <  2:2261011/150‑1 (MQ=255)
                                                                                                              tCAAGCGCTGAGCCCGTAATGAATTAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACTTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCt                                        <  1:2635219/150‑1 (MQ=255)
                                                                                                                  gcgcTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGACGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTaa                                    <  2:5567295/150‑1 (MQ=255)
                                                                                                                  gcgcTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTaa                                    <  1:994196/150‑1 (MQ=255)
                                                                                                                        aGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGt                                      >  2:5154664/1‑142 (MQ=255)
                                                                                                                         gCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATAcc                             <  1:3560686/150‑1 (MQ=255)
                                                                                                                          cccGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCt                            >  2:5241110/1‑150 (MQ=255)
                                                                                                                           ccGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTa                           >  2:4500064/1‑150 (MQ=255)
                                                                                                                            cGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTAc                          >  2:410552/1‑150 (MQ=255)
                                                                                                                            cGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTAc                          <  1:5189881/150‑1 (MQ=255)
                                                                                                                               aatgaatgAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTa                       <  1:1040833/150‑1 (MQ=255)
                                                                                                                                      gagaGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGc                <  1:131624/150‑1 (MQ=255)
                                                                                                                                      aagagGCCCGTCACATGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGc                >  2:3032061/2‑150 (MQ=255)
                                                                                                                                      aagagGCCCGTCACATGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGc                >  1:4131027/2‑150 (MQ=255)
                                                                                                                                        gagGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCg              <  1:1023728/150‑1 (MQ=255)
                                                                                                                                          ggCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCt            <  1:710619/150‑1 (MQ=255)
                                                                                                                                           gCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGa                                           <  2:2565942/118‑1 (MQ=255)
                                                                                                                                           gCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGa                                           >  1:2565942/1‑118 (MQ=255)
                                                                                                                                           gCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTc           >  1:2774643/1‑150 (MQ=255)
                                                                                                                                           gCCCGTCACATGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTc           >  1:343063/1‑150 (MQ=255)
                                                                                                                                           gCCCGTCACATGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTc           >  1:3307026/1‑150 (MQ=255)
                                                                                                                                              cGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACg        >  1:4703116/1‑150 (MQ=255)
                                                                                                                                              cGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGACGCTCACg        >  2:2227843/1‑150 (MQ=255)
                                                                                                                                               gTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCActg       <  2:4373243/150‑3 (MQ=255)
                                                                                                                                               gTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACgg       <  2:80016/150‑1 (MQ=255)
                                                                                                                                               gTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACgg       <  2:3242412/150‑1 (MQ=255)
                                                                                                                                               gTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACgg       <  2:5473651/150‑1 (MQ=255)
                                                                                                                                               gTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACgg       >  2:3105210/1‑150 (MQ=255)
                                                                                                                                                  aTGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACGGGcc    >  2:1487349/1‑150 (MQ=255)
                                                                                                                                                    gtgACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACGGGCCAg  >  1:3980328/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATCCTGGCAAGGATCTTCCCCGACCATGAGGTGGTGATGATCCCGGGCCGCGAGCTGCTCTTGGGTGGCGGCAACATTCACTGCCTGACCCAACAGCAGCCGGCGCCGGTCAAGCGCTGAGCCCGTAATGAATGAGAGGCCCGTCATGTGACGGGCCTCTTTTTTGCGCGGGCGACGCCCGGAAGCCGAGCTGGCATACGTTTTGTATTGTTTTTCAAGCAGATAGCAAGATGGCACGGTCAGGCCGCGCCTCTGATCTGTAACAATACCTACGTAAATTAGCCGCTCACGGGCCAG  >  NC_002947/321595‑321892

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: