Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,525,773 T→G 9.2% intergenic (+63/‑60) PP_4858 → / → PP_4859 hypothetical protein/TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,525,7730TG 9.2% 191.5 / 8.9 64intergenic (+63/‑60)PP_4858/PP_4859hypothetical protein/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base G (2/4);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

TCCGATCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGTGGCCGTGGTTTCG  >  NC_002947/5525630‑5525920
                                                                                                                                                |                                                                                                                                                   
tCCGATCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCg                                                                                                                                                <  7:195744/149‑1 (MQ=255)
   gATCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCAATGGACCGGGc                                                                                                                                             <  4:49638/149‑1 (MQ=255)
       gTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCAATGGACCGGGcctg                                                                                                                                          <  2:196571/148‑4 (MQ=255)
        tCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGt                                                                                                                                        >  6:19618/1‑149 (MQ=255)
            aTCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGtt                                                                                                                                       <  1:37439/147‑1 (MQ=255)
            aTCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGttttt                                                                                                                                    >  5:45603/1‑149 (MQ=255)
               gAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTa                                                                                                                                 >  7:60227/1‑149 (MQ=255)
               gAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTa                                                                                                                                 >  1:312638/1‑149 (MQ=255)
                 aCAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATc                                                                                                                               >  5:58763/1‑149 (MQ=255)
                 aCAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCAATGGACCGGGc                                                                                                                                             >  3:295934/1‑135 (MQ=255)
                 aCAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCAATGGACCGGGc                                                                                                                                             <  4:295934/135‑1 (MQ=255)
                     gCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTaga                                                                                                                           <  6:58763/149‑1 (MQ=255)
                     gCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAAGCCCGGTCAATGGACCGGGcct                                                                                                                                           <  6:86110/133‑3 (MQ=255)
                     gCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAAGCCCGGTCAATGGACCGGGcct                                                                                                                                           >  5:86110/1‑131 (MQ=255)
                           cGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGccc                                                                                                                     <  6:45603/149‑1 (MQ=255)
                               tCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGc                                                                                                                  <  2:312638/148‑1 (MQ=255)
                                   tGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGAc                                                                                                             >  4:121367/1‑149 (MQ=255)
                                    gAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCCCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGAcc                                                                                                            >  2:171507/1‑149 (MQ=255)
                                     aaGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCg                                                                                                           >  3:207393/1‑149 (MQ=255)
                                      aGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGt                                                                                                          <  4:207393/149‑1 (MQ=255)
                                                gATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGttttt                                                                                                                                    >  6:215819/1‑113 (MQ=255)
                                                gATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGttttt                                                                                                                                    <  5:215819/113‑1 (MQ=255)
                                                    gCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGca                                                                                            >  2:1267/1‑149 (MQ=255)
                                                    gCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGc                                                                                             <  8:60227/148‑1 (MQ=255)
                                                    gCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGc                                                                                             <  5:19618/148‑1 (MQ=255)
                                                     ccTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTa                                                                                                                             <  5:48601/115‑1 (MQ=255)
                                                     ccTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTa                                                                                                                             >  6:48601/1‑115 (MQ=255)
                                                              aCTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCCCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGcc                                                                                  <  1:171507/149‑1 (MQ=255)
                                                                tCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCcgc                                                                                >  7:76835/1‑149 (MQ=255)
                                                                tCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCAAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCcgc                                                                                >  5:287783/1‑149 (MQ=255)
                                                                    ccGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGcatcat                                                                                        >  1:18897/1‑137 (MQ=255)
                                                                    ccGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGcatcat                                                                                        <  2:18897/137‑1 (MQ=255)
                                                                    ccGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCgcgcgc                                                                            <  4:235663/149‑1 (MQ=255)
                                                                              tGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAg                                                                  >  3:43377/1‑149 (MQ=255)
                                                                                         gtccaaaaGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGccc                                                                                               <  1:289722/105‑1 (MQ=255)
                                                                                         gCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGccc                                                                                               >  2:289722/1‑109 (MQ=255)
                                                                                          ccGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCg                                                                          >  3:227354/1‑129 (MQ=255)
                                                                                          ccGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCg                                                                          <  4:227354/129‑1 (MQ=255)
                                                                                                  aacaacTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCt                                              >  6:327873/1‑149 (MQ=255)
                                                                                                    caacTTCTGTCACAACGCTTTCTCCCCTTCATGGCCGGTTCATTGACCGGGGTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGc                                                 <  8:124248/144‑1 (MQ=255)
                                                                                                    caacTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGc                                            >  6:271768/1‑149 (MQ=255)
                                                                                                    caacTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAAAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGATGCCGCGCGCCGACCAGAAGCTACAGACCCGCCAGcccag                                              >  7:124248/1‑145 (MQ=255)
                                                                                                         ccgagctcaACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTGTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGa                                        <  3:253670/141‑1 (MQ=255)
                                                                                                           tGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgc                                     >  7:8627/1‑149 (MQ=255)
                                                                                                            gTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTa                                                               >  2:13400/1‑122 (MQ=255)
                                                                                                             tCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGATTTCCCAGCCTGACCGTTACAGAACCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCAGCACGCTgccg                                   >  8:117942/1‑149 (MQ=255)
                                                                                                             tCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgccg                                   >  5:184144/1‑149 (MQ=255)
                                                                                                             tCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgc                                     <  4:43377/147‑1 (MQ=255)
                                                                                                               acaACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgccgc                                  <  8:8627/148‑1 (MQ=255)
                                                                                                               acaACGCTTGCTCCCCTCCAGGCCAGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgccgcc                                 >  4:319043/1‑149 (MQ=255)
                                                                                                                 aaCGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCgcctg                               >  2:7835/1‑149 (MQ=255)
                                                                                                                 aaCGCTTGCTCCCCTCCAGGCCCGGACCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGAATCATCCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGcccggagcga                                        >  4:253670/1‑133 (MQ=255)
                                                                                                                    gCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCgcctgcct                            >  8:300118/1‑149 (MQ=255)
                                                                                                                             cccTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgccgc                                  <  8:6180/135‑1 (MQ=255)
                                                                                                                             cccTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATgccgc                                  >  7:6180/1‑135 (MQ=255)
                                                                                                                             cccTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATg                     <  5:271768/148‑1 (MQ=255)
                                                                                                                              ccTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGga                   <  3:121367/149‑1 (MQ=255)
                                                                                                                              ccTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGga                   <  8:76835/149‑1 (MQ=255)
                                                                                                                              ccTCCAGGCCAGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGga                   <  3:319043/149‑1 (MQ=255)
                                                                                                                                 ccAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGagag                >  8:182168/1‑149 (MQ=255)
                                                                                                                                  cAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGgaga                 <  7:117942/147‑1 (MQ=255)
                                                                                                                                    ggCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGTgg             <  5:327873/149‑1 (MQ=255)
                                                                                                                                           tCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACACGCCAGGCCCTGCTCGATGCCgcctgccag                           >  3:253620/1‑126 (MQ=255)
                                                                                                                                             cATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGTGGCCGTGGttt    <  1:7835/149‑1 (MQ=255)
                                                                                                                                               ttGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGTGGCCGTGGTTTCg  >  5:112442/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TCCGATCGTCTGATCGAACAAGCTGGCCGTGTCGCTGAAGGTGGCTCGGATCGCCTGGTTGAACTCAGCCGTGTGAGCTGATCGCCATGGCCGAAAAGAACAACTTCTGTCACAACGCTTGCTCCCCTCCAGGCCCGGTCCATTGACCGGGCTTTGTTTTTTTATCTAGAGTGCCCAGCCTGACCGTTACAGAAGCCCGCATCATGCTGCCGCGCGCCGAACAGAAGCTACAGACCCGCCAGGCCCTGCTCGATGCCGCCTGCCTGCTCATGGAGAGTGGCCGTGGTTTCG  >  NC_002947/5525630‑5525920

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: