Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,249,753 A→T 6.9% intergenic (+26/+45) bvgA → / ← PP_1091 virulence factors transcription regulator BvgA/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,249,7530AT 6.9% 174.5 / 3.3 58intergenic (+26/+45)bvgA/PP_1091virulence factors transcription regulator BvgA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/28);  new base T (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTT  >  NC_002947/1249606‑1249891
                                                                                                                                                   |                                                                                                                                          
ccGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTccagccag                                                                                                                                           >  4:315481/1‑149 (MQ=255)
    cAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCaggagg                                                                                                                                       <  5:295224/149‑1 (MQ=255)
       cATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCtt                                                                                                                                    <  2:6566/149‑1 (MQ=255)
              cTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCttt                                                                                                                                   <  4:30287/143‑1 (MQ=255)
              cTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCttt                                                                                                                                   >  3:30287/1‑143 (MQ=255)
              cTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAgg                                                                                                                             <  6:3544/149‑1 (MQ=255)
               tCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAggg                                                                                                                            >  3:88477/1‑149 (MQ=255)
               tCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAggg                                                                                                                            >  3:15381/1‑149 (MQ=255)
                 aGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCTTTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGc                                                                                                                           >  8:12309/1‑148 (MQ=255)
                 aGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTTATCGCCAAGCCGTGAAGCCTCCAGCCAGGAGGCTTTGCCAgggat                                                                                                                          <  4:154785/149‑3 (MQ=255)
                 aGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAggg                                                                                                                            <  3:285559/147‑1 (MQ=255)
                   tAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTAt                                                                                                                        >  8:224721/1‑149 (MQ=255)
                    aaCAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATc                                                                                                                       >  1:149576/1‑149 (MQ=255)
                     aCAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCg                                                                                                                      >  2:212740/1‑149 (MQ=255)
                        aGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTc                                                                                                                   <  7:224721/149‑1 (MQ=255)
                             gTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGa                                                                                                              <  5:272832/149‑1 (MQ=255)
                             gTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGa                                                                                                              <  4:88477/149‑1 (MQ=255)
                              tCAGTACCTTCAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGAc                                                                                                             >  7:78432/1‑149 (MQ=255)
                                         aaGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCTGGAGGCttt                                                                                                                                   >  3:195004/1‑116 (MQ=255)
                                         aaGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCTGGAGGCttt                                                                                                                                   <  4:195004/116‑1 (MQ=255)
                                         aaGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCg                                                                                                  >  8:188879/1‑149 (MQ=255)
                                         aaGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCg                                                                                                  >  3:184154/1‑149 (MQ=255)
                                                ggCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACa                                                                                           <  2:149576/149‑1 (MQ=255)
                                                 gctgctgctCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCTGGCTGGAGGCttcaa                                                                                                                                 <  4:226108/110‑4 (MQ=255)
                                                 gctgctgctCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCTGGCTGGAGGCttcaa                                                                                                                                 >  3:226108/1‑107 (MQ=255)
                                                  ctgctgctCACCCTCAATCCCCGGTCATAGGTTGACCTGATCGAGCGGCCGCAACGGCAAGGTCAGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGTGACCCGCACTACAgg                                                                                         <  4:304545/149‑1 (MQ=255)
                                                      tgctCAAGCTCAATGCACGGTCATTGGTTGACCTGAACGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCg                                                                                      >  5:45768/1‑148 (MQ=255)
                                                          caagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGATCCATGACTACGGGACCCGCACTACAGGTCGTACTc                                                                                 <  1:31332/149‑1 (MQ=255)
                                                             gctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGaa                                                                              <  8:78432/149‑1 (MQ=255)
                                                                 aaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATc                                                                           <  3:315481/148‑1 (MQ=255)
                                                                 aaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCg                                                                          <  4:184154/149‑1 (MQ=255)
                                                                 aaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGCAATCGAAATCg                                                                          >  7:66973/1‑149 (MQ=255)
                                                                     cACGGTCTTTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGcc                                                                       <  7:12309/148‑1 (MQ=255)
                                                                            aTTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCtt                                                               >  8:187578/1‑149 (MQ=255)
                                                                            aTTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCt                                                                <  5:42177/148‑1 (MQ=255)
                                                                            aTTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCt                                                                <  5:47051/148‑1 (MQ=255)
                                                                                     ccTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGgcg                                                      >  6:137575/1‑149 (MQ=255)
                                                                                      cTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGgcgc                                                     >  2:240088/1‑149 (MQ=255)
                                                                                         aTCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGcc                                                    <  7:188879/147‑1 (MQ=255)
                                                                                              gCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttcttc                                             >  6:179139/1‑149 (MQ=255)
                                                                                              gCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttcttc                                             >  7:119272/1‑149 (MQ=255)
                                                                                              gCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGCAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttcttc                                             <  8:66973/149‑1 (MQ=255)
                                                                                                  gcgcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCa                                           <  7:187578/147‑1 (MQ=255)
                                                                                                  gcgcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCa                                           <  1:212740/147‑1 (MQ=255)
                                                                                                    gcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttctt                                              >  6:274523/1‑142 (MQ=255)
                                                                                                    gcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGttctt                                              <  5:274523/142‑1 (MQ=255)
                                                                                                    gcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAt                                       >  2:244382/1‑149 (MQ=255)
                                                                                                    gcAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAt                                       <  8:91311/149‑1 (MQ=255)
                                                                                                      aaCGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCAtgt                                     >  8:307570/1‑149 (MQ=255)
                                                                                                               gTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATgg                            >  5:226505/1‑149 (MQ=255)
                                                                                                                    aaCTGATCGCCAAGCCTTGAAGTCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGa                       >  8:178372/1‑149 (MQ=255)
                                                                                                                    aaCTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCa                                           <  2:26693/129‑1 (MQ=255)
                                                                                                                    aaCTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCa                                           >  1:26693/1‑129 (MQ=255)
                                                                                                                              cAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCa             <  5:179139/149‑1 (MQ=255)
                                                                                                                               aaGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAg            >  1:307717/1‑149 (MQ=255)
                                                                                                                                  ccTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGtt         >  3:209282/1‑149 (MQ=255)
                                                                                                                                   cTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTg        <  8:93517/149‑1 (MQ=255)
                                                                                                                                    ttGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTgg       >  3:111117/1‑149 (MQ=255)
                                                                                                                                      gAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGTGACTCGCACTACAGGTCGTAATCGATATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTTCGACGTTTGCTCAGGTTgg       <  6:45768/147‑1 (MQ=255)
                                                                                                                                         gCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCttt  <  7:307570/149‑1 (MQ=255)
                                                                                                                                                gccagGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGc                <  4:294348/128‑1 (MQ=255)
                                                                                                                                                gccagGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGc                >  3:294348/1‑128 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
CCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCTTCTGCAGGCGCCGTTCTTCCAGCATGTTGTCGATGGTGCGACGTTTGCTCAGGTTGGTCTTT  >  NC_002947/1249606‑1249891

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