Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,508,546 A→T 5.8% intergenic (+35/‑164) PP_4839 → / → cycA iron‑regulated membrane protein/D‑alanine/beta‑alanine/D‑serine/glycine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,508,5460AT 5.8% 263.2 / 2.3 83intergenic (+35/‑164)PP_4839/cycAiron‑regulated membrane protein/D‑alanine/beta‑alanine/D‑serine/glycine permease
Reads supporting (aligned to +/- strand):  ref base A (38/40);  new base T (2/3);  total (40/43)
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

CCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTG  >  NC_002947/5508399‑5508691
                                                                                                                                                   |                                                                                                                                                 
ccGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTtaa                                                                                                                                                 >  1:4151630/1‑150 (MQ=255)
 cGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAc                                                                                                                                                >  2:509665/1‑150 (MQ=255)
 cGTGATGCTGGTGCTAGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTtaaa                                                                                                                                                >  1:1020724/1‑149 (MQ=255)
  gTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAcc                                                                                                                                               >  1:594751/1‑150 (MQ=255)
  gTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAcc                                                                                                                                               <  2:2012468/150‑1 (MQ=255)
    gATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAcacc                                                                                                                                             >  1:528493/1‑147 (MQ=255)
      tgctggtgctggGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCt                                                                                                                                           >  2:2968333/1‑150 (MQ=255)
          ggtgctggGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTATGCTGGCCTCTTCGCGGGTTTACCCGcgaaga                                                                                                                                       <  2:148725/150‑6 (MQ=255)
          ggtgctggGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTATGCTGGCCTCTTCGCGGGTTTACCCGcgaaga                                                                                                                                       <  2:148719/150‑6 (MQ=255)
           gtgctggGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTTTACCCGcgaagag                                                                                                                                      <  2:1745953/150‑7 (MQ=255)
                gggggTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAACCCCACCGACGAggg                                                                                                                                 >  1:1268362/1‑150 (MQ=255)
                 ggggAGCTGGTCCCGCTGGTGCGGGTGGCGATGGTGCGGCGGGGGGGGGTGGAGAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCGTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGa                                                                                                                                <  1:4955532/150‑1 (MQ=37)
                       aTGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGt                                                                                                                          >  1:4781079/1‑150 (MQ=255)
                        tGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGtt                                                                                                                         <  2:3138231/150‑1 (MQ=255)
                        cgTTCCCGCTGGGGGGGGTGTCGATGGTGGTGATGGGGGGGGGGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGtt                                                                                                                         <  2:2810216/149‑1 (MQ=255)
                         gTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCCCTCCTACGAGGGACCGCGTTg                                                                                                                        >  2:5062078/1‑150 (MQ=255)
                          ttCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGt                                                                                                                       >  2:4841347/1‑150 (MQ=255)
                          ttCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTGGGTGGTGAGGCGCCCCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGt                                                                                                                       >  2:4841176/1‑150 (MQ=255)
                           tCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTa                                                                                                                      <  2:1735763/150‑1 (MQ=255)
                           tCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTa                                                                                                                      <  2:4477705/150‑1 (MQ=255)
                           tCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCATTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCACCTACGAGGGACCGCGGTGTa                                                                                                                      >  2:2527765/1‑150 (MQ=255)
                            cccGCTGGTGGGGGTGTCGCTGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAg                                                                                                                     <  2:3006042/150‑1 (MQ=255)
                            cccGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCGAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGCCCGCGTTGTAg                                                                                                                     >  2:3514431/1‑150 (MQ=255)
                            cccGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAg                                                                                                                     <  2:1149000/150‑1 (MQ=255)
                            cccGCTGGTGGGGGTGTCGATGGTGGTGAGGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGGTGTAg                                                                                                                     >  1:1745639/1‑150 (MQ=255)
                             ccGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAgg                                                                                                                    <  2:4781079/150‑1 (MQ=255)
                             ccGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAgg                                                                                                                    >  2:1172966/1‑150 (MQ=255)
                             ccGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAgg                                                                                                                    >  2:1228529/1‑150 (MQ=255)
                                cTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGc                                                                                                                 <  2:3749097/150‑1 (MQ=255)
                                cTGGGGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCGTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGc                                                                                                                 <  2:4389528/150‑1 (MQ=255)
                                 gggtggGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCg                                                                                                                <  2:1185049/149‑1 (MQ=255)
                                  ggtggGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCgg                                                                                                               <  2:2874377/150‑1 (MQ=255)
                                   gtggGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCggg                                                                                                              <  1:1228529/150‑1 (MQ=255)
                                   gtggGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCggg                                                                                                              <  1:1172966/150‑1 (MQ=255)
                                      ggggTGTCGAGGGTGGTGATGGGGGGGGTGGAGAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGttt                                                                                                           <  2:3764352/150‑1 (MQ=255)
                                         gtgtCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTAcc                                                                                                        >  2:873848/1‑150 (MQ=255)
                                           gtCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCggg                                                                                                              >  2:1647468/1‑142 (MQ=255)
                                           gtCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCggg                                                                                                              <  1:1647468/142‑1 (MQ=255)
                                              gATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGa                                                                                                   >  2:3884017/1‑150 (MQ=255)
                                              gATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCGCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGa                                                                                                   >  2:2014229/1‑150 (MQ=255)
                                                  gtggtgATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTTAAACCCCACCTACCAGGGACCCCGGTGTAGGAGCGGGGTTAACCCCCAagag                                                                                               >  1:4748935/1‑150 (MQ=255)
                                                   tggtgATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCgg                                                                                                               <  1:3039749/133‑1 (MQ=255)
                                                   tggtgATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCgg                                                                                                               >  2:3039749/1‑133 (MQ=255)
                                                   tggtgATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGGTGTAGGAGCGGGTTTACCCGCGAAGAgg                                                                                              >  2:3547662/1‑150 (MQ=255)
                                                    ggtgCTGGGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTGGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTCCGAGGGACCGCGTTGTCGGAGCGGGTTTACCCGCGAAGAGGc                                                                                             <  2:2218928/150‑1 (MQ=38)
                                                           tgGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACGGGGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACa                                                                                      <  2:3862137/150‑1 (MQ=255)
                                                             gggggTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCGGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGTACAgg                                                                                    <  2:5080571/150‑1 (MQ=255)
                                                                  tGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGaaa                                                                               >  2:2484493/1‑150 (MQ=255)
                                                                    gATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTTAAAACCCTCCTAAGAGGGGCCGCGGTGTAGGAGGGGGTTTTACCCCGGAGAAGCCCGGACAGGCGCAATg                                                                             >  1:604322/1‑150 (MQ=255)
                                                                    gATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTACTACGAGGGACCGCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATg                                                                             >  2:1453561/1‑150 (MQ=255)
                                                                      ttgTTGGGGGGTGCGGCGCCGCCGGGCCTGGGCCAGCGCTTGCGCGTTGTGTTGCCGGTGCGGGCCTCTCCGGGGGAAACCCCGCCCCCCCGCGGGACCGGTTTGTCGGCGGGGTTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCt                                                                           <  2:1472408/148‑1 (MQ=255)
                                                                      tcgTTGGGTGGTGCGGCGCCGCCGGGTCTTGGCCCGCGCTTGAGCCTTGTGTGGCCTGTGCGGGCCTCTTCGGGGGAAAACCCGCGCCTCCGCGGGCCCGGGTTGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCt                                                                           <  2:1472393/148‑1 (MQ=37)
                                                                        gTTGGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGCCATCTTCGCGGGAAAACCCGCACCTACGAGGGACCGCGTTGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGa                                                                         <  2:4683645/150‑1 (MQ=38)
                                                                          ttGGTGGTGAGGCGCCGCCGGGTGTGGGCCTGCGCTTGCGACTGGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGGTCCTACTAGGGACCGCGTTGTAGGAGGGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGc                                                                       <  2:5205909/150‑1 (MQ=37)
                                                                           gggtggtGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGcc                                                                      <  2:3591225/149‑1 (MQ=255)
                                                                               ggtgTGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGTAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCCGAGCCGaaa                                                                  <  2:3127984/150‑1 (MQ=255)
                                                                                gtgAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAACCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGaaaa                                                                 >  2:2899722/1‑150 (MQ=255)
                                                                                   aGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGc                                                                                                                                            <  1:77223/72‑1 (MQ=255)
                                                                                   aGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGc                                                                                                                                            <  1:77224/72‑1 (MQ=255)
                                                                                   aGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGc                                                                                                                                            >  2:77224/1‑72 (MQ=255)
                                                                                   aGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGc                                                                                                                                            >  2:77223/1‑72 (MQ=255)
                                                                                    ggCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAACCGCTCCAAAGAGGGACCGCGTTGTAGGAGCGGGTTTACCCCCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAAAGCa                                                             >  1:5335824/1‑150 (MQ=38)
                                                                                     gcgcCGCCGGGTCTTGGCCTGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCCCTCCTAACAGGGACCGCGGTGTAGGAGCGGGGTTACCCCCGAAGAGGCCGGGACAGGCGCAATGCTGAGCCGAAAATTCag                                                            >  2:2322755/1‑150 (MQ=37)
                                                                                     gcgcCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGCGTTGTAGGAGCGGGTTTAACCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAAAGCag                                                            >  2:4306887/1‑150 (MQ=255)
                                                                                      cgccgcGGGGTCTGGGCCCGCGCTTGTGCGTGGTGTGGCCTGTGCGGGCCGCTTCGCGGTTAAACCCGCTCCTACGTGGGACGGCTTTTTAGGAGGGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCaga                                                           <  1:3547662/150‑1 (MQ=25)
                                                                                               gTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCCCACCTACCCGGGGCCGCGTTGTTGGAGCGGGTTTTACCGCGCAGAGGCCCGTAAAGGCGCAAAGCTGGGCCGAAAAAGCAGAGTCGATCTg                                                  >  1:4392240/1‑150 (MQ=25)
                                                                                                        agggcTTGCGCCTTGGGTTGCCGGGGCGGGCCCCTGCGCGGGTAACCCCGCTCCTCCGCGGGCCCGCGGTGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGcccc                                         <  2:4151630/147‑1 (MQ=25)
                                                                                                        aGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTAAGAGGGACCCCGGTGTTGGATCGGGTTTAACCGCGAAGAGGCCGGTAAAGGCGAAATTCGGAGCCGAAAATGCCGAGTCGATGTGGCGCGcccc                                         >  2:4556447/1‑150 (MQ=25)
                                                                                                              tGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCCATCTg                                   >  2:3266694/1‑150 (MQ=255)
                                                                                                                      gtgtTGCCTGTGCTGGCCTCTCCGCGGTTAAACCCGCTCCTACGAGGGACCGCGTTTTAGTAGCGGTTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAAcgcgc                           <  2:2298548/150‑1 (MQ=38)
                                                                                                                      gtgtTGCCTGGGCTGGCCTCTTCGCGGGTAACCCCGCTCCTACGCGGGACCGCGTTGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAAcgcgc                           <  2:1654129/150‑1 (MQ=255)
                                                                                                                       tgtTGCCTGTGCTGGCCTCTTCGCGGGTAAAACCGCTCCTACGAGGGACCGCGTTGGAGGAGCGGGTTTAACCGCGAAGAGGCCGGTTCAGGCGAAATGCGGTGCCGAAAATGCCGAGTCGATCTGTCGCGCCCCGATCTGGGAcgcgcg                          >  1:3621385/1‑150 (MQ=37)
                                                                                                                            ccggtgCTGGCCACTCCGCGGGTAACCCCGCTCCTACGCGGGCCGGGGTTTTTGGAGGGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTta                     <  2:1948069/147‑1 (MQ=255)
                                                                                                                            ccTGGGGTGGCCTCTTCGCGGGTAACCCCGCGCCTCCGAGGGACCGCGTGGTAGGCGGGGTTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTta                     <  2:1948022/150‑1 (MQ=255)
                                                                                                                                   tGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCGGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGAGCGGGACAGCGCGTGTtaggctagc              >  2:4205209/1‑150 (MQ=37)
                                                                                                                                     gCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCg            <  2:4450024/150‑1 (MQ=255)
                                                                                                                                       ctctTCGCGGGTAAACCCCCTCCTACCAGGGACCGCGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCGGAGCCGAAAATGCAGAGTCGAACTGTCGCGCCCCGAACTGGAACCCCCCTGTTAGCCAAGCCCGc          >  2:4657983/1‑150 (MQ=37)
                                                                                                                                         ctTCGCGGTTAACCCCGCGCCCACGCGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTa        <  2:1020724/150‑1 (MQ=39)
                                                                                                                                         ctTCGCGGGTAAACCCGCTCCTCCGAGGGACCGCGTTGTAGGAGGGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTa        <  2:1268362/150‑1 (MQ=255)
                                                                                                                                          ttCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTAc       >  2:5383837/1‑150 (MQ=255)
                                                                                                                                             ggggtaaaaCCCGCGCTTACGTGGCACCGGTTTGGTGGTGGGGTTTCCCCCGCGAAGAGGCCGGTACGGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACctc    <  1:509665/144‑1 (MQ=255)
                                                                                                                                             gcgGGTAAACCCGCTCCTACGAGGGACCGCGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCAGTCGCGCCCCGATCTGGAACGCGCGTGGTAGCCTAGCCGGCGACctc    >  1:1985908/1‑150 (MQ=255)
                                                                                                                                              cgGGTAAACCCGCCCCTCCGAGGGACGGGGTTGTTGGTGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACctct   <  2:270677/150‑1 (MQ=255)
                                                                                                                                               gggTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGtt                      >  1:2886959/1‑130 (MQ=255)
                                                                                                                                               gggTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGtt                      <  2:2886959/130‑1 (MQ=255)
                                                                                                                                               gggTAAACCCGCTCCTACGAGGGACCGCGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCATAGCCGGCTACCTCTg  >  1:4017265/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTG  >  NC_002947/5508399‑5508691

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