Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,996,900 A→C 5.3% M88R (ATG→AGG)  PP_3525 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,996,9000AC 5.3% 237.9 / 4.5 75M88R (ATG→AGG) PP_3525hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/39);  new base C (2/2);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

CAGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTAGGCCTG  >  NC_002947/3996752‑3997048
                                                                                                                                                    |                                                                                                                                                    
cAGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCa                                                                                                                                                      >  7:299119/1‑149 (MQ=255)
cAGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCa                                                                                                                                                      >  2:118396/1‑149 (MQ=255)
  gCACGCCCTGCTGCAGTTGCAACCTCTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATg                                                                                                                                                    >  7:229612/1‑149 (MQ=255)
        ccTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTg                                                                                                                                               >  1:149617/1‑148 (MQ=255)
        ccTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTg                                                                                                                                               >  5:334071/1‑148 (MQ=255)
         ctgctgCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCATGCCTGCGCTGat                                                                                                                                             <  5:160734/149‑2 (MQ=255)
            ctgcAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCATGCCTGCGCTGatc                                                                                                                                            <  4:236301/147‑3 (MQ=255)
            ctgcAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGt                                                                                                                                          <  7:14410/149‑1 (MQ=255)
              gcAGTTGCAACCATTGCACCTGGCCCTCCTGCCCGGGGTACCGACTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCg                                                                                                                                        >  8:43639/1‑149 (MQ=255)
                  ttGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCa                                                                                                                                    >  7:23656/1‑149 (MQ=255)
                  ttGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCa                                                                                                                                    <  6:240828/149‑1 (MQ=255)
                    gCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCAc                                                                                                                                   <  7:101338/148‑1 (MQ=255)
                    gCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCAc                                                                                                                                   <  4:101365/148‑1 (MQ=255)
                     cAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCAc                                                                                                                                   >  1:167039/1‑147 (MQ=255)
                     cAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCAc                                                                                                                                   <  2:167039/147‑1 (MQ=255)
                        ccACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATa                                                                                                                              <  5:245023/149‑1 (MQ=255)
                         cACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCGGGCATGCGCTGCATGTCGGCCACGAATa                                                                                                                              >  7:146425/1‑148 (MQ=255)
                             gCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAgcg                                                                                                                         >  6:329705/1‑149 (MQ=255)
                              cACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAgcgg                                                                                                                        >  8:160499/1‑149 (MQ=255)
                              cACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAgc                                                                                                                          <  8:229612/147‑1 (MQ=255)
                               aCCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAgcggc                                                                                                                       >  1:392090/1‑149 (MQ=255)
                                 cTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCgg                                                                                                                                       >  6:239719/1‑131 (MQ=255)
                                 cTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCgg                                                                                                                                       <  5:239719/131‑1 (MQ=255)
                                    gCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCg                                                                                                                                        >  8:84782/1‑127 (MQ=255)
                                    gCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCg                                                                                                                                        <  7:84782/127‑1 (MQ=255)
                                      cctcctGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCt                                                                                                                >  2:222660/1‑149 (MQ=255)
                                             cccGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAg                                                                                                           <  2:392090/147‑1 (MQ=255)
                                              ccGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTa                                                                                                        <  7:154469/149‑1 (MQ=255)
                                                 gggTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGcc                                                                                                     <  6:367538/149‑1 (MQ=255)
                                                      ccGGCTGGCCACTGGCGCATCACGGTCGAGCAGCACGAGGTAGCTGTAATCGCGCATGGCCGGGACCGCGATCAGTTTGCTGATGAGCGCAGGCATGCGCTGAATGTCGGCCACGAATAGAGCGGCGCGTGCTTCCAGGTAGCCCTTg                                                                                                 <  7:120492/148‑1 (MQ=255)
                                                       cGGCTGGCCACTCGCCCATCACTGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTggg                                                                                               >  3:324410/1‑149 (MQ=255)
                                                       cGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTggg                                                                                               >  7:23520/1‑149 (MQ=255)
                                                            ggCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGgctgct                                                                                          >  3:371744/1‑149 (MQ=255)
                                                              ccACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTggg                                                                                               <  5:296552/142‑1 (MQ=255)
                                                              ccACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTggg                                                                                               >  6:296552/1‑142 (MQ=255)
                                                              ccACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCg                                                                                        <  3:36182/149‑1 (MQ=255)
                                                               cACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCAc                                                                                                                                   <  2:243575/105‑1 (MQ=255)
                                                                aCTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGc                                                                                      >  5:299650/1‑149 (MQ=255)
                                                                aCTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGc                                                                                      >  3:209037/1‑149 (MQ=255)
                                                                    gCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGt                                                                                  <  4:371744/149‑1 (MQ=255)
                                                                           aCGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCt                                                                           <  8:23656/149‑1 (MQ=255)
                                                                             ggTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCtt                                                                          <  7:243368/148‑1 (MQ=255)
                                                                             ggTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCtt                                                                          >  8:243368/1‑148 (MQ=255)
                                                                             ggTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTg                                                                         <  7:43639/149‑1 (MQ=255)
                                                                              gTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGa                                                                        <  2:157250/149‑1 (MQ=255)
                                                                                 gAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAgg                                                                                                          <  6:234436/112‑1 (MQ=255)
                                                                                 gAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAgg                                                                                                          >  1:227831/1‑112 (MQ=255)
                                                                                 gAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAgg                                                                                                          <  2:227831/112‑1 (MQ=255)
                                                                                 gAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAgg                                                                                                          >  5:234436/1‑112 (MQ=255)
                                                                                      gcaCCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGTACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCttt                                                                <  6:339495/149‑2 (MQ=255)
                                                                                      gcaCCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTg                                                                <  5:329705/149‑1 (MQ=255)
                                                                                          ccGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCAc                                                            >  5:158373/1‑149 (MQ=255)
                                                                                           cGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTgg                                                               <  3:373466/145‑1 (MQ=255)
                                                                                           cGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTgg                                                               >  4:373466/1‑145 (MQ=255)
                                                                                             gTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGt                                                         <  2:149617/149‑1 (MQ=255)
                                                                                             gTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGt                                                         <  3:268700/149‑1 (MQ=255)
                                                                                                gCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCa                                                      >  8:20821/1‑149 (MQ=255)
                                                                                                  tGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGt                                                                                  >  7:36307/1‑119 (MQ=255)
                                                                                                  tGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGt                                                                                  <  8:36307/119‑1 (MQ=255)
                                                                                                      aTCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATcgc                                                >  5:363019/1‑149 (MQ=255)
                                                                                                         gcgcATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATcgcgc                                              <  1:118396/148‑1 (MQ=255)
                                                                                                         gcgcATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATcgcgc                                              <  6:334071/148‑1 (MQ=255)
                                                                                                             aTGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGccac                                         <  7:296801/149‑1 (MQ=255)
                                                                                                                     gACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTaaaa                                 <  7:160499/149‑1 (MQ=255)
                                                                                                                          cgAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGccacca                                       <  6:359070/138‑1 (MQ=255)
                                                                                                                          cgAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGccacca                                       >  5:359070/1‑138 (MQ=255)
                                                                                                                           gAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTg                             <  8:23520/147‑1 (MQ=255)
                                                                                                                           gAACAGTTTGCTGATCAGAGCAGGCATGCGCTGAATGTCGCCCACGAATACAGCGGAGCGTGCTTCCACGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGCCCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCg                           >  2:46526/1‑149 (MQ=255)
                                                                                                                            aaCAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTGCATATCGCGCGCCACCAGTAAAATCTGc                            <  1:228287/147‑1 (MQ=255)
                                                                                                                            aaCAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGc                            >  2:228287/1‑147 (MQ=255)
                                                                                                                                gTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCg                      <  6:158373/149‑1 (MQ=255)
                                                                                                                                 tttGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCgg                     >  1:220498/1‑149 (MQ=255)
                                                                                                                                     acattCAGCGCATGCCTGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCAc                  >  3:39691/5‑148 (MQ=255)
                                                                                                                                         tCAGCGCATGCCTGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAgg             >  3:153981/1‑149 (MQ=255)
                                                                                                                                             cgcAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACGAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGt         <  6:363019/149‑1 (MQ=255)
                                                                                                                                             cgcAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGt         <  1:222660/149‑1 (MQ=255)
                                                                                                                                                  gCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTAGGcc    >  6:123832/1‑149 (MQ=255)
                                                                                                                                                    aTGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTAGGCCTg  >  6:104031/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTAGGCCTG  >  NC_002947/3996752‑3997048

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