Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,882,009 T→A 6.2% intergenic (+26/+79) rpiA → / ← PP_5151 ribose‑5‑phosphate isomerase A/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,882,0090TA 6.2% 1056.1 / 10.0 350intergenic (+26/+79)rpiA/PP_5151ribose‑5‑phosphate isomerase A/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (140/178);  new base A (18/10);  total (160/190)
Fisher's exact test for biased strand distribution p-value = 4.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 373 total aligned reads displayed.

GGAGCTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCAT  >  NC_002947/5881860‑5882158
                                                                                                                                                     |                                                                                                                                                      
ggAGCTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCT                                                                                                                                                        <  1:2723752‑M1/150‑1 (MQ=255)
   gCTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACA                                                                                                                                                     <  2:268250‑M1/150‑1 (MQ=255)
     tGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAgg                                                                                                                                                   <  2:7855061‑M1/150‑3 (MQ=255)
      ggAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAgg                                                                                                                                                  >  2:5950060/1‑150 (MQ=255)
       gAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCGCCTACAAGGt                                                                                                                                                 >  2:6748201/1‑150 (MQ=255)
        taGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTATACCCGCTCCTACAAGGTa                                                                                                                                                <  2:1865008/149‑1 (MQ=255)
            gcAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTAcccc                                                                                                                                            >  2:2089523/1‑150 (MQ=255)
                  cAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGc                                                                                                                                      <  1:7608259/150‑1 (MQ=255)
                    aCGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCt                                                                                                                                   <  1:6522145/150‑1 (MQ=255)
                    aCGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCt                                                                                                                                   <  1:6189061/150‑1 (MQ=255)
                    aCGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCt                                                                                                                                   <  2:2795134/150‑1 (MQ=255)
                     cGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGATGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTg                                                                                                                                  >  1:8753335/1‑150 (MQ=255)
                              cgGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTgg                                                                                                                         <  1:710972/150‑1 (MQ=255)
                              cgGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTgg                                                                                                                         <  1:4577988/150‑1 (MQ=255)
                               ggCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGc                                                                                                                        <  2:2990265/150‑1 (MQ=255)
                                         ccAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgta                                                                                                              >  2:561478‑M2/1‑146 (MQ=255)
                                         ccAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgta                                                                                                              >  2:6944581‑M2/1‑146 (MQ=255)
                                         ccAACGGCCTGGTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgta                                                                                                              >  2:5999155‑M2/1‑146 (MQ=255)
                                          cAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtag                                                                                                             >  1:3450961‑M2/1‑145 (MQ=255)
                                              ggCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagc                                                                                                         >  1:7303299‑M2/1‑141 (MQ=255)
                                               gCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcg                                                                                                         <  2:3865931‑M1/150‑45 (MQ=255)
                                                ccTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgg                                                                                                        <  2:3079402‑M1/150‑46 (MQ=255)
                                                ccTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgg                                                                                                       >  2:4088879‑M2/1‑139 (MQ=255)
                                                ccTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgg                                                                                                       >  2:1667385‑M2/1‑139 (MQ=255)
                                                ccTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCGCTCCTACAACGTACCCCACCGGCCAGGAACCCTGGCGATCCCagtaggagcgg                                                                                                       <  2:5967639‑M2/150‑12 (MQ=255)
                                                 cTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcggg                                                                                                       <  2:5169028‑M1/150‑47 (MQ=255)
                                                 cTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcggg                                                                                                       <  2:4834285‑M1/150‑47 (MQ=255)
                                                 cTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcggg                                                                                                       <  2:7424862‑M1/150‑47 (MQ=255)
                                                 cTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcggg                                                                                                      >  2:3873732‑M2/1‑138 (MQ=255)
                                                 cTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcggg                                                                                                      >  2:7241386‑M2/1‑138 (MQ=255)
                                                  tGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggt                                                                                                     >  1:8254155‑M2/1‑137 (MQ=255)
                                                  tGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACACGCGCCAACAAAGTAACCCAACGGCCAGGGAAACGGGCGATCCCGGAAGGGGCgggg                                                                                                     >  2:147309/1‑149 (MQ=255)
                                                   gTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtt                                                                                                    >  2:6095654‑M2/1‑136 (MQ=255)
                                                   gTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgta                                                                                                              >  2:6486355‑M2/1‑136 (MQ=255)
                                                       gcggcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGGGAACCCCCGCCAACAAGGTACCCCCCCCGCCGGGGACCCGGGCGATCCCcgtgggagcggggttacc                                                                                                <  1:6399728‑M2/150‑19 (MQ=255)
                                                        cggcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttaccc                                                                                               >  2:2186723‑M2/1‑131 (MQ=255)
                                                        cggcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttaccc                                                                                               >  1:3824502‑M2/1‑131 (MQ=255)
                                                               gccggcGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCCCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACa                                                                                        >  1:3279813/1‑150 (MQ=255)
                                                                          tgctgctgctgGGCACCGCTGACGGCGTGAAGAGCCTGAAGGCCGATTAAACCCCTTCGCGGGTAAACCCGCACCTACACGGTACCCCACGGGCCTGGAACGGTGGCGATCCCTGTCGGAGGGGTTTTACCCGCGAACAGGCCCTGACAg                                                                             <  1:5194513/150‑1 (MQ=255)
                                                                          tgctgctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAg                                                                             >  2:2609721/1‑150 (MQ=255)
                                                                            cggctgctgGGCACCGCTGACGGCTTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGGGGGTAAACCCGCTCCTACACGGTACCCCACGGGCCTGTAACCGTGGCGATCCCTGTTGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGc                                                                           <  1:182363/148‑1 (MQ=255)
                                                                             tgctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGc                                                                                                                                      >  2:1245360/1‑91 (MQ=255)
                                                                              gctggtgGGCACCGCTGATGGGGTTAAGAGCCTGATGGCCGATTAAAACCCTTCGGGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGACCGGTGGCGATCCCTGTTGGAGGGGGTTTACCCGCGAACAGGCCCTGACAGGCGa                                                                         <  1:7449047/150‑1 (MQ=255)
                                                                               gtgctgGGCACCGCTTATGGCGTTAAGCGCCTGAAGGCCGATAAAAACCCTCCGGGGTTAAACCCGCTCCTAAACGTTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTTGGCGGGGTTTTACCCGCGAACAGGCCCTGACAGGCGAt                                                                        <  1:3111098/149‑1 (MQ=255)
                                                                               ctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAt                                                                        <  2:3563784/150‑1 (MQ=255)
                                                                                 ggtgGGCACCGCTGACGGCGTTAAGCGCCTGACGGCCGTGTAAACCCCTTCGGGGTAAAACCCGCTCCTACAAGGTACCCCACCGGCCTGTAACCGTGGCGATCCCTGTAGGCGGGGTTTTACCCGCGAACAGGCCCTGACAGGCGAtca                                                                      <  1:1123264/148‑1 (MQ=255)
                                                                                   tgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagt                                                                    <  2:3770947/150‑1 (MQ=255)
                                                                                     ggCACCGCTGACGGCTTGAAGCGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAACCCCGCTCCTACAAGGTACCCCACCGGCCTGGAACGGTGGCGACCCCGGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagtca                                                                  <  1:5546941/150‑1 (MQ=38)
                                                                                     ggCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGTGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagtca                                                                  <  1:591689/150‑1 (MQ=255)
                                                                                     ggCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagtca                                                                  <  1:7188086/150‑1 (MQ=255)
                                                                                     ggCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagtca                                                                  <  1:5760586/150‑1 (MQ=255)
                                                                                      gCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCGCCTAAAAGGTACCCCCCCGGCCGGGGACCGGGGCGATCCCGGGTGGAGCCGGTTTACCCGCCAACAGGCCCTGACCGGCGAtcagtcag                                                                 >  2:2406149/1‑150 (MQ=255)
                                                                                          cGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCAGACAGGCGATCCGTCAGTCgg                                                             >  2:6571635/1‑150 (MQ=255)
                                                                                            cTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGACCCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGgct                                                           <  1:1253141/150‑1 (MQ=255)
                                                                                               aaGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGgctgct                                                        <  2:6868703/150‑1 (MQ=255)
                                                                                                   cGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGaa                                                    <  1:2089523/150‑1 (MQ=255)
                                                                                                   cGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCATGGAACCGTGGCTATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGaa                                                    <  1:1964743/150‑1 (MQ=255)
                                                                                                    gTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGc                                                                                                                        <  1:4475050/81‑1 (MQ=255)
                                                                                                         gagCCTGATGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCACCTACAAGTTACCCCACCGGCCTGTAACGGTGGCGATCCCTGTAGGTGGGGTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATa                                              <  1:7758945/150‑1 (MQ=255)
                                                                                                          agCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAg                                             <  2:4597568/150‑1 (MQ=255)
                                                                                                                aaGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAgg                                       <  2:4973056/150‑1 (MQ=255)
                                                                                                                  ggCCGAGTAAAACCCTTCGCGGGTTAACCCGCTCCAACAAGGTACCCCACCGGCCGGGAGACGTGGGGATGCCGGTAGGAGCGGGTTTACCCGCGAACAGGCCCAGACAGGCCCTCAGTCAGTCGGCGGCGTGAACACATAGAACAGGtt                                     >  2:4177897/1‑150 (MQ=255)
                                                                                                                      gAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGc                                 <  2:3452240/150‑1 (MQ=255)
                                                                                                                       aGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCt                                <  1:4214536/150‑1 (MQ=255)
                                                                                                                          aaaaCCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCAc                             <  2:622803/150‑1 (MQ=255)
                                                                                                                          aaaaCCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCAc                             <  2:5674159/150‑1 (MQ=255)
                                                                                                                              cccTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGAc                         <  1:3079007/150‑1 (MQ=255)
                                                                                                                               ccTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACAAAGAACAGGTTCGGCTCACTGAcc                        >  2:8779127/1‑150 (MQ=255)
                                                                                                                                cTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCa                       <  1:4929189/150‑1 (MQ=255)
                                                                                                                                 ttCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgg                                                                                                       <  1:4039326‑M2/69‑12 (MQ=255)
                                                                                                                                    gcgGGTAACCCCGCTCCTACAgggatcgccacggttccaggcggggggggtaccttgttggagggggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgc                    <  2:580631‑M1/150‑130 (MQ=255)
                                                                                                                                      gggTAAACCCGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccg                  >  2:6131255‑M1/1‑19 (MQ=255)
                                                                                                                                      gggTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTaaa                 >  1:8052697/1‑150 (MQ=255)
                                                                                                                                      gggTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTaaa                 <  1:7874814/150‑1 (MQ=255)
                                                                                                                                       ggTAAACCCGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgc                 >  1:3587013‑M1/1‑18 (MQ=255)
                                                                                                                                        gTAAACCCGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcg                >  2:1303524‑M1/1‑17 (MQ=255)
                                                                                                                                        gTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATg               <  2:545418/150‑1 (MQ=255)
                                                                                                                                        gTAAACCCGCACCTACCgggatcgccagggtccctggccggtggggtaccttgttggagcggttttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcg                <  2:1065162‑M1/150‑134 (MQ=255)
                                                                                                                                            aCCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTcc           <  2:116001/150‑1 (MQ=255)
                                                                                                                                             cccGCTCCTACAgggatcgccacggttccaggccggtggggtatcttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacag           >  2:3726847‑M1/1‑12 (MQ=255)
                                                                                                                                             cccGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacag           >  2:5651466‑M1/1‑12 (MQ=255)
                                                                                                                                             cccGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacag          <  1:3902044‑M2/150‑105 (MQ=255)
                                                                                                                                              ccGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacagg          >  1:4338332‑M1/1‑11 (MQ=255)
                                                                                                                                              ccGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaa                                                                                           >  2:4039326‑M1/1‑11 (MQ=255)
                                                                                                                                              ccGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacagg         <  2:5062959‑M2/150‑106 (MQ=255)
                                                                                                                                              ccGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTcccc         >  2:4440153/1‑150 (MQ=255)
                                                                                                                                               cGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggc         <  2:4554304‑M1/150‑141 (MQ=255)
                                                                                                                                               cGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCg        >  2:1228556/1‑150 (MQ=255)
                                                                                                                                               cGCCCCTACTgggatcgccacggttccaggccggtggggtaccttttaggagcggttttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggc         <  2:2529466‑M1/150‑141 (MQ=255)
                                                                                                                                                gCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcggagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggcc        >  1:8524719‑M1/1‑9 (MQ=255)
                                                                                                                                                 cTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccg      <  2:8319173‑M2/150‑109 (MQ=255)
                                                                                                                                                  tCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgt     <  1:6259984‑M2/150‑110 (MQ=255)
                                                                                                                                                   ccTACCgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgtt     <  2:7673499‑M1/150‑145 (MQ=255)
                                                                                                                                                    cTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgttg    >  2:4354149‑M1/1‑5 (MQ=255)
                                                                                                                                                    cTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgttg    >  2:1565748‑M1/1‑5 (MQ=255)
                                                                                                                                                    cTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTca   >  2:7140297/1‑150 (MQ=255)
                                                                                                                                                     tACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgttgg  <  1:1981‑M2/150‑113 (MQ=255)
                                                                                                                                                     tACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTcat  >  1:6619997/1‑150 (MQ=255)
                                                                                                                                                     tACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTcat  >  2:8451332/1‑150 (MQ=255)
                                                                                                                                                     tACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGGAAATGGTCCCCGCTTcat  >  1:2190587/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GGAGCTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCAT  >  NC_002947/5881860‑5882158

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