Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,124,124 T→C 5.5% H3R (CAT→CGT)  PP_4510 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,124,1240TC 5.5% 220.5 / 2.0 73H3R (CAT→CGT) PP_4510hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/33);  new base C (2/2);  total (38/35)
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.75e-01

GCCCACCAGATAGACGCGATCAGCCTGGTAGCCCTTGAGAAACTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGCCTGCCCTGCGAGGAAGCCCCCATGACCGCCGGCC  >  NC_002947/5123979‑5124271
                                                                                                                                                  |                                                                                                                                                   
gCCCACCAGATAGACGCGATCAGCCTGGTAGCCCTTGAGAAACTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATgg                                                                                                                                                   <  3:325176/148‑1 (MQ=255)
 cccACCAGATAGACGCGATCAGCCTGGTAGCCCTTGAGAAACTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGc                                                                                                                                                  >  5:180211/1‑148 (MQ=255)
                     aGCCTGGTAGCCCTTGAGAAACTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGgcgc                                                                                                                              <  1:96895/148‑1 (MQ=255)
                      gCCTGGTAGCCCTTGAGAAACTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAg                                                                                                                            >  2:316690/1‑149 (MQ=255)
                      gCCTGGTAGCCCTTGAGAAACTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTAGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAg                                                                                                                            >  3:259040/1‑149 (MQ=255)
                               cccTTGAGAAACTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTg                                                                                                                     >  8:89657/1‑147 (MQ=255)
                               cccTTGAGAAACTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGga                                                                                                                   <  6:204051/149‑1 (MQ=255)
                                  ttGAGAAACTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGAc                                                                                                                >  6:2083/1‑149 (MQ=255)
                                    gagaAACTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGAACGTCGGAGAT‑CCACTGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACt                                                                                                               <  7:248476/148‑1 (MQ=255)
                                      gaAACTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTgg                                                                                                             >  6:156913/1‑148 (MQ=255)
                                      gaAACTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTgg                                                                                                             >  4:189115/1‑148 (MQ=255)
                                       aaaCTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTggca                                                                                                           <  5:156913/149‑1 (MQ=255)
                                       aaaCTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTgg                                                                                                             >  6:162831/1‑147 (MQ=255)
                                         aCTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGga                                                                                                                   <  4:71134/139‑1 (MQ=255)
                                         aCTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGga                                                                                                                   >  3:71134/1‑139 (MQ=255)
                                         aCTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTggcagg                                                                                                         <  6:180211/149‑1 (MQ=255)
                                          cTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTggcaggc                                                                                                        <  1:316690/149‑1 (MQ=255)
                                          cTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTggcaggc                                                                                                        <  4:12837/149‑1 (MQ=255)
                                             cGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGcac                                                                                                      <  5:2083/148‑1 (MQ=255)
                                                   gtgtTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGCCTGGCAGGCACACATGa                                                                                                >  5:280425/1‑148 (MQ=255)
                                                   gtgtTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGa                                                                                                >  5:287894/1‑148 (MQ=255)
                                                   gtgtTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCACGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGTCACACAt                                                                                                  <  3:263837/146‑1 (MQ=255)
                                                   gtgtTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCACGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACAt                                                                                                  >  4:263837/1‑146 (MQ=255)
                                                   gtgtTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCACGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACAt                                                                                                  >  4:259325/1‑146 (MQ=255)
                                                   gtgtTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCACGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACAt                                                                                                  <  3:259325/146‑1 (MQ=255)
                                                      ttCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAgc                                                                                            >  1:199684/1‑149 (MQ=255)
                                                           ccTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCAt                                                                                       <  1:273542/149‑1 (MQ=255)
                                                             tGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATgg                                                                                     >  1:280774/1‑149 (MQ=255)
                                                             tGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATgg                                                                                     >  5:153599/1‑149 (MQ=255)
                                                               gCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGc                                                                                    <  7:89657/148‑1 (MQ=255)
                                                                   tCGCGGGTACCCAGGTGCACGTCGGAGGTGCCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCaaaa                                                                                >  1:198384/1‑149 (MQ=255)
                                                                   tCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAAc                                                                               >  8:275518/1‑149 (MQ=255)
                                                                   tCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGACGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAAc                                                                               >  1:266995/1‑149 (MQ=255)
                                                                    cgcgGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAATCg                                                                              <  6:287894/149‑1 (MQ=255)
                                                                         gTACCCAGGGGCACTTCGGAGAT‑CCACAGGGTACGCACGCGTTGGTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGAc                                                                         <  2:199684/149‑1 (MQ=255)
                                                                             ccAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGc                                                                                                                       <  8:275331/99‑1 (MQ=255)
                                                                             ccAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGc                                                                                                                       >  7:275331/1‑99 (MQ=255)
                                                                             ccAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGt                                                                     >  4:300413/1‑149 (MQ=255)
                                                                             ccAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGt                                                                     >  6:50746/1‑149 (MQ=255)
                                                                                    cACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGc                                                               <  6:153599/148‑1 (MQ=255)
                                                                                    cACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGGGACGAGTGATGGCg                                                              >  6:169647/1‑149 (MQ=255)
                                                                                          ggAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAg                                                         <  3:189115/148‑1 (MQ=255)
                                                                                           gagaT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAgg                                                        >  4:226854/1‑148 (MQ=255)
                                                                                              tt‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGAtt                                                    <  2:198384/148‑1 (MQ=255)
                                                                                              aT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGAt                                                     <  7:275518/148‑1 (MQ=255)
                                                                                                  cacaGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGa                                                  <  5:162831/148‑1 (MQ=255)
                                                                                                  cacaGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGa                                                  <  5:239540/148‑1 (MQ=255)
                                                                                                  cacaGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAg                                                 >  3:164301/1‑149 (MQ=255)
                                                                                                  cacaGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAg                                                 <  4:259040/149‑1 (MQ=255)
                                                                                                   acaGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAgg                                                >  7:257977/1‑149 (MQ=255)
                                                                                                       ggTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGcacac                                                                                                    <  4:62295/93‑1 (MQ=255)
                                                                                                           cgcacgcGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATg                                         <  5:50746/148‑1 (MQ=255)
                                                                                                              acgcGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGgcg                                     >  7:302263/1‑149 (MQ=255)
                                                                                                                gcgTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGcc                                   >  2:49383/1‑149 (MQ=255)
                                                                                                                 cgTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGcct                                  <  2:280774/149‑1 (MQ=255)
                                                                                                                 cgTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGcct                                  <  4:164301/149‑1 (MQ=255)
                                                                                                                     gcttgcGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGcctgccc                              >  6:229421/1‑149 (MQ=255)
                                                                                                                      atctcgcgAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGcctgccct                             <  3:273927/145‑1 (MQ=255)
                                                                                                                       ttgcGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGcctgccctg                            >  2:107609/1‑149 (MQ=255)
                                                                                                                         gcgcGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACg                                                                              <  2:57364/97‑1 (MQ=255)
                                                                                                                          cgcgAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGCCTGCCCTGCGa                         <  3:300413/149‑1 (MQ=255)
                                                                                                                           gcgAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGcctgccctgc                           >  1:219456/1‑146 (MQ=255)
                                                                                                                           gcgAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGcctgccct                             >  4:273927/1‑144 (MQ=255)
                                                                                                                           gcgAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGCCTGCCCTGCGa                         <  2:219456/148‑1 (MQ=255)
                                                                                                                            cgAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGCCTGCCCTGCGAgg                       <  2:243528/149‑1 (MQ=255)
                                                                                                                               ggggCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGCCTGCCCTGCGAGGAAg                    >  7:194912/1‑149 (MQ=255)
                                                                                                                                  gCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGCCTGCCCTGCGAGGAAGccc                 <  3:226854/149‑1 (MQ=255)
                                                                                                                                      gACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGCCTGCCCTGCGAGGAAGCCCCCAt             >  7:27971/1‑149 (MQ=255)
                                                                                                                                      gACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGCCTGCCCTGCGAGGAAGCCCCCAt             >  8:239263/1‑149 (MQ=255)
                                                                                                                                                cATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGGGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGCCTGCCCTGCGAGGAAGCCCCCATGACCGCCGGc   <  6:280425/149‑1 (MQ=255)
                                                                                                                                                 aTGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGCCTGCCCTGCGAGGAAGCCCCCATGACCGCCGGcc  >  6:154190/1‑149 (MQ=255)
                                                                                                                                                 aTGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGCCTGCCCTGCGAGGAAGCCCCCATGACCGCCGGcc  >  5:19496/1‑149 (MQ=255)
                                                                                                                                                 aTGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGCCTGCCCTGCGAGGAAGCCCCCATGACCGCCGGcc  >  8:313446/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCCCACCAGATAGACGCGATCAGCCTGGTAGCCCTTGAGAAACTGCGACAGGTGTTCGGCCTGGCAGTCGCGGGTACCCAGGTGCACGTCGGAGAT‑CCACAGGGTACGCACGCGTTGCTTGCGCGAGGGGCGGGACAGTTCGGCATGGCTCATGGGCAGGCTCCGGCGCAGTTTGCTGGAGACTGGCAGGCACACATGACAGCGCCATGGCAAAACGGTGACGAGTGATGGCGTGCAGGGATTGAGGCACAATGGGCGCCTGCCCTGCGAGGAAGCCCCCATGACCGCCGGCC  >  NC_002947/5123979‑5124271

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