Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,248,921 G→T 7.5% intergenic (‑184/‑183) PP_1089 ← / → bvgA hypothetical protein/virulence factors transcription regulator BvgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,248,9210GT 7.5% 187.4 / 2.2 60intergenic (‑184/‑183)PP_1089/bvgAhypothetical protein/virulence factors transcription regulator BvgA
Reads supporting (aligned to +/- strand):  ref base G (32/23);  new base T (3/2);  total (35/25)
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.00e+00

TTACTGCTGTACTCAACTAAAGGCTGGATTGTTTCCCGGGGGAGGAATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCTTTCGCTGGTAAATCGCCGCGCCAGCAAGCGTCCTGTCGGTTATCG  >  NC_002947/1248773‑1249067
                                                                                                                                                    |                                                                                                                                                  
ttACTGCTGTACTCAACTAAAGGCTGGATTGTTTCCCGGGGGAGGAATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAgtgg                                                                                                                                                    <  1:242633/149‑1 (MQ=255)
              aaCTAAAGTCTGGTTTGTTTACCGGGGGTGGAATGATCGTCAGAGTACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGTGTGGTACtttt                                                                                                                                         <  7:20076/146‑1 (MQ=255)
              aaCTAAAGGCTGGATTGTTTCCCGGGGGAGGAATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACtttttt                                                                                                                                       >  5:243976/1‑148 (MQ=255)
              aaCTAAAGGCTGGATTGTTTCCCGGGGGAGGAATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACtttttt                                                                                                                                       >  6:212431/1‑148 (MQ=255)
              aaCTAAAGGCTGGATTGTTTCCCGGGGGAGGAATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACtttt                                                                                                                                         >  1:231350/1‑146 (MQ=255)
              aaCTAAAGGCTGGATTGTTTCCCGGGGGAGGAATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACtttt                                                                                                                                         <  2:231350/146‑1 (MQ=255)
              aaCTAAAGGCTGGATTGTTTACCGGGTGAGGAATGATCGTCAGAGAAACTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAATATCAGTATTAAACTGGTACGTTTATGTCAGTTTGTGTTACtttt                                                                                                                                         >  8:20076/1‑146 (MQ=255)
                 tAAAGGCTGGATTGTTTCCCGGGGGAGGAATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCgg                                                                                                                                   >  1:221283/1‑149 (MQ=255)
                      gCTGGATTGTTTCCCGGGGGAGGAATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGtaat                                                                                                                               >  1:300758/1‑148 (MQ=255)
                          gATTGTTTCCCGGGGGAGGAATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGt                                                                                                                           >  7:21048/1‑148 (MQ=255)
                                 tCCCGGGGGAGGAATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAatttatt                                                                                                                   >  5:93709/1‑149 (MQ=255)
                                   ccGGGGGAGGAATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTg                                                                                                                 >  5:5946/1‑149 (MQ=255)
                                     gggggAGGAATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGtt                                                                                                               <  2:221283/149‑1 (MQ=255)
                                          aggaATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAAc                                                                                                          <  5:212431/149‑1 (MQ=255)
                                          aggaATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAAc                                                                                                          <  6:93709/149‑1 (MQ=255)
                                           cttccgATCTTCAGAGAACCTCGTTGTGTTGCTTTTTCGACTTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACaa                                                                                                           <  3:173223/142‑1 (MQ=255)
                                            gaATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCa                                                                                                        >  5:83471/1‑149 (MQ=255)
                                              aTGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTATGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCg                                                                                                                                    <  2:115380/119‑1 (MQ=255)
                                              aTGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTATGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCg                                                                                                                                    >  1:115380/1‑119 (MQ=255)
                                               tGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAg                                                                                                     >  3:258209/1‑149 (MQ=255)
                                                gATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGt                                                                                                    >  7:176123/1‑149 (MQ=255)
                                                    gTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAATATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCAtt                                                                                                <  8:21048/149‑1 (MQ=255)
                                                      cAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACaa                                                                                                           >  4:173223/1‑136 (MQ=255)
                                                         agaACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATtgt                                                                                              <  5:193937/146‑1 (MQ=255)
                                                         agaACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATtgt                                                                                              >  6:193937/1‑146 (MQ=255)
                                                              cTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGtt                                                                                                               <  8:234810/124‑1 (MQ=255)
                                                              cTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGtt                                                                                                               >  7:234810/1‑124 (MQ=255)
                                                               tCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGaa                                                                                     >  1:217345/1‑149 (MQ=255)
                                                                    gtgtTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGttt                                                                                 <  6:5946/148‑1 (MQ=255)
                                                                       ttGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGTGTGGTACCTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTg                                                                              >  1:174882/1‑148 (MQ=255)
                                                                       ttGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTg                                                                              >  8:247384/1‑148 (MQ=255)
                                                                       ttGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTg                                                                              <  4:258209/148‑1 (MQ=255)
                                                                         gCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGt                                                                             >  2:140522/1‑147 (MQ=255)
                                                                          cTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTaa                                                                           <  1:140522/148‑1 (MQ=255)
                                                                           tttttCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCATTATTCAACTGGTACGTTTATGTAAGTGTGTGGTACTTTTTTTCGTTAATAAGTAATTTATTTTTTACAAACATAGTCATTGTGTGTAGGAAGTTTCTGTAAGt                                                                         >  7:14157/1‑149 (MQ=255)
                                                                           tttttCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGt                                                                         >  7:49389/1‑149 (MQ=255)
                                                                                 gACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTg                                                                     >  4:147063/1‑147 (MQ=255)
                                                                                 gACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTa                                                                   <  8:176123/149‑1 (MQ=255)
                                                                                          gAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCAtt                                                                                                <  6:39769/111‑1 (MQ=255)
                                                                                            gTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCt                                                         <  2:217345/148‑1 (MQ=255)
                                                                                               ggAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTc                                                      >  5:54044/1‑148 (MQ=255)
                                                                                                   gTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCttttt                                                           >  1:185056/1‑139 (MQ=255)
                                                                                                   gTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCttttt                                                           <  2:185056/139‑1 (MQ=255)
                                                                                                     tttCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCt                                               >  3:33495/1‑149 (MQ=255)
                                                                                                      ttCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCtt                                              >  7:52650/1‑149 (MQ=255)
                                                                                                          tttGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCTTTCGc                                          >  3:107947/1‑149 (MQ=255)
                                                                                                           ttGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCt                                                                               <  5:10966/111‑1 (MQ=255)
                                                                                                           ttGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCt                                                                               >  6:10966/1‑111 (MQ=255)
                                                                                                             gCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATg                                                  <  2:208814/138‑1 (MQ=255)
                                                                                                             gCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATg                                                  >  1:208814/1‑138 (MQ=255)
                                                                                                             gCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCTTTCGCTg                                        <  6:243976/148‑1 (MQ=255)
                                                                                                             gCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCTTTCGCTg                                        <  4:107947/148‑1 (MQ=255)
                                                                                                                   gATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCTTTCGCTGGTAAATc                                 <  3:185879/149‑1 (MQ=255)
                                                                                                                   gATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCTTTCGCTGGTAAATc                                 >  6:161367/1‑149 (MQ=255)
                                                                                                                        gTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCTTTCGCTGGTAAATcgccgc                            <  6:54044/149‑1 (MQ=255)
                                                                                                                         tATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCTTTCGCTGGTAAATCGCcgcg                           >  4:2074/1‑149 (MQ=255)
                                                                                                                           ttCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCTTTCGCTGGTAAATCGCCgcgc                          >  5:198551/1‑148 (MQ=255)
                                                                                                                           ttCAACTGGTACGTTTATGTCAGTGGGCGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCTTTCGCTGGTAAATCGCCgcgc                          >  6:456/1‑148 (MQ=255)
                                                                                                                              aaCTGGTACGTTTATGTCAGTGTGTGGTACCTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCTTTCGCTGGTAAATCGCCGCGCCAg                       <  2:174882/148‑1 (MQ=255)
                                                                                                                                 tGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCGATGGCTTTCGCTGGTAAATCGCCGCGCGAgagag                   >  4:234666/1‑145 (MQ=255)
                                                                                                                                            tGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCTTTCGCTGGTAAATCGCCGCGCCAGCAAGCGTCCTGTCgg        <  3:234666/149‑1 (MQ=255)
                                                                                                                                                  tgggtggTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCTTTCGCTGGTAAATCGCCGCGCCAGCAAGCGTCCTGTCGGTTATCg  >  2:169017/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTACTGCTGTACTCAACTAAAGGCTGGATTGTTTCCCGGGGGAGGAATGATCGTCAGAGAACCTCGTTGTGTTGCTTTTTCGACGTGTTTGAGTCGGAAGTTTTCCTTTGCTGAAGATCAGTATTCAACTGGTACGTTTATGTCAGTGGGTGGTACTTTTTTTCGGTAATAAGTAATTTATTTGTTACAACCATAGTCATTGTGTGTAGGAAGTTTCTGTAAGTTCTGTAGGCTTTTTCTGTCCATGGCTTTCGCTGGTAAATCGCCGCGCCAGCAAGCGTCCTGTCGGTTATCG  >  NC_002947/1248773‑1249067

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: ‑