Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,687,086 C→T 6.3% D173N (GAC→AAC)  ydcV ← polyamine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,687,0860CT 6.3% 191.6 / 2.2 63D173N (GAC→AAC) ydcVpolyamine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (32/27);  new base T (2/2);  total (34/29)
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.25e-01

AGATCGGCAGGGTGCGCTCGTGGCCGGCGGTGAACGTGGTGACGATGATTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGCCGACTACGATGGTGAACACCCCGGGCTCTATCCCCAGCGTCTTGAA  >  NC_002947/1686942‑1687233
                                                                                                                                                |                                                                                                                                                   
aGATCGGCAGGGTGCGCTCGTGGCCGGCGGTGAACGTGGTGACGATGATTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGgcg                                                                                                                                                 >  8:86454/1‑149 (MQ=255)
aGATCGGCAGGGTGCGCTCGTGGCCGGCGGTGAACGTGGTGACGATGATTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGgcg                                                                                                                                                 >  6:34346/1‑149 (MQ=255)
       cAGGGTGCGCTCGTGGCCGGCGGTGAACGTGGTGACGATGATTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGt                                                                                                                                          <  5:335012/149‑1 (MQ=255)
          ggTGCGCTCGTGGCCGGCGGTGAACGTGGTGACGATGATTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTTGGCGCCAAGGTCCa                                                                                                                                       >  1:248518/1‑149 (MQ=255)
          ggTGCGCTCGTGGCCGGCGGTGAACGTGGTGACGATGATTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTTGGCGCCAAGGTCCa                                                                                                                                       >  5:55922/1‑149 (MQ=255)
           gTGCGCTCGTGGCCGGCGGTGAACGTGGTGACGATGATTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCAt                                                                                                                                      >  4:333620/1‑149 (MQ=255)
            tGCGCTCGTGGCCGGCGGTGAACGTGGTGACGATGATTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATc                                                                                                                                     <  7:86454/149‑1 (MQ=255)
                cTCGTGGCCGGCGGTGAACGTGGTGACGATGATTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAgg                                                                                                                                 <  5:34346/149‑1 (MQ=255)
                    tGGCCGGCGGTGAACGTGGTGACGATGATTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTTGGCGCCAAGGTCCATCGAGGCTTc                                                                                                                             <  6:55922/149‑1 (MQ=255)
                    tGGCCGGCGGTGAACGTGGTGACGATGATTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTTGGCGCCAAGGTCCATCGAGGCTTc                                                                                                                             <  2:248518/149‑1 (MQ=255)
                         ggcggTGAACGTGGTGACGATGATTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCa                                                                                                                        <  2:62043/149‑1 (MQ=255)
                                           gatgatTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAgg                                                                                                                                           <  4:131658/112‑1 (MQ=255)
                                           gatgatTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAgg                                                                                                                                           >  3:131658/1‑112 (MQ=255)
                                           gatgatTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAg                                                                                                             <  5:211324/142‑1 (MQ=255)
                                           gatgatTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAg                                                                                                             >  6:211324/1‑142 (MQ=255)
                                            atgatTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTAcgccg                                                                                                     >  7:305589/1‑149 (MQ=255)
                                            atgatTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAAGcct                                                                                                                              <  4:130187/124‑3 (MQ=255)
                                            atgatTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAAGcct                                                                                                                              >  3:130187/1‑122 (MQ=255)
                                              gatTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTAcgccg                                                                                                     >  6:120914/1‑147 (MQ=255)
                                                tttCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCTGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAgg                                                                                                 >  1:18277/1‑149 (MQ=255)
                                                    aTCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAggcggg                                                                                             >  2:14683/1‑149 (MQ=255)
                                                    aTCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAAGCTTTGCGAGGTACGCCGCAggcggg                                                                                             >  6:30528/1‑149 (MQ=255)
                                                     tCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCa                                                                                                   <  5:386241/142‑1 (MQ=255)
                                                     tCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCa                                                                                                   >  6:386241/1‑142 (MQ=255)
                                                     tCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAggcgggc                                                                                            >  7:49255/1‑149 (MQ=255)
                                                      cGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAggcgggc                                                                                            >  5:193360/1‑148 (MQ=255)
                                                        aaCGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCAt                                                                                                                                      >  6:170824/1‑104 (MQ=255)
                                                        aaCGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCAt                                                                                                                                      <  5:170824/104‑1 (MQ=255)
                                                         aCGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATc                                                                                        <  3:333620/149‑1 (MQ=255)
                                                           gACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATca                                                                                       >  6:326930/1‑148 (MQ=255)
                                                             cAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTTCGAGGTACGCCGCAGGCGGGCGATCACAt                                                                                    >  3:216515/1‑149 (MQ=255)
                                                                ggCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACAttgt                                                                                 >  3:406145/1‑149 (MQ=255)
                                                                 gCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCTGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACAttgtt                                                                                <  2:18277/149‑1 (MQ=255)
                                                                 gCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACAttgtt                                                                                <  8:49255/149‑1 (MQ=255)
                                                                  cAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACACCGCAGGCGGGCGATCACAttgttg                                                                               >  5:59424/1‑149 (MQ=255)
                                                                         cAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACg                                                                         >  8:144830/1‑148 (MQ=255)
                                                                         cAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGa                                                                        >  5:343242/1‑149 (MQ=255)
                                                                          aGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGAt                                                                       >  1:70175/1‑149 (MQ=255)
                                                                          aGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGAt                                                                       <  8:305589/149‑1 (MQ=255)
                                                                           gCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATc                                                                      <  2:70175/149‑1 (MQ=255)
                                                                            cATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATca                                                                     <  1:242131/149‑1 (MQ=255)
                                                                            cATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATca                                                                     >  2:165153/1‑149 (MQ=255)
                                                                                ccgcccgcCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACACCGCAGGCGGGCGATCACATTGTTGTAGACGATcacca                                                                  <  6:59424/148‑1 (MQ=255)
                                                                                 cgcccgccAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACg                                                                <  5:120914/149‑1 (MQ=255)
                                                                                   cccgccAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACg                                                                >  5:7793/1‑147 (MQ=255)
                                                                                          gcagcGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAAcgtg                                                       <  1:14683/149‑1 (MQ=255)
                                                                                                  gagccgagGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGcc                                               <  3:285319/149‑1 (MQ=255)
                                                                                                  gagccgagGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGcc                                               <  4:406145/149‑1 (MQ=255)
                                                                                                       gagGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGCCGACTa                                          >  2:79793/1‑149 (MQ=255)
                                                                                                       gagGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGCCGACTa                                          <  6:7793/149‑1 (MQ=255)
                                                                                                          gTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGCCGACTACGa                                       >  4:338559/1‑149 (MQ=255)
                                                                                                             gggCAGGACGATGTAGCGGAAGGTCTGCCAACGGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGCCGACTACGATg                                     <  4:216515/148‑1 (MQ=255)
                                                                                                              ggCAGGACGAAGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGCTCACCACGCAGAACGTGGCGTGGCCGACTACGATGGt                                   <  5:326930/149‑1 (MQ=255)
                                                                                                               gCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGCCGACTACGATGGt                                   >  3:349347/1‑148 (MQ=255)
                                                                                                                   gACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGCCGACTACGATGGTGAaca                              >  6:322179/1‑149 (MQ=255)
                                                                                                                   gACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGCCGACTACGATGGTGAaca                              >  8:312160/1‑149 (MQ=255)
                                                                                                                                     cTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGCCGACTACGATGGTGAACACCCCGGGCTCTATCCCCa            <  6:343242/149‑1 (MQ=255)
                                                                                                                                      tGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACACAGAACGTGGCGTGGCCGACTACGATGGTGAACACCCCGGGCTCTATCCCCAg           >  1:133391/1‑149 (MQ=255)
                                                                                                                                       gCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGCCGACTACGATGGTGAACACCCCGGGCTCTATCCCCAGc          >  4:110030/1‑149 (MQ=255)
                                                                                                                                            ccGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGCCGACTACGATGGTGAACACCCCGGGCTCTATCCCCAGCGTCtt     <  3:338559/149‑1 (MQ=255)
                                                                                                                                              gTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGCCGACt                                           >  7:142391/1‑109 (MQ=255)
                                                                                                                                              gTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGCCGACt                                           <  8:142391/109‑1 (MQ=255)
                                                                                                                                               tCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGCCGACTACGATGGTGAACACCCCGGGCTCTATCCCCAGCGTCTTGaa  <  5:322179/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AGATCGGCAGGGTGCGCTCGTGGCCGGCGGTGAACGTGGTGACGATGATTTCATCGAACGACAAGGCAAAGGCCAGCATGCCGCCCGCCAGCAGCGCCGAGCCGAGGTTGGGCAGGACGATGTAGCGGAAGGTCTGCCAACCGTCGGCGCCAAGGTCCATCGAGGCTTCGATCAGGCTTTGCGAGGTACGCCGCAGGCGGGCGATCACATTGTTGTAGACGATCACCACGCAGAACGTGGCGTGGCCGACTACGATGGTGAACACCCCGGGCTCTATCCCCAGCGTCTTGAA  >  NC_002947/1686942‑1687233

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: