Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,421,572 T→C 8.6% L128P (CTG→CCG)  moaB‑I → molybdenum cofactor biosynthesis protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,421,5720TC 8.6% 159.1 / 4.9 59L128P (CTG→CCG) moaB‑Imolybdenum cofactor biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base T (27/26);  new base C (2/3);  total (30/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 = 1.31e-01

GTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGCCTGTGA  >  NC_002947/2421425‑2421716
                                                                                                                                                    |                                                                                                                                                 
gTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTg                                                                                                                                                  <  4:199906/149‑1 (MQ=255)
 tACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTgg                                                                                                                                                 <  6:21826/149‑1 (MQ=255)
          tACTGGTCGCGACAGCACCCCCGCGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCGTTGGCCAGCCCGGCCaaggc                                                                                                                                          <  5:34388/148‑4 (MQ=255)
             tGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGc                                                                                                                                     >  3:114647/1‑149 (MQ=255)
               gTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCGTTGGCCAGCCCGGCCAAGGCACGcga                                                                                                                                   <  1:302963/149‑3 (MQ=255)
                 cgcgACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGgt                                                                                                                                 <  1:344991/149‑1 (MQ=255)
                  gcgACAGCACCCCCGAGGCGGTGGACTGCCTGCTGGACAATCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGAACggttgtt                                                                                                                                >  7:57737/1‑143 (MQ=255)
                    gACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCTAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGtgtg                                                                                                                              >  5:374312/1‑149 (MQ=255)
                          aCCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGctgcct                                                                                                                        >  7:5216/1‑149 (MQ=255)
                            ccccGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGctgcctgc                                                                                                                      >  1:326330/1‑149 (MQ=255)
                             cccGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCGTTGGCCAGCCCGGCCaag                                                                                                                                            >  4:158030/1‑126 (MQ=255)
                              ccGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCa                                                                                                                    <  8:55818/149‑1 (MQ=255)
                              ccGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAAAGGGCCGCTGGTGTGctgcctgc                                                                                                                      >  3:195895/1‑147 (MQ=255)
                                     ggTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTcca                                                                                                             >  5:213650/1‑149 (MQ=255)
                                        ggCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGGCCGCTGGTGTGGTGCCTGCCCGGCTccacc                                                                                                           >  1:383089/1‑148 (MQ=255)
                                           ctgcctgcTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCGTTGGCCAGCCCGGCCaaggca                                                                                                                                         >  2:209379/1‑112 (MQ=38)
                                           ctgcctgcTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCGTTGGCCAGCCCGGCCaaggca                                                                                                                                         <  1:209379/116‑5 (MQ=38)
                                             gcctgcTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCt                                                                                                     >  4:231679/1‑149 (MQ=255)
                                                 gctgGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGc                                                                                                 >  4:98209/1‑149 (MQ=255)
                                                 gctgGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGc                                                                                                 >  3:108052/1‑149 (MQ=255)
                                                     gACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTg                                                                                                                                                  <  2:27396/96‑1 (MQ=255)
                                                     gACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTg                                                                                                                                                  >  1:27396/1‑96 (MQ=255)
                                                     gACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCATGCTCCACCAACGCTGTGCGTAc                                                                                             >  2:239925/1‑149 (MQ=255)
                                                     gACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTAc                                                                                             >  1:231193/1‑149 (MQ=255)
                                                     gACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTAc                                                                                             >  4:81055/1‑149 (MQ=255)
                                                          gcaagTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGtt                                                                                        >  3:311270/1‑149 (MQ=255)
                                                              gTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGa                                                                                    >  8:222836/1‑149 (MQ=255)
                                                                 gACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGGCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATg                                                                                  <  5:85495/148‑1 (MQ=255)
                                                                 gACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATg                                                                                  <  7:222836/148‑1 (MQ=255)
                                                                     gCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATc                                                                             <  2:326330/149‑1 (MQ=255)
                                                                       tttGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCt                                                                           >  5:93215/1‑149 (MQ=255)
                                                                           gTGAGCTGTCCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGc                                                                        <  6:213650/148‑1 (MQ=255)
                                                                            tGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCg                                                                      <  4:108052/149‑1 (MQ=255)
                                                                              aGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGag                                                                    <  5:210384/149‑1 (MQ=255)
                                                                              aGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGa                                                                     <  4:124599/148‑1 (MQ=255)
                                                                                   ttCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGagcagc                                                                <  2:385297/148‑1 (MQ=255)
                                                                                   ttCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTGTGCCGGGCGGGCCAAACGCACCCGGGTGTGCTTCCAGCCCCGCTCCACCAACGCTGTGGGTTCCGGGTGGGGTTGCCTGCTGGCCGCGCCGCt                                                               >  1:300476/1‑149 (MQ=255)
                                                                                     ccGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTc                                                              <  3:231679/148‑1 (MQ=255)
                                                                                      cGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGa                                                            >  3:67840/1‑149 (MQ=255)
                                                                                        gCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATg                                                         <  6:374312/149‑1 (MQ=255)
                                                                                           gATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGTCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGccc                                                      <  4:311270/149‑1 (MQ=255)
                                                                                             tttCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGt                                                    <  2:130834/149‑1 (MQ=255)
                                                                                               tCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTc                                                   <  4:114647/148‑1 (MQ=255)
                                                                                                ccGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCAt                                                 >  3:315297/1‑149 (MQ=255)
                                                                                                      ccGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCAGCCATGCTCCACCAACCCTTTGCGTACCGGTTGGGATTGCATCCAGGCCGACAAGCACGAAGCCCGTCATCGccca                                           >  7:232975/1‑148 (MQ=255)
                                                                                                               gCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTc                                  <  4:195895/149‑1 (MQ=255)
                                                                                                                   cTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTgg                              <  1:239925/149‑1 (MQ=255)
                                                                                                                      cACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGgcg                            >  6:327271/1‑148 (MQ=255)
                                                                                                                          gTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCt                       >  5:35296/1‑149 (MQ=255)
                                                                                                                             cAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGGTCCACCAACGCTGTGGGTTCCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTg                      >  6:287393/1‑147 (MQ=255)
                                                                                                                                     tGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAgcaggcag            >  1:359623/1‑149 (MQ=255)
                                                                                                                                      gCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGc           <  3:81055/149‑1 (MQ=255)
                                                                                                                                       ccTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATCGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGc           <  1:280829/149‑1 (MQ=255)
                                                                                                                                       ccTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCg          <  5:327271/149‑1 (MQ=255)
                                                                                                                                       ccTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCg          <  2:231193/149‑1 (MQ=255)
                                                                                                                                        cTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCgg         >  6:134092/1‑149 (MQ=255)
                                                                                                                                        cTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCgg         >  3:116310/1‑149 (MQ=255)
                                                                                                                                            ggCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGc        <  2:359623/147‑1 (MQ=255)
                                                                                                                                                gggCTGGCCAACGTCACGCTGGTGTGCTGCCTGCCATGATCCACCAACGCTGTGCGTACCGGTTTGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTTATTGCTCTTGCAACTTTGTGGTGCATCTGAAGTATGTAGCGGCCTGTGa  <  8:57737/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GTACCGGGTTTACTGGTCGCGACAGCACCCCCGAGGCGGTGGCCTGCCTGCTGGACAAGCAAGTAGACGGCTTTGGTGAGCTGTTCCGGCAGATTTCCGTGGCCGACATTGGCACCTCCACCGTGCAGTCGCGTGCCTTGGCCGGGCTGGCCAACGGCACGCTGGTGTGCTGCCTGCCAGGCTCCACCAACGCTGTGCGTACCGGTTGGGATGGCATCCTGGCCGAGCAGCTCGATGCCCGTCATCGCCCTTGCAACTTCGTGGCGCACCTGAAGCAGGCAGCGGCCTGTGA  >  NC_002947/2421425‑2421716

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