Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,467,390 T→G 5.0% intergenic (+16/+18) PP_4798 → / ← PP_4799 membrane‑bound lytic murein transglycosylase/muramoyltetrapeptide carboxypeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,467,3900TG 5.0% 932.9 / 25.8 270intergenic (+16/+18)PP_4798/PP_4799membrane‑bound lytic murein transglycosylase/muramoyltetrapeptide carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base T (133/123);  new base G (10/4);  total (143/127)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 273 total aligned reads displayed.

GATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGCACACC  >  NC_002947/5467243‑5467539
                                                                                                                                                   |                                                                                                                                                     
gatgCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTgc                                                                                                                                                     <  1:1826549/150‑1 (MQ=255)
  tgCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTgcgc                                                                                                                                                   <  2:2320790/150‑1 (MQ=255)
  tgCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTgcgc                                                                                                                                                   <  1:60349/150‑1 (MQ=255)
    ccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCgg                                                                                                                                                 >  2:2335772/5‑150 (MQ=255)
       gTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGccc                                                                                                                                              <  2:1445191/150‑1 (MQ=255)
         aaTGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCg                                                                                                                                            <  2:5557591/150‑1 (MQ=255)
          aTGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCgg                                                                                                                                           <  1:3830947/150‑1 (MQ=255)
           tGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGt                                                                                                                                          >  1:2202108/1‑150 (MQ=255)
              cGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTatg                                                                                                                                       <  1:4927862/150‑3 (MQ=255)
                  ggCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCa                                                                                                                                   <  1:315358/150‑1 (MQ=255)
                  ggCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCa                                                                                                                                   <  1:3841066/150‑1 (MQ=255)
                   gCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCAc                                                                                                                                  <  2:5159653/150‑1 (MQ=255)
                   gCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCAc                                                                                                                                  <  1:5333427/150‑1 (MQ=255)
                    cATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTCATCTCACg                                                                                                                                 >  2:2084832/1‑150 (MQ=255)
                     aTTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGc                                                                                                                                >  2:1650687/1‑150 (MQ=255)
                     aTTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGc                                                                                                                                >  2:701212/1‑150 (MQ=255)
                      ttATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTCATCTCACGcc                                                                                                                               >  1:3944957/1‑150 (MQ=255)
                        aTCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCt                                                                                                                             <  2:3195243/150‑1 (MQ=255)
                           ggggCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGaa                                                                                                                          >  1:4912448/1‑150 (MQ=255)
                           ggggCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGaa                                                                                                                          <  2:2170218/150‑1 (MQ=255)
                           ggggCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGaa                                                                                                                          <  2:2097526/150‑1 (MQ=255)
                           ggggCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGaa                                                                                                                          <  1:2240053/150‑1 (MQ=255)
                           ggggCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGaa                                                                                                                          <  2:627210/150‑1 (MQ=255)
                            gggCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAAc                                                                                                                         >  1:4110073/1‑150 (MQ=255)
                            gggCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAAc                                                                                                                         >  1:4168089/1‑150 (MQ=255)
                              gCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAg                                                                                                                       >  1:2353767/1‑150 (MQ=255)
                                 aaTACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCTACACGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATc                                                                                                                    <  2:2025233/150‑1 (MQ=255)
                                 aaTACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATc                                                                                                                    <  1:453385/150‑1 (MQ=255)
                                      ccGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGTCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCt                                                                                                               >  1:2913349/1‑150 (MQ=255)
                                        gCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTcc                                                                                                             <  2:4187440/150‑1 (MQ=255)
                                         cAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCa                                                                                                            <  1:641633/150‑1 (MQ=255)
                                         cAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCa                                                                                                            >  2:1154472/1‑150 (MQ=255)
                                            ggCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCa                                                                                                         >  2:3288386/1‑150 (MQ=255)
                                              cAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCAcc                                                                                                       >  1:3085470/1‑150 (MQ=255)
                                                 tCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACACACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAgt                                                                                                    >  2:1901847/1‑150 (MQ=255)
                                                      aaTGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTg                                                                                               <  2:266921/150‑1 (MQ=255)
                                                        tGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCt                                                                                             <  2:3674246/150‑1 (MQ=255)
                                                           ccAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGaa                                                                                          <  1:1060011/150‑1 (MQ=255)
                                                           ccAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGaa                                                                                          <  2:4401397/150‑1 (MQ=255)
                                                            cAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAac                                                                                         >  2:2496377/1‑150 (MQ=255)
                                                                      tGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGc                                                                               <  2:4986639/150‑1 (MQ=255)
                                                                             gCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCAccc                                                                        <  1:2839785/150‑1 (MQ=255)
                                                                             gCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCAccc                                                                        <  1:3722213/150‑1 (MQ=255)
                                                                              ccGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCGGTGTTTGCTGAACACTGTCCAGCCTCACCCg                                                                       >  1:5393492/1‑150 (MQ=255)
                                                                               cggcggATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGa                                                                      <  2:4419507/150‑1 (MQ=255)
                                                                               cggcggATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGa                                                                      >  2:3417267/1‑150 (MQ=255)
                                                                               cggcggATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGa                                                                      <  1:4347268/150‑1 (MQ=255)
                                                                               cggcggATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGACCAGCCTCACCCGa                                                                      >  2:2823806/1‑150 (MQ=255)
                                                                                       ggCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAg                                                              <  1:2704219/150‑1 (MQ=255)
                                                                                         cTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAggg                                                            <  1:3305517/150‑1 (MQ=255)
                                                                                         cTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAggg                                                            >  1:307306/1‑150 (MQ=255)
                                                                                          tACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAgggg                                                           >  1:170529/1‑150 (MQ=255)
                                                                                            cccGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCa                                                         >  2:294139/1‑150 (MQ=255)
                                                                                            cccGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCa                                                         >  2:5417547/1‑150 (MQ=255)
                                                                                             ccGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCaa                                                        >  1:5566875/1‑150 (MQ=255)
                                                                                             ccGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCaa                                                        <  1:1757941/150‑1 (MQ=255)
                                                                                              cGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAg                                                       <  2:882822/150‑1 (MQ=255)
                                                                                               gACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGc                                                      >  2:4258308/1‑150 (MQ=255)
                                                                                                  ccACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATa                                                   >  2:2376475/1‑150 (MQ=255)
                                                                                                  ccACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATa                                                   <  2:5076908/150‑1 (MQ=255)
                                                                                                  ccACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATa                                                   >  1:814241/1‑150 (MQ=255)
                                                                                                    acaAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACa                                                 >  2:4686243/1‑150 (MQ=255)
                                                                                                         ctgctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCgg                                            <  1:2011716/150‑1 (MQ=255)
                                                                                                         ctgctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCgg                                            <  1:5098933/150‑1 (MQ=255)
                                                                                                         ctgctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCgg                                            <  2:2841317/150‑1 (MQ=255)
                                                                                                          tgctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCggg                                           <  1:3081456/150‑1 (MQ=255)
                                                                                                            ctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTc                                         >  1:1365793/1‑150 (MQ=255)
                                                                                                             ttGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTca                                        <  1:4687418/150‑1 (MQ=255)
                                                                                                               gagaGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTcaca                                      <  2:4111816/150‑1 (MQ=255)
                                                                                                                agagCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAg                                     >  1:4461466/1‑150 (MQ=255)
                                                                                                                   gCCTGCGCCAGCATTAAGCATTCCGGGGCCGCTGCTCGGCACCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTTTTTTCTGAACACTGTCCAGCCTCAACCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGc                                  <  1:3966032/150‑1 (MQ=255)
                                                                                                                   gCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGTGTCACAGTGc                                  <  2:4577101/150‑1 (MQ=255)
                                                                                                                      tGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCg                               >  2:4588128/1‑150 (MQ=255)
                                                                                                                      tGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCg                               >  1:449860/1‑150 (MQ=255)
                                                                                                                       gcgcCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGc                              <  2:4792068/150‑1 (MQ=255)
                                                                                                                       gcgcCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGc                              <  1:4957041/150‑1 (MQ=255)
                                                                                                                         gcCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTg                            >  2:4039682/1‑150 (MQ=255)
                                                                                                                           cagcagTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTgcg                          >  2:1220541/1‑150 (MQ=255)
                                                                                                                             gcagTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGcc                        <  1:48679/150‑1 (MQ=255)
                                                                                                                               agTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCaa                      <  2:3065729/150‑1 (MQ=255)
                                                                                                                                  gagattaccGGGGCCGCGCAGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAAccc                   >  1:3628472/8‑150 (MQ=255)
                                                                                                                                   aGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCg                  >  2:1310829/1‑150 (MQ=255)
                                                                                                                                    gattaccGGGGCCGCGCAGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGc                 >  2:4768223/6‑150 (MQ=255)
                                                                                                                                    gCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGc                 <  1:3650305/150‑1 (MQ=255)
                                                                                                                                    gCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGc                 <  1:5605781/150‑1 (MQ=255)
                                                                                                                                      aTTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCa               <  2:4343808/150‑1 (MQ=255)
                                                                                                                                        tCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGc             <  1:5217878/150‑1 (MQ=255)
                                                                                                                                        accGGGGCCGCGCAGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGc             >  1:901394/2‑150 (MQ=255)
                                                                                                                                         ccGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCAAACCCGCCAGCa            >  2:1836061/1‑150 (MQ=255)
                                                                                                                                         ccGGGGCCGCGCAGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCa            >  1:5000791/1‑150 (MQ=255)
                                                                                                                                         ccGGGGCCGCGCAGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCa            >  1:2476267/1‑150 (MQ=255)
                                                                                                                                             ggCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCgg        <  1:2335708/150‑1 (MQ=255)
                                                                                                                                             ggCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCgg        <  1:4486893/150‑1 (MQ=255)
                                                                                                                                               ccGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGca      >  1:141087/1‑150 (MQ=255)
                                                                                                                                               ccGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGca      >  1:2267266/1‑150 (MQ=255)
                                                                                                                                                cGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGcac     <  1:1154472/150‑1 (MQ=255)
                                                                                                                                                cGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGcac     <  2:5590448/150‑1 (MQ=255)
                                                                                                                                                 gCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGcaca    >  2:5151034/1‑150 (MQ=255)
                                                                                                                                                 gCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGcaca    >  2:1487997/1‑150 (MQ=255)
                                                                                                                                                 gCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGcaca    >  2:1787513/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGCACACC  >  NC_002947/5467243‑5467539

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: