Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 494,056 T→A 5.7% S164T (TCC→ACC)  PP_0406 → nucleotidyltransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947494,0560TA 5.7% 218.0 / 2.1 70S164T (TCC→ACC) PP_0406nucleotidyltransferase family protein
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base A (2/2);  total (36/34)
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.73e-01

CGCGTTTGCAGGCTCCCCTGCAAGGCCTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTGGGCACTGGGTGGATGTTGGCACGCTTGAACGCCT  >  NC_002947/493910‑494201
                                                                                                                                                  |                                                                                                                                                 
cgcgTTTGCAGGCTCCCCTGCAAGGCCTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGCCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTcc                                                                                                                                                 >  5:188998/1‑149 (MQ=255)
 gcgTTTGCAGGCTCCCCTGCAAGGCCTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCg                                                                                                                                                >  7:40781/1‑149 (MQ=255)
    tttGCAGGCTCCCCTGCAAGGCCTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCTGCGGCATTTCCGtgc                                                                                                                                             <  7:228523/149‑1 (MQ=255)
    tttGCAGGCTCCCCTGCAAGGCCTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGtgc                                                                                                                                             <  3:275805/149‑1 (MQ=255)
             tCCCCTGCAAGGCCTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGATGACCTTCAGCGGCATTTCCGTGCTGCACCCgc                                                                                                                                    >  6:150861/1‑148 (MQ=255)
               cccTGCAAGGCCTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGgcg                                                                                                                                  >  2:327981/1‑149 (MQ=255)
                ccTGCAAGGCCTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGTCATTTCCGtgc                                                                                                                                             >  5:128314/1‑137 (MQ=255)
                ccTGCAAGGCCTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGtgc                                                                                                                                             <  6:128314/137‑1 (MQ=255)
                ccTGCAAGGCCTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGgcgc                                                                                                                                 <  6:267118/149‑1 (MQ=255)
                 cTGCAAGGCCTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTACCGTGCTGCACCCGGCGCt                                                                                                                                >  4:44676/1‑149 (MQ=255)
                 cTGCAAGGCCTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTACCGTGCTGCACCCGGCGCt                                                                                                                                >  4:40474/1‑149 (MQ=255)
                     aaGGCCTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTc                                                                                                                            >  1:31487/1‑149 (MQ=255)
                        gCCTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCg                                                                                                                           >  8:1610/1‑147 (MQ=255)
                          cTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAgg                                                                                                                       >  4:187054/1‑149 (MQ=255)
                             gCTCACCTGGTACTGGTCGACAACCCTGCCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGtt                                                                                                                     <  5:292424/148‑1 (MQ=255)
                                 aCCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAg                                                                                                                <  8:40781/149‑1 (MQ=255)
                                 aCCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAg                                                                                                                >  5:132817/1‑149 (MQ=255)
                                 aCCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAg                                                                                                                >  3:35351/1‑149 (MQ=255)
                                      gTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTACCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCggg                                                                                                           <  3:44676/149‑1 (MQ=255)
                                      gTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTACCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCggg                                                                                                           <  3:40474/149‑1 (MQ=255)
                                       tACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGt                                                                                                          >  4:250916/1‑149 (MQ=255)
                                          tGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGcc                                                                                                       >  5:52278/1‑149 (MQ=255)
                                            gTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCtt                                                                                                     <  1:77401/149‑1 (MQ=255)
                                            gTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCtt                                                                                                     >  2:77401/1‑149 (MQ=255)
                                            gTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCt                                                                                                      <  2:31487/148‑1 (MQ=255)
                                            gTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCt                                                                                                      <  6:252303/148‑1 (MQ=255)
                                            gTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCt                                                                                                      <  3:187054/148‑1 (MQ=255)
                                                 cAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAAc                                                                                                <  1:47564/149‑1 (MQ=255)
                                                         gCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCAc                                                                                          <  3:250916/147‑1 (MQ=255)
                                                          ccATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGt                                                                                       <  6:132817/149‑1 (MQ=255)
                                                           cATCTCGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGtt                                                                                      >  5:122937/1‑149 (MQ=255)
                                                           cATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGtt                                                                                      <  4:96026/149‑1 (MQ=255)
                                                             tCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTtgc                                                                                    <  1:327981/149‑1 (MQ=255)
                                                              cACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCtt                                                                                                     <  6:328971/131‑1 (MQ=255)
                                                              cACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCtt                                                                                                     >  5:328971/1‑131 (MQ=255)
                                                                   ccGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTgcgc                                                                               >  3:262924/1‑148 (MQ=255)
                                                                   ccGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTgcgc                                                                               >  7:72520/1‑148 (MQ=255)
                                                                       ggtgACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCg                                                                                                                                                <  4:316386/79‑1 (MQ=255)
                                                                       ggtgACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCg                                                                                                                                                >  3:316386/1‑79 (MQ=255)
                                                                            cTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATgg                                                                     <  6:188998/149‑1 (MQ=255)
                                                                             ttCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATggc                                                                    >  7:259151/1‑149 (MQ=255)
                                                                                cGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATggcggc                                                                 >  1:23659/1‑149 (MQ=255)
                                                                                  ccTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGcggccg                                                               <  8:72520/149‑1 (MQ=255)
                                                                                      gtggGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCaa                                                           <  6:52278/149‑1 (MQ=255)
                                                                                      gtggACGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGcggccggc                                                             <  7:1610/147‑1 (MQ=255)
                                                                                          gCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTa                                                       >  3:311907/1‑149 (MQ=255)
                                                                                              gCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTg                                                   >  2:307098/1‑149 (MQ=255)
                                                                                                 ggTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGc                                                 >  2:77278/1‑148 (MQ=255)
                                                                                                     gTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAaca                                            >  6:180706/1‑149 (MQ=255)
                                                                                                            gCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGt                                     <  6:122937/149‑1 (MQ=255)
                                                                                                             cGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTg                                    <  4:311907/149‑1 (MQ=255)
                                                                                                                     cgcCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTggg                                  >  3:153197/1‑143 (MQ=255)
                                                                                                                     cgcCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTggg                                  <  4:153197/143‑1 (MQ=255)
                                                                                                                      gcCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTGGGCACTggg                           <  8:259151/149‑1 (MQ=255)
                                                                                                                        ccGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGaa                                              <  1:56802/128‑1 (MQ=255)
                                                                                                                        ccGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGaa                                              >  2:56802/1‑128 (MQ=255)
                                                                                                                        ccGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTGGGCACTGGGTg                         >  3:92940/1‑149 (MQ=255)
                                                                                                                        ccGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTGGGCACTGGGTg                         <  5:180706/149‑1 (MQ=255)
                                                                                                                        ccGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGCTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTGGGCACTGGGTg                         >  8:278926/1‑149 (MQ=255)
                                                                                                                          ggCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTGGGCACTGGGTg                         <  4:250092/147‑1 (MQ=255)
                                                                                                                          ggCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTGGGCACTGGGTg                         >  3:250092/1‑147 (MQ=255)
                                                                                                                           gCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTGGGCACTGGGTGGAt                      <  1:307098/149‑1 (MQ=255)
                                                                                                                           gCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTGGGCACTGGGTGGAt                      <  1:77278/149‑1 (MQ=255)
                                                                                                                               gCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGaa                                              >  6:212098/1‑121 (MQ=255)
                                                                                                                               gCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGaa                                              <  5:212098/121‑1 (MQ=255)
                                                                                                                                cTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTGGGCACTGGGTGGATGTTgg                 >  2:280086/1‑149 (MQ=255)
                                                                                                                                     ctacAGCGGCATTTCCGAGCTGCACCCGGAGCAGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCCCCGTTGCTGCGGCAGGCGATGGCGGCCGGCAAGGTATCTGGCGCACACAACCGTGGGCAATGGGTGGATGTTGGCTCGc            >  6:268134/4‑149 (MQ=255)
                                                                                                                                             gCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACAGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTGGGCACTGGGTGGATGTTGGCACGCTTGAACGc    >  4:29679/1‑149 (MQ=255)
                                                                                                                                               aTTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTGGGCACTGGGTGGATGTTGGCACGCTTGAACGCCt  <  4:35351/149‑1 (MQ=255)
                                                                                                                                               aTTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTGGGCACTGGGTGGATGTTGGCACGCTTGAACGCCt  >  3:196361/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGCGTTTGCAGGCTCCCCTGCAAGGCCTGGCTCACCTGGTACTGGTCGACAACCCTGGCCATCACGGCCGTGGTGACTTCCGCCTGGTGGGCGAGCAGGTCGTTGACGGCGATGACGCGCCCGGCACGCTGACCTTCAGCGGCATTTCCGTGCTGCACCCGGCGCTGTTCGAAGGTTGCCAGGCGGGTGCCTTCAAACTGGCACCGTTGCTGCGCCAGGCGATGGCGGCCGGCAAGGTATCTGGCGAACACTACCGTGGGCACTGGGTGGATGTTGGCACGCTTGAACGCCT  >  NC_002947/493910‑494201

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