Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 429,115 C→G 6.9% G137A (GGC→GCC)  PP_0352 ← ECF family RNA polymerase sigma‑70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947429,1150CG 6.9% 225.7 / 7.4 73G137A (GGC→GCC) PP_0352ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base C (33/35);  new base G (3/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-01

ATCGCGGCTTCCAGGACCCGGGACGAAACCGGGAGGTTGCTCACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGCTTCAGGGATCAGTGCCA  >  NC_002947/428966‑429259
                                                                                                                                                     |                                                                                                                                                  
aTCGCGGCTTCAAGGACCCGGGACGGAACCGGGAGGTTGCTCACGTCGGCTCCCCGTACAAGGCCAGATCTGACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGc                                                                                                                                                    <  1:1380729/150‑1 (MQ=255)
      gCTTCCAGGACCCGGGACGAAACCGGGAGGCTGCTCACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGc                                                                                                                                              <  1:1005656/150‑1 (MQ=255)
            aGGACCCGGGACGAAACCGGGAGGTTGCTCACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCa                                                                                                                                        >  1:1369676/1‑150 (MQ=255)
               aCCCGGGACGAAACCGGGAGGTTGCTCACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGAc                                                                                                                                     <  2:898439/150‑1 (MQ=255)
               aCCCGGGACGAAACCGGGAGGTTGCTCACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGAc                                                                                                                                     <  1:846612/150‑1 (MQ=255)
                      aCGAAACCGGGAGGTTGCTCACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCtgtg                                                                                                                              >  2:144429/1‑150 (MQ=255)
                        gAACCCGGGAGGGTGCGCACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCa                                                                                                                            >  2:379741/1‑150 (MQ=255)
                         aaaCCGGGAGGTTGCTCACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAg                                                                                                                           >  2:545508/1‑150 (MQ=255)
                          aaCCGGGAGGTTGCTCACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACCGCTGTGCAGa                                                                                                                          >  2:243376/1‑150 (MQ=255)
                           aCCGGGAGGTTGCTCACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCAt                                                                                                                                                <  2:1314371/127‑1 (MQ=255)
                           aCCGGGAGGTTGCTCACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCAt                                                                                                                                                >  1:1314371/1‑127 (MQ=255)
                           aCCGGGAGGTTGCTCACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGaa                                                                                                                         >  1:703105/1‑150 (MQ=255)
                             cGGGAGGTTGCTCACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAAc                                                                                                                       <  1:144429/150‑1 (MQ=255)
                               ggAGGTTGCTCACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAAcgc                                                                                                                     >  1:509886/1‑150 (MQ=255)
                                         cACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATACCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGcc                                                                                                           >  2:447995/1‑150 (MQ=255)
                                          aCGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGACCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCt                                                                                                          >  1:307071/1‑150 (MQ=255)
                                                  tCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGAc                                                                                                  >  2:868181/1‑150 (MQ=255)
                                                  tCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCACGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGAc                                                                                                  <  2:703105/150‑1 (MQ=255)
                                                     ccGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGc                                                                                               >  1:184305/1‑150 (MQ=255)
                                                        tACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTg                                                                                            >  2:50855/1‑150 (MQ=255)
                                                         aCAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCa                                                                                                                                        <  1:1332878/105‑1 (MQ=255)
                                                         aCAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCa                                                                                                                                        >  2:1332878/1‑105 (MQ=255)
                                                         aCAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCt                                                                                                      >  2:1086669/1‑139 (MQ=255)
                                                         aCAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCt                                                                                                      <  1:1086669/139‑1 (MQ=255)
                                                         aCAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGa                                                                                                   <  1:1259524/142‑1 (MQ=255)
                                                         aCAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGa                                                                                                   >  2:1259524/1‑142 (MQ=255)
                                                           aaGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑Cg                                                                                       >  1:265911/1‑150 (MQ=255)
                                                                             cgcATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGCGCAGGCGGTCGATGGCt                                                                     <  2:950467/150‑1 (MQ=255)
                                                                                aTGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTg                                                                  >  1:699725/1‑150 (MQ=255)
                                                                                 tGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGa                                                                 <  2:265911/150‑1 (MQ=255)
                                                                                  gCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAg                                                                <  2:853786/150‑1 (MQ=255)
                                                                                  gCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGCCCAGACGGCTGTCCAGAAACGCCGCGCGAGCCTTGCGCGACCGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAg                                                                <  1:174932/150‑1 (MQ=255)
                                                                                    cTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATgg                                                                       <  2:1287998/141‑1 (MQ=255)
                                                                                    cTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATgg                                                                       >  1:1287998/1‑141 (MQ=255)
                                                                                      tGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGTGATTTCGGCATGGCCCATACCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCt                                                            <  1:447995/150‑1 (MQ=255)
                                                                                          cGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCaa                                                        <  2:606856/150‑1 (MQ=255)
                                                                                             tGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGa                                                     <  2:1369676/150‑1 (MQ=255)
                                                                                               tACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATc                                                   <  1:868181/150‑1 (MQ=255)
                                                                                               tACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATc                                                   >  2:854080/1‑150 (MQ=255)
                                                                                                    gcgcACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGtg                                              <  2:184305/150‑1 (MQ=255)
                                                                                                     cgcACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCCAGGATCAGGtgc                                             >  2:521704/1‑150 (MQ=255)
                                                                                                      gcACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGtgct                                            <  1:50855/150‑1 (MQ=255)
                                                                                                      gcACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGtgct                                            >  2:1401141/1‑150 (MQ=255)
                                                                                                         cccGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCa                                         >  2:125830/1‑150 (MQ=255)
                                                                                                          ccGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCAt                                        <  2:149083/150‑1 (MQ=255)
                                                                                                              gATACCGAAACGCCAAGGCGCTTGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGATCCTTGCTCGACAGCTTGC‑‑CGAGCAGGTGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTcc                                    <  2:2772/150‑2 (MQ=255)
                                                                                                                 aCCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCAAGCCGGCCAGACGGCTGTGACAAAACGCCGCGCGAGACTT‑C‑CGGTGCCTTGCAGCTAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCCAGGGCTGGAGTTCTGGc                                 >  1:324574/1‑150 (MQ=255)
                                                                                                                    gAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAg                              <  2:96303/150‑1 (MQ=255)
                                                                                                                    gAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAg                              <  1:125830/150‑1 (MQ=255)
                                                                                                                         gCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGGCGTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTg                         >  1:75842/1‑150 (MQ=255)
                                                                                                                            aaGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTc                      >  2:1176434/1‑150 (MQ=255)
                                                                                                                                 gcTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTTGCGAGGGCTGTTCGGCtt                 >  1:433494/1‑150 (MQ=255)
                                                                                                                                   tCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGc                                                                                    >  1:685755/1‑81 (MQ=255)
                                                                                                                                   tCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGc                                                                                    <  2:685755/81‑1 (MQ=255)
                                                                                                                                   tCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGCTTCa               <  1:243376/150‑1 (MQ=255)
                                                                                                                                    cGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGCTTCAg              <  2:730642/150‑1 (MQ=255)
                                                                                                                                    cGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGCTTCAg              >  1:341740/1‑150 (MQ=255)
                                                                                                                                     ggCAATTTCGGCATGGCCCATGCCGTCCCGACGTCAGTCCATAAACCCCTCTCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAATGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGCTTCAgg             <  2:866582/150‑1 (MQ=255)
                                                                                                                                     ggCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAATGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGCTTCAgg             >  1:866582/1‑150 (MQ=255)
                                                                                                                                     ggCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGCTTCAgg             <  2:1256767/150‑1 (MQ=255)
                                                                                                                                     ggCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGCTTCAgg             <  1:1176434/150‑1 (MQ=255)
                                                                                                                                     ggCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGCTTCAgg             <  2:509886/150‑1 (MQ=255)
                                                                                                                                      gCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGATTCAggg            >  2:219300/1‑150 (MQ=255)
                                                                                                                                          ggacGGCATGGGCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGCTTCAGGGATCa        >  2:1406702/4‑150 (MQ=255)
                                                                                                                                           gacGGCATGGGCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCt                                    >  1:371309/3‑121 (MQ=255)
                                                                                                                                           gacGGCATGGGCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCt                                    <  2:371309/119‑1 (MQ=255)
                                                                                                                                           gacGGCATGGGCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAg                              >  1:28648/3‑127 (MQ=255)
                                                                                                                                           gacGGCATGGGCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAg                              <  2:28648/125‑1 (MQ=255)
                                                                                                                                            tCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGCTTCAGGGATCAGt      <  1:545508/150‑1 (MQ=255)
                                                                                                                                               gCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGc                                                                                           >  1:481152/1‑64 (MQ=255)
                                                                                                                                               gCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGc                                                                                           <  2:481152/64‑1 (MQ=255)
                                                                                                                                                cATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGCTTCAGGGATCAGTGCCa  <  2:701052/150‑1 (MQ=255)
                                                                                                                                                 aTGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGCTTCAg              >  1:1109149/1‑137 (MQ=255)
                                                                                                                                                 aTGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGCTTCAg              <  2:1109149/137‑1 (MQ=255)
                                                                                                                                                     cccATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTg                                   <  1:1289507/112‑1 (MQ=255)
                                                                                                                                                     cccATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTg                                   >  2:1289507/1‑112 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATCGCGGCTTCCAGGACCCGGGACGAAACCGGGAGGTTGCTCACGTCGGCTCCCCGTACAAGGCCACGTAGCACTGGCGCATGCCTTGGGCGATGTACTGGCGCACCCGCGATACCGAAACGCCAAGGCGCTCGGCAATTTCGGCATGGCCCATGCCGTCCAGACGGCTGTGCAGAAACGCCGCGCGAGCCTTGCTCGACAGCTTGC‑‑CGAGCAGGCGGTCGATGGCTTTGAGGTCTTCAAGGATCAGGTGCTGCATTTCTGGCGAGGGCTGTTCGGCTTCAGGGATCAGTGCCA  >  NC_002947/428966‑429259

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