Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 693,590 C→G 8.2% 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.2% 200.0 / 9.8 73G151A (GGC→GCC) PP_0595LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (35/32);  new base G (3/3);  total (38/35)
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.49e-01

TGCTGGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGCTCGCCGCG  >  NC_002947/693449‑693735
                                                                                                                                              |                                                                                                                                                 
tgctgGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTc                                                                                                                                            >  1:198626/1‑149 (MQ=255)
  ctgGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGGCGTGCTCGa                                                                                                                                          <  3:341300/149‑1 (MQ=255)
  ctgGTCGTCGGCGTAGAACAGCGGGGGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGTTTGATACTCAAGCCCCGATAACGGCGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTg                                                                                                                                               <  1:279376/145‑1 (MQ=255)
     gtcgtcGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGGCGTGCTCGAc                                                                                                                                         <  6:336580/147‑1 (MQ=255)
     ggcgtcgGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGAc                                                                                                                                         <  3:275931/145‑1 (MQ=255)
       cgtcgGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACctgg                                                                                                                                     >  7:224620/1‑149 (MQ=255)
                 gAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGa                                                                                                                           >  2:207795/1‑149 (MQ=255)
                  aaCAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGTCCTGGCTGGGGGc                                                                                                                             >  6:197237/1‑146 (MQ=255)
                  aaCAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACctgg                                                                                                                                     >  8:188766/1‑138 (MQ=255)
                  aaCAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACctgg                                                                                                                                     <  7:188766/138‑1 (MQ=255)
                   aCAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGAt                                                                                                                          >  1:84842/1‑148 (MQ=255)
                       cGGGTGGCCAACGGCGCATTACAGCAGCTAGCGTTCGCTGTACAAGGGTTGATACTCATGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGc                                                                                                                     <  6:51731/149‑1 (MQ=255)
                       cGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGc                                                                                                                     <  1:207795/149‑1 (MQ=255)
                        gggTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCg                                                                                                                    >  7:117245/1‑149 (MQ=255)
                             gCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCt                                                                                                               <  4:57003/149‑1 (MQ=255)
                               cAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGAcc                                                                                                                                        <  1:184565/122‑1 (MQ=255)
                               cAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAAAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTg                                                                                                              >  3:244690/1‑148 (MQ=255)
                                aaCCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGGCGTGCTCGa                                                                                                                                          <  8:200640/119‑1 (MQ=255)
                                aaCCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGGCGTGCTCGa                                                                                                                                          >  7:200640/1‑119 (MQ=255)
                                 aCCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGAt                                                                                                           <  6:181332/149‑1 (MQ=255)
                                     cgcAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGACAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATgcgc                                                                                                       <  3:98771/149‑1 (MQ=255)
                                      gcAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGca                                                                                                      <  6:73345/149‑1 (MQ=255)
                                         gTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACAt                                                                                                   <  5:299234/149‑1 (MQ=255)
                                            cagcagCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAg                                                                                                >  6:29477/1‑149 (MQ=255)
                                             agcagcGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGa                                                                                                                           >  7:249548/1‑121 (MQ=255)
                                             agcagcGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGa                                                                                                                           <  8:249548/121‑1 (MQ=255)
                                               cagcGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTCCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGc                                                                                             >  5:5331/1‑149 (MQ=255)
                                                 gcgagcgTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGccc                                                                                           <  8:224620/149‑1 (MQ=255)
                                                   gagcgTTCGCTGTACAACGGTTGATAATCAAGCCCCGATAACGGGGTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCg                                                                                         <  1:147386/149‑1 (MQ=255)
                                                   gagcgTTCGCTGTAAAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGCGCACGCCGTGCGCGACCTGGCTGGGGGAGCTCATGCGGATCTGGATGCGCCCCTAAGGGcccc                                                                                          >  4:340446/1‑148 (MQ=255)
                                                      cgTGCGCTGAACAACGGTTTATACTCACCCCCCGATACGGGGGGGGTCTGGGGTACGCCGCCAAGGTGCCGGGGGCCGTCGAGCACGCGGGGCTCCACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATc                                                                                      <  6:248869/149‑1 (MQ=255)
                                                       gTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATcc                                                                                     <  6:316585/149‑1 (MQ=255)
                                                        ttCGCTGTACAACGGTTGATACTCAAGCCCCTATACCGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCt                                                                                    <  3:340446/149‑1 (MQ=255)
                                                            cTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAg                                                                                >  2:354000/1‑149 (MQ=255)
                                                              gTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGtt                                                                              >  2:231135/1‑149 (MQ=255)
                                                                  aaCGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGAcc                                                                                                                                        <  5:43203/87‑1 (MQ=255)
                                                                      gTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATGTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAgg                                                                      >  5:159413/1‑149 (MQ=255)
                                                                        tGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGTCCTGGCTGGGGGCCATCATGCGGATCTGGATGCGCACAACAggg                                                                                              >  5:109528/1‑123 (MQ=255)
                                                                        tGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAggg                                                                                              <  6:109528/123‑1 (MQ=255)
                                                                          aTACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGCTCCTTGAGTTCCGCCAGGGCAt                                                                  >  7:10827/1‑149 (MQ=255)
                                                                             cTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGGGATCATGCGGTTCTGGTTGGGCACCTCCGGGCCCCGGTGCTTGAGTTCCGCCAGGGCATggg                                                               >  4:11652/1‑149 (MQ=255)
                                                                              tCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCCCGCCGGGGGAGACCTGGCGGGGGGGGGTCATTCGGGTGTTGGTGGGCACCTAAGGGGCCCGGTCCGGGAGTTGCGCCGCGGCAGGGGt                                                              >  8:301607/1‑149 (MQ=255)
                                                                              tCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATg                                                                 <  3:228830/146‑1 (MQ=255)
                                                                              tCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATg                                                                 >  4:228830/1‑146 (MQ=255)
                                                                              tCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGt                                                              >  5:228807/1‑149 (MQ=255)
                                                                              tCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGt                                                              >  6:296604/1‑149 (MQ=255)
                                                                                aaGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGa                                                            >  6:170208/1‑149 (MQ=255)
                                                                                    cccGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCCCGCCGTGCGCGACCGGGCTGGGGGCGATCATGCGGCTCGGGCTGCGCACATCAGGGCCCCAATCCTTGACTTCCGCCCGGGCAAGGGGGATACg                                                        >  3:208712/1‑149 (MQ=255)
                                                                                       gATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCAt                                                     >  2:144198/1‑149 (MQ=255)
                                                                                          aaCGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCAtgtg                                                  <  5:29477/149‑1 (MQ=255)
                                                                                                gCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCa                                             <  1:231135/149‑1 (MQ=255)
                                                                                                  tGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAg                                                                                >  4:104805/1‑112 (MQ=255)
                                                                                                  tGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAg                                                                                <  3:104805/112‑1 (MQ=255)
                                                                                                  tGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAgg                                           <  3:82688/149‑1 (MQ=255)
                                                                                                  tGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAgg                                           <  2:84842/149‑1 (MQ=255)
                                                                                                    gTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAggg                                          <  2:198626/148‑1 (MQ=255)
                                                                                                    gTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAgg                                           <  1:354000/147‑1 (MQ=255)
                                                                                                       tgtgGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGAt                                                           <  4:362993/128‑1 (MQ=255)
                                                                                                       tgtgGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATAATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGAt                                                           >  3:362993/1‑128 (MQ=255)
                                                                                                              cgacgCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCAGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTg                               >  4:109535/1‑149 (MQ=255)
                                                                                                                acgCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACTTCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTc                             <  7:301607/149‑1 (MQ=255)
                                                                                                                   ccAACTTGCAGGCTGCCGTCTCTCTCGCCGTTATAGTCCTGGCTGGGGGCGATCATGCGGATCTGGAGGCGCACATCAGGGCCCCTATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGGTGTCg                            <  4:208712/147‑1 (MQ=255)
                                                                                                                     aaCGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCa                        <  8:117245/149‑1 (MQ=255)
                                                                                                                        gTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGc                     >  6:163039/1‑149 (MQ=255)
                                                                                                                            aGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGAt                 <  6:228807/149‑1 (MQ=255)
                                                                                                                                   cGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGc          >  7:99009/1‑149 (MQ=255)
                                                                                                                                   cGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGc          >  4:173777/1‑149 (MQ=255)
                                                                                                                                    gTCGAGCACGGCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGCt         >  8:287276/1‑149 (MQ=255)
                                                                                                                                    gTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGCt         >  7:329300/1‑149 (MQ=255)
                                                                                                                                    gTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCGTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGCt         >  2:237266/1‑149 (MQ=255)
                                                                                                                                    gTCGAGCACGCCGTGCTCGACCTGGATGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCCGCt         >  8:234953/1‑149 (MQ=255)
                                                                                                                                     tCGAGCACGGCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGCTGTTCAGCTc        >  6:17100/1‑149 (MQ=255)
                                                                                                                                            cGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGCTCGCcgcg  >  1:82737/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
TGCTGGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCATGCGGATCTGGATGCGCACATCAGGGCCCCGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGCTCGCCGCG  >  NC_002947/693449‑693735

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