Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 693,590 C→G 8.7% G151A (GGC→GCC)  PP_0595 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947693,5900CG 8.7% 207.5 / 9.8 69G151A (GGC→GCC) PP_0595LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base G (3/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

GTCTATCTGCTGGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGCTCGCCGCGC  >  NC_002947/693442‑693736
                                                                                                                                                    |                                                                                                                                                  
gtcTATCTGCTGGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGc                                                                                                                                                    <  2:59785/149‑1 (MQ=255)
 tcTATCTGCTGGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGcc                                                                                                                                                   >  7:163277/1‑149 (MQ=255)
  cTATCTGCTGGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGGCg                                                                                                                                                  >  8:291578/1‑149 (MQ=255)
     tCTGCTGGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTg                                                                                                                                                <  8:235048/148‑1 (MQ=255)
         ctgGTCGTCGGGGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGGCGTGCTCGa                                                                                                                                           <  4:96286/149‑1 (MQ=255)
         ctgGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGTTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGc                                                                                                                                               <  1:81348/145‑1 (MQ=255)
         ctgGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGGCGTGCTCGa                                                                                                                                           <  1:357093/149‑1 (MQ=255)
         ctgGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGc                                                                                                                                               <  1:12754/145‑1 (MQ=255)
            gtcgtcGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGGCGTGCTCGAc                                                                                                                                          <  7:291578/147‑1 (MQ=255)
                tcgGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCTTCGAGCACGCCGTGCTCGACctggct                                                                                                                                    <  2:305294/149‑1 (MQ=255)
                 cgGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACct                                                                                                                                        >  4:13925/1‑144 (MQ=255)
                 cgGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACct                                                                                                                                        <  3:13925/144‑1 (MQ=255)
                         aaCAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGa                                                                                                                            <  6:21642/148‑1 (MQ=255)
                         aaCAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCAGGCTGGGGGCGa                                                                                                                            >  5:21642/1‑148 (MQ=255)
                          aCAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGAt                                                                                                                           >  5:153543/1‑148 (MQ=255)
                           cAGCGGGGGGCAAACCGCGCAGCACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATACCGGGCTGGTCTGCGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCa                                                                                                                         <  7:53732/149‑1 (MQ=255)
                            aGCGGGTGGCCAACCGCGCAGTACAGAAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGGTCGACCGGTCTGGGGGCGATCAt                                                                                                                        >  3:237263/1‑149 (MQ=255)
                                 gTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGa                                                                                                                   >  8:50970/1‑149 (MQ=255)
                                   ggCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATc                                                                                                                 >  2:132103/1‑149 (MQ=255)
                                    gCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCt                                                                                                                <  2:311300/149‑1 (MQ=255)
                                     ccAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCCAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTg                                                                                                               >  3:44209/1‑149 (MQ=255)
                                      cAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTgg                                                                                                              <  6:346419/149‑1 (MQ=255)
                                      cAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTgg                                                                                                              >  5:11749/1‑149 (MQ=255)
                                      cAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTg                                                                                                               >  2:294769/1‑148 (MQ=255)
                                            cgcAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATgcgc                                                                                                        <  5:125670/149‑1 (MQ=255)
                                                  aCAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACAt                                                                                                    >  4:40539/1‑147 (MQ=255)
                                                   cagcagCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCAGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAg                                                                                                 >  5:23890/1‑149 (MQ=255)
                                                    agcagcGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAgg                                                                                                <  6:153543/149‑1 (MQ=255)
                                                    agcagcGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAgg                                                                                                <  7:50970/149‑1 (MQ=255)
                                                        gcgagcgTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGccc                                                                                            <  3:40539/149‑1 (MQ=255)
                                                         cgagcgTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGcccc                                                                                           <  6:11749/149‑1 (MQ=255)
                                                                 cGCTGTACAACGGTTGATACTCAAGCACCGATAACGGGCTTGTATGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGTCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATACtt                                                                                    <  4:237263/148‑1 (MQ=255)
                                                                      tACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTc                                                                              >  7:140860/1‑149 (MQ=255)
                                                                         aaCGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCGGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGc                                                                           <  4:44209/149‑1 (MQ=255)
                                                                            ggTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTCCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAg                                                                        <  8:163277/149‑1 (MQ=255)
                                                                              ttGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAggg                                                                      >  8:152509/1‑149 (MQ=255)
                                                                                 aTACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCa                                                                    <  1:132103/148‑1 (MQ=255)
                                                                                  tACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGc                                                                           <  1:62113/140‑1 (MQ=255)
                                                                                  tACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGc                                                                           >  2:62113/1‑140 (MQ=255)
                                                                                    cTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATgg                                                                 <  8:140860/148‑1 (MQ=255)
                                                                                     tCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGt                                                               >  5:332013/1‑149 (MQ=255)
                                                                                      cAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCGTGGGt                                                               >  3:117537/1‑148 (MQ=255)
                                                                                       aaGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTg                                                              >  7:74517/1‑148 (MQ=255)
                                                                                              gATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCAt                                                      <  7:152509/149‑1 (MQ=255)
                                                                                                  aCGGGCTGGTCTGTGGTACGACGCCAACGTCCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGc                                                  >  1:31296/1‑149 (MQ=255)
                                                                                                        tGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAgg                                            <  6:23890/149‑1 (MQ=255)
                                                                                                                       cgCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAgg                                            >  8:270156/1‑134 (MQ=255)
                                                                                                                       cgCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAgg                                            <  7:270156/134‑1 (MQ=255)
                                                                                                                        gCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCg                             <  8:74517/148‑1 (MQ=255)
                                                                                                                           aaCGTGCAGGCTGCCGTCGCGCACGCCGTGCTCGCCCTGGCTGGGGGCGCTCATGCGGGTCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCggg                           >  3:183337/1‑146 (MQ=255)
                                                                                                                           aaCGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAggg                                                                      <  3:220194/104‑1 (MQ=255)
                                                                                                                           aaCGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAggg                                                                      >  4:220194/1‑104 (MQ=255)
                                                                                                                                 cAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCGTGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGa                   <  4:117537/149‑1 (MQ=255)
                                                                                                                                   ggCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTcc                                                                             <  5:15837/89‑1 (MQ=255)
                                                                                                                                      ccaggTCGAGCACGGCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGtt               >  8:314673/5‑148 (MQ=255)
                                                                                                                                          gTCGAGCACGGCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGCt          >  5:380711/1‑149 (MQ=255)
                                                                                                                                           tCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGCTc         >  5:113827/1‑149 (MQ=255)
                                                                                                                                           tCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTAAGGCCGATGTTCAGCTc         >  7:230361/1‑149 (MQ=255)
                                                                                                                                             gAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAgg                       <  2:261063/133‑1 (MQ=255)
                                                                                                                                             gAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAgg                       <  2:255942/133‑1 (MQ=255)
                                                                                                                                             gAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAgg                       >  1:255942/1‑133 (MQ=255)
                                                                                                                                             gAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAgg                       >  1:261063/1‑133 (MQ=255)
                                                                                                                                             gAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGCTcgc       >  5:187852/1‑149 (MQ=255)
                                                                                                                                             gAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGCTcg        <  6:113827/148‑1 (MQ=255)
                                                                                                                                               gCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGc                                                        >  5:190853/1‑98 (MQ=255)
                                                                                                                                               gCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGc                                                        <  6:190853/98‑1 (MQ=255)
                                                                                                                                                 aCGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGCTCGCcgcg   <  6:332013/149‑1 (MQ=255)
                                                                                                                                                   gCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCTGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGCTCGCCgcgc  >  8:160093/1‑148 (MQ=255)
                                                                                                                                                   gCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATAGGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGCTCGCCgcgc  >  6:28225/1‑148 (MQ=255)
                                                                                                                                                   gCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCAcc                                      >  2:244523/1‑112 (MQ=255)
                                                                                                                                                   gCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCAcc                                      <  1:244523/112‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTCTATCTGCTGGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGCTCGCCGCGC  >  NC_002947/693442‑693736

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑