Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,677,630 G→C 10.2% intergenic (+23/‑108) thrC → / → PP_5471 threonine synthase/two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,677,6300GC10.2% 198.0 / 12.3 69intergenic (+23/‑108)thrC/PP_5471threonine synthase/two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base C (3/4);  total (31/38)
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.09e-01

AGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCACCAGGCTCGCCCATC  >  NC_002947/1677486‑1677778
                                                                                                                                                |                                                                                                                                                    
aGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTgcagcg                                                                                                                                                    >  6:209291/1‑147 (MQ=255)
  cTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTACATGGCGTTTGGGGCTgcagcgcagc                                                                                                                                                <  6:121818/149‑1 (MQ=255)
     cccGCGCACCTCAGCGACGTTTTCAGCCGTTAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGcccc                                                                                                                                             <  1:356648/149‑1 (MQ=255)
     cccGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCGCTGCAGcccc                                                                                                                                             <  3:43266/149‑1 (MQ=255)
     cccGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCGCTGCAGccc                                                                                                                                              <  7:316526/148‑1 (MQ=255)
           cACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGtc                                                                                                                                       >  3:163338/1‑149 (MQ=255)
              cTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGtct                                                                                                                                      >  3:248668/1‑147 (MQ=255)
              cTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGtct                                                                                                                                      <  4:248668/147‑1 (MQ=255)
              cTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCa                                                                                                                                    <  6:16676/149‑1 (MQ=255)
                cAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTAGGGGCTGCAGCGCAGCCCCAACGTCTCAg                                                                                                                                   >  7:116256/1‑148 (MQ=255)
                 aGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGcc                                                                                                                                 >  6:180380/1‑149 (MQ=255)
                 aGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTATAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGcc                                                                                                                                 >  6:369219/1‑149 (MQ=255)
                   cGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCtt                                                                                                                               <  7:149017/149‑1 (MQ=255)
                    gACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCttt                                                                                                                              <  7:172347/149‑1 (MQ=255)
                     aCCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGCCCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCtttt                                                                                                                             <  8:15242/149‑1 (MQ=255)
                     aCCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCtttt                                                                                                                             <  2:295368/149‑1 (MQ=255)
                      ccTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTc                                                                                                                            <  7:130837/149‑1 (MQ=255)
                          ttCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTccc                                                                                                                          <  1:307725/147‑1 (MQ=255)
                          ttCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCcctc                                                                                                                        >  1:198112/1‑149 (MQ=255)
                          ttCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCcct                                                                                                                         <  5:180380/148‑1 (MQ=255)
                              gCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGt                                                                                                                    <  5:156217/149‑1 (MQ=255)
                              gCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTATAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGt                                                                                                                    <  5:369219/149‑1 (MQ=255)
                                 gTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTa                                                                                                                 >  6:59278/1‑149 (MQ=255)
                                  tGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTa                                                                                                                 >  6:62116/1‑148 (MQ=255)
                                   gAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTAt                                                                                                                >  1:135557/1‑148 (MQ=255)
                                      gagCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCcct                                                                                                                         >  5:11195/1‑136 (MQ=255)
                                      gagCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCcct                                                                                                                         <  6:11195/136‑1 (MQ=255)
                                                        gCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCt                                                                                           <  5:59278/148‑1 (MQ=255)
                                                          cAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGcc                                                                                        >  3:179488/1‑149 (MQ=255)
                                                          cAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGcc                                                                                        >  6:151978/1‑149 (MQ=255)
                                                          cAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGcc                                                                                        >  5:276870/1‑149 (MQ=255)
                                                                ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGgcgc                                                                                  >  8:118605/1‑149 (MQ=255)
                                                                ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGgcgc                                                                                  >  8:123650/1‑149 (MQ=255)
                                                                 cTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGgcgc                                                                                  >  8:210750/1‑148 (MQ=255)
                                                                   gAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGcc                                                                                 >  4:66035/1‑147 (MQ=255)
                                                                    aaGGCTGTACAGGGCCTTGTCAGCCAGCAAGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCaa                                                                               <  6:296832/148‑1 (MQ=255)
                                                                     aGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGc                                                                             <  5:238822/149‑1 (MQ=255)
                                                                       gCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCAt                                                                           <  6:276870/149‑1 (MQ=255)
                                                                           tACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATaa                                                                         <  5:62116/147‑1 (MQ=255)
                                                                               gggCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGc                                                                   <  7:221308/149‑1 (MQ=255)
                                                                               gggCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGc                                                                   <  6:129790/149‑1 (MQ=255)
                                                                               gggCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGc                                                                   <  2:198112/149‑1 (MQ=255)
                                                                                ggCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGcc                                                                  >  1:12481/1‑149 (MQ=255)
                                                                                 gCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCg                                                                 <  5:209291/149‑1 (MQ=255)
                                                                                       tCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTa                                                                     >  1:225623/1‑139 (MQ=255)
                                                                                       tCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTa                                                                     <  2:225623/139‑1 (MQ=255)
                                                                                         agccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCGCTGCAGcccc                                                                                                                                             >  5:343399/1‑65 (MQ=38)
                                                                                         agccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCGCTGCAGcccc                                                                                                                                             <  6:343399/65‑1 (MQ=38)
                                                                                           ccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCAt                                                       >  8:73427/1‑149 (MQ=255)
                                                                                           caagcACGGCAACAGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCACCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCAt                                                       >  4:301296/3‑149 (MQ=255)
                                                                                               cACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGcc                                                                                        <  4:209291/111‑1 (MQ=255)
                                                                                               cACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGcc                                                                                        >  3:209291/1‑111 (MQ=255)
                                                                                                  ggCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCa                                                     <  3:168750/144‑1 (MQ=255)
                                                                                                  ggCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCa                                                     >  4:168750/1‑144 (MQ=255)
                                                                                                     aaCCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCATAAACATCCTATTCAGCCTGCCGTGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGc                                             <  7:210750/149‑1 (MQ=255)
                                                                                                         gcgGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCAt                                         <  4:179488/149‑1 (MQ=255)
                                                                                                          cgGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATg                                        >  5:326962/1‑149 (MQ=255)
                                                                                                                    cTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCt                              <  6:326962/149‑1 (MQ=255)
                                                                                                                    cTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCt                              <  2:12481/149‑1 (MQ=255)
                                                                                                                    cTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCt                              <  7:123650/149‑1 (MQ=255)
                                                                                                                    cTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCt                              <  7:118605/149‑1 (MQ=255)
                                                                                                                        aaTTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCGATTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTAtgt                          >  1:259873/1‑149 (MQ=255)
                                                                                                                            gCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGtt                       <  7:73427/148‑1 (MQ=255)
                                                                                                                             cATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGa                     >  1:346806/1‑149 (MQ=255)
                                                                                                                                ggCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGc                                                                   >  6:333235/1‑100 (MQ=255)
                                                                                                                                ggCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGc                                                                   <  5:333235/100‑1 (MQ=255)
                                                                                                                                   gTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTAt                            >  2:310779/1‑136 (MQ=255)
                                                                                                                                   gTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTAt                            <  1:310779/136‑1 (MQ=255)
                                                                                                                                     ttGGGGCTGCGCTGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGcc                                >  2:61738/1‑130 (MQ=255)
                                                                                                                                     ttGGGGCTGCGCTGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGcc                                <  1:61738/130‑1 (MQ=255)
                                                                                                                                       ggggCTGCGCTGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCACCAGGc           >  6:6734/1‑149 (MQ=255)
                                                                                                                                                gcgcagcCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCACCAGGCTCGCCCATc  >  5:332277/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCACCAGGCTCGCCCATC  >  NC_002947/1677486‑1677778

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