Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,810,366 T→A 6.0% K76M (AAG→ATG)  rpmE ← 50S ribosomal protein L31

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,810,3660TA 6.0% 203.8 / 2.7 66K76M (AAG→ATG) rpmE50S ribosomal protein L31
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base A (2/2);  total (34/32)
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.32e-01

AAAAAAGGCGCCCACAAGGGACGCCTTCTTCAGCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTCGTAGTTCGGGTGGATACCTTCTTTCATGGTCAC  >  NC_002947/5810219‑5810511
                                                                                                                                                   |                                                                                                                                                 
aaaaaaGGCGCCCACAAGGGACGCCTTCTTCAGCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGGACATGCCGAAGCGATCGGCCAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCt                                                                                                                                                   >  5:179001/1‑148 (MQ=255)
 aaaaaGGCGCCCACAAGGGACGCCTTCTTCAGCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTa                                                                                                                                                 >  5:87821/1‑149 (MQ=255)
  aaaaGGCGCCCACAAGGGACGCCTTCTTCAGCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTAc                                                                                                                                                >  7:18038/1‑149 (MQ=255)
          cccACAAGGGACGCCTTCTTCAGCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGgtg                                                                                                                                           >  5:197687/1‑146 (MQ=255)
            cacaAGGGACGCCTTCTTCAGCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGaa                                                                                                                                      >  3:76633/1‑149 (MQ=255)
             acaAGGGACGCCTTCTTCAGCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCATACCGGTGTAGAAc                                                                                                                                     >  3:88170/1‑149 (MQ=255)
                  ggACGCCTTCTTCAGCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTg                                                                                                                                <  6:197687/149‑1 (MQ=255)
                       ccTTCTTCAGCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGt                                                                                                                                 <  5:307080/143‑1 (MQ=255)
                       ccTTCTTCAGCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGt                                                                                                                                 >  6:307080/1‑143 (MQ=255)
                         ttcttcAGCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAgg                                                                                                                         >  5:883/1‑149 (MQ=255)
                              cAGCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACTCGTTGCACACGttaa                                                                                                                     >  4:154204/1‑146 (MQ=255)
                              cAGCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTg                                                                                                                      <  3:154204/147‑1 (MQ=255)
                              cAGCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGc                                                                                                                     <  8:18038/148‑1 (MQ=255)
                                gCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCaaa                                                                                                                  >  7:150861/1‑149 (MQ=255)
                                   acaaTGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGTCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGt                                                                                                               <  6:208254/149‑1 (MQ=255)
                                   acaaTGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGgtgt                                                                                                                                          >  5:182350/1‑122 (MQ=255)
                                   acaaTGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGgtgt                                                                                                                                          <  6:182350/122‑1 (MQ=255)
                                   acaaTGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGt                                                                                                                                 <  7:150713/131‑1 (MQ=255)
                                   acaaTGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGt                                                                                                                                 >  8:150713/1‑131 (MQ=255)
                                    caaTGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGtc                                                                                                              >  6:51724/1‑149 (MQ=255)
                                      aTGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGtcga                                                                                                            >  3:32981/1‑149 (MQ=255)
                                      aTGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGtcga                                                                                                            <  5:51724/149‑1 (MQ=255)
                                       tGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGtcgat                                                                                                           >  8:76046/1‑149 (MQ=255)
                                       tGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGtcgat                                                                                                           <  6:184286/149‑1 (MQ=255)
                                          gAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGgtgt                                                                                                                                          <  6:125148/115‑1 (MQ=255)
                                          gAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGgtgt                                                                                                                                          >  5:125148/1‑115 (MQ=255)
                                           aaaCCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGAc                                                                                                       <  6:81846/149‑1 (MQ=255)
                                            aaCCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGtcgatcga                                                                                                        >  2:140363/1‑147 (MQ=255)
                                                cGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGAc                                                                                                  >  5:284923/1‑149 (MQ=255)
                                                       ttCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAAAGTCGATCGACAGGACGTTGcc                                                                                            >  2:73166/1‑148 (MQ=255)
                                                              tGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGCCAGGACGTTGCCCAGGGTc                                                                                     >  4:232771/1‑148 (MQ=255)
                                                                cgcATGATTACTTCTTGGCACCGGACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGaa                                                                                  <  6:179001/149‑1 (MQ=255)
                                                                   aTGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAAc                                                                                 >  1:289879/1‑147 (MQ=255)
                                                                    tGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAg                                                                              <  7:292981/149‑1 (MQ=255)
                                                                       ttACTTCTTTGCACCGAACATTTCGAAGCGATCGGCGAACTTCTGTACGCGAGCACCGGTGTCCAGATCTTTCTTGATACCGGTGTATAAAGGTTGGGACTGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGtt                                                                            <  4:88170/148‑1 (MQ=255)
                                                                          cttcttGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTc                                                                                     <  6:191410/136‑1 (MQ=255)
                                                                             cttGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGaa                                                                       <  7:76046/147‑1 (MQ=255)
                                                                             cttGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGaa                                                                       <  6:87821/147‑1 (MQ=255)
                                                                                ggCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAg                                                                                                                          >  1:3758/1‑93 (MQ=255)
                                                                                ggCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAg                                                                                                                          <  2:3758/93‑1 (MQ=255)
                                                                                ggCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGt                                                                  >  8:97958/1‑149 (MQ=255)
                                                                                   aCCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTAc                                                               <  1:140363/149‑1 (MQ=255)
                                                                                       aaCATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCg                                                             >  1:900/1‑147 (MQ=255)
                                                                                        aCATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAg                                                          >  5:202949/1‑149 (MQ=255)
                                                                                                      tCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGCGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGCACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTg                                             >  2:300274/1‑148 (MQ=255)
                                                                                                        ggCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCaa                                          <  4:76633/149‑1 (MQ=255)
                                                                                                               ttCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTc                                   <  6:284923/149‑1 (MQ=255)
                                                                                                               ttCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTc                                   >  7:186189/1‑149 (MQ=255)
                                                                                                                 cTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTCg                                  <  3:3283/148‑1 (MQ=255)
                                                                                                                    tACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTCGTAGt                              >  3:305040/1‑149 (MQ=255)
                                                                                                                           ccaccGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTCGTAGTTCGGGTg                       <  8:150861/149‑1 (MQ=255)
                                                                                                                           acaccGGTAAGCAGAAAGTTCTGGACACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTCg                                  >  5:166260/2‑138 (MQ=38)
                                                                                                                           acaccGGTAAGCAGAAAGTTCTGGACACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTCg                                  <  6:166260/137‑1 (MQ=38)
                                                                                                                              ccGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGa                                                                               >  3:31498/1‑90 (MQ=255)
                                                                                                                              ccGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGa                                                                               <  4:31498/90‑1 (MQ=255)
                                                                                                                              ccGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTCGTAGTTCGGGTGGAt                    <  8:186189/149‑1 (MQ=255)
                                                                                                                                  tgtCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTCGTAGTTCGGGTGGATACct                <  4:305040/149‑1 (MQ=255)
                                                                                                                                      cAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTCGTAGTTCGGGTGGATAcc                 <  5:62055/144‑1 (MQ=255)
                                                                                                                                      cAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTCGTAGTTCGGGTGGATAcc                 >  6:62055/1‑144 (MQ=255)
                                                                                                                                       aGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCaa                                          <  4:166720/118‑1 (MQ=255)
                                                                                                                                       aGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCaa                                          >  3:166720/1‑118 (MQ=255)
                                                                                                                                       aGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAac                                         <  3:102655/119‑1 (MQ=255)
                                                                                                                                       aGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAac                                         >  4:102655/1‑119 (MQ=255)
                                                                                                                                        gAACTTTCTGCTTACCGGTGTAGAACGCGTGGCACAGCATACAAACCTCGATCGACACTACGTAGCCCAGGCTCGAACGAGTTACGAACTTCTTACCGCAGCAGCAGCTCACTGCAACAACTACGTAGTACGGGTGGATAACATCttt           >  8:141112/1‑148 (MQ=255)
                                                                                                                                           cTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTCGTAGTTCGGGTGGATACCTTCTTTCATg       <  7:97958/149‑1 (MQ=255)
                                                                                                                                             ttCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTCGTAGTTCGGGTGGATACCTTCTTTCATgg      >  5:210461/1‑148 (MQ=255)
                                                                                                                                                  cTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTCGTAGTTCGGGTGGATACCTTCTTTCATGGTCAc  <  6:202949/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AAAAAAGGCGCCCACAAGGGACGCCTTCTTCAGCAACAATGCGAAACCCGAAGGATTCGCCATGCGCATGATTACTTCTTGGCACCGAACATGCCGAAGCGATCGGCGAACTTCTGTACGCGACCACCGGTGTCCAGAACTTTCTGCTTACCGGTGTAGAACGGGTGGCACAGGTTGCAAACGTCGATCGACAGGACGTTGCCCAGGGTCGAACGAGTTTCGAACTTGTTACCGCAGCTGCAGGTGACTGCAACTACTTCGTAGTTCGGGTGGATACCTTCTTTCATGGTCAC  >  NC_002947/5810219‑5810511

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