Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,483,006 T→A 5.8% Q973L (CAG→CTG)  PP_3091 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,483,0060TA 5.8% 208.8 / 2.8 69Q973L (CAG→CTG) PP_3091membrane protein
Reads supporting (aligned to +/- strand):  ref base T (27/38);  new base A (2/2);  total (29/40)
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.40e-01

GGCTGAACAGCCCCTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACACCCGCGCCACCGGCTCGAAGCGCTTGGCCAGGCTGC  >  NC_002947/3482860‑3483147
                                                                                                                                                  |                                                                                                                                             
ggCTGAACAGCCCCTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGc                                                                                                                                             >  8:160698/1‑149 (MQ=255)
 gCTGAACAGCCCCTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCg                                                                                                                                            >  2:368044/1‑149 (MQ=255)
  cTGAACAGACCCTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGGCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGa                                                                                                                                           <  1:130136/149‑1 (MQ=255)
         gCCCCTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCAGATCCTGCGAGCGCTgg                                                                                                                                    <  5:371916/149‑1 (MQ=255)
           cccTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGAt                                                                                                                                  <  3:125704/149‑1 (MQ=255)
                cAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGAtgttgt                                                                                                                             <  2:176410/149‑1 (MQ=255)
                cAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCAGCGCTCGCAGGAtgtgg                                                                                                                              <  6:320658/148‑3 (MQ=255)
                    gagcgCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCAt                                                                                                                         <  5:352751/149‑1 (MQ=255)
                    gagcgCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCAt                                                                                                                         <  7:160698/149‑1 (MQ=255)
                      gcgcGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATgc                                                                                                                       >  5:156595/1‑149 (MQ=255)
                             aGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCtt                                                                                                                <  4:102549/149‑1 (MQ=255)
                               gCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCtt                                                                                                                >  5:175325/1‑147 (MQ=255)
                               gCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGa                                                                                                              <  4:248475/149‑1 (MQ=255)
                               gCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGa                                                                                                              <  4:197574/149‑1 (MQ=255)
                                ccATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAg                                                                                                             <  3:363305/149‑1 (MQ=255)
                                 cATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCAGCgctcgca                                                                                                                                      >  4:178665/1‑117 (MQ=255)
                                 cATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCAGCgctcgca                                                                                                                                      <  3:178665/123‑7 (MQ=255)
                                     ctggctggTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTgc                                                                                                        <  2:171396/149‑1 (MQ=255)
                                     ctggctggTGTCGACGCTGCTTTCCACGTAGGTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTgc                                                                                                        <  7:222719/149‑1 (MQ=255)
                                        gctggTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTgcgc                                                                                                      >  1:325426/1‑148 (MQ=255)
                                               tCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAggcggc                                                                                              >  3:157690/1‑149 (MQ=255)
                                                 gACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAggcggc                                                                                              >  1:38450/1‑147 (MQ=255)
                                                 gACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAggcggc                                                                                              >  3:173721/1‑147 (MQ=255)
                                                 gACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAg                                                                                                   >  4:350593/1‑142 (MQ=255)
                                                 gACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAg                                                                                                   <  3:350593/142‑1 (MQ=255)
                                                   cGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAgg                                                                                          <  6:156595/149‑1 (MQ=255)
                                                   cGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAgg                                                                                          >  1:200887/1‑149 (MQ=255)
                                                     ctgctTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTa                                                                                        <  8:184472/149‑1 (MQ=255)
                                                        ctTTCCCCGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAg                                                                                       <  2:308324/147‑1 (MQ=255)
                                                        ctTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAg                                                                                       <  1:36634/147‑1 (MQ=255)
                                                             cACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCa                                                                                <  3:126441/149‑1 (MQ=255)
                                                             cACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCa                                                                                <  1:368044/149‑1 (MQ=255)
                                                             cACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCa                                                                                <  4:173721/149‑1 (MQ=255)
                                                               cGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCat                                                                               >  6:323674/1‑148 (MQ=255)
                                                                   gTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCatgatg                                                                           <  4:189115/148‑1 (MQ=255)
                                                                   gTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATgg                                                                          >  2:35901/1‑149 (MQ=255)
                                                                        gcACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTg                                                                     >  6:210259/1‑149 (MQ=255)
                                                                         cACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGt                                                                    <  5:323674/149‑1 (MQ=255)
                                                                              tggtgATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGgcg                                                               >  6:105784/1‑149 (MQ=255)
                                                                                      tcgctcgGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGc                                                       >  3:123748/1‑149 (MQ=255)
                                                                                      tcgctcgGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGc                                                       >  5:151662/1‑149 (MQ=255)
                                                                                        gctcgGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCgg                                                        <  7:268510/146‑1 (MQ=255)
                                                                                         ctcgGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCgtt                                                    >  3:46515/1‑149 (MQ=255)
                                                                                            ggCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCgttg                                                   <  2:38450/147‑1 (MQ=255)
                                                                                                 gCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTc                                              <  2:200887/147‑1 (MQ=255)
                                                                                                 gCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCg                                             <  1:35901/148‑1 (MQ=255)
                                                                                                   cTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCg                                          <  2:325426/149‑1 (MQ=255)
                                                                                                       cAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAaca                                      <  8:279943/149‑1 (MQ=255)
                                                                                                            tctcGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACACCcgc                                 <  6:296026/149‑1 (MQ=255)
                                                                                                                 cTGATCAAGTGTTTGCTGCTGTTTCCCAGATCCTTCGAGCTCTGGTTGTTGTTCATGAGCATCTTGATTTTGCGCAAGGCGTGCATGTAGCGATGCATGATGGTGCTGTCGGCGCCTTCGGGGTTGTTTTCGGGGAACTCCCGCGCCa                             <  6:151662/148‑1 (MQ=255)
                                                                                                                         cTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATc                                                                                  >  1:159921/1‑87 (MQ=255)
                                                                                                                         cTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATc                                                                                  <  2:159921/87‑1 (MQ=255)
                                                                                                                            tttGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACACCCGCGCCACCgg                         <  5:370080/141‑1 (MQ=255)
                                                                                                                            tttGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACACCCGCGCCACCgg                         >  6:370080/1‑141 (MQ=255)
                                                                                                                            tttGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTAGGCGAACACCCGCGCCCCCGGCTCGAAgc                 >  2:260931/1‑149 (MQ=255)
                                                                                                                                ctgctCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTGGTTTTCGGCGAACACCCGCGCCACCGGCTCGAAGCGCtt             <  4:123748/149‑1 (MQ=255)
                                                                                                                                   ctctTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACACCCGCGCCACCGGCTCGAAGCGCTTg            <  4:157690/147‑1 (MQ=255)
                                                                                                                                   ctctTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACACCCGCGCCACCGGCTCGAAGCGCTTGGc          <  4:46515/149‑1 (MQ=255)
                                                                                                                                    tctTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGGTGTTTTAGGCGAAAACCCGCGCCACCGGCGCGCAGCGCTTGGcc         >  7:93022/1‑149 (MQ=255)
                                                                                                                                        gCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTCGTTTTCGGCGAACACCCGCGCCACCGGCTCGAAGCGCTTGGCCAg       <  1:260931/147‑1 (MQ=255)
                                                                                                                                         cccACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCATTCGGCGTTGTTTTCGGCGAACACCCGCGCCACCGGCTCGAAGCGCTTGGCCAGGCt    >  3:147819/1‑149 (MQ=255)
                                                                                                                                         acaacaTCCAGCGCTCGCAGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCGTAGGCGGTGTTTTAGGCGAACACCCCCGCCAACCGCGCGAAGCGCTTGGCCAGGCt    >  5:334116/4‑149 (MQ=255)
                                                                                                                                          ccACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACACCCGCGCCACCGGCTCGAAGCGCTTGGCCAGGCTg   >  6:374585/1‑149 (MQ=255)
                                                                                                                                          ccACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACACCCGCGCCACCGGCTCGAAGCGCTTGGCCAGGCTg   >  3:137265/1‑149 (MQ=255)
                                                                                                                                           cacaTCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCAGGATGGTGCTGTCGGCGCCTTCGGCGGTGGGTTCGGCGAACACCCGCGCCCCCGGCGCGAAGGGCTTGGCCAGGCTGc  >  2:258663/1‑149 (MQ=255)
                                                                                                                                              aTCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCatgatg                                                                           >  3:12585/1‑73 (MQ=255)
                                                                                                                                              aTCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCatgatg                                                                           <  4:12585/73‑1 (MQ=255)
                                                                                                                                              aTCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACACCCGCGCCACCGGCTCGAAGCGCtt             <  1:324508/135‑1 (MQ=255)
                                                                                                                                              aTCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCGAGGGCGGTGTTTTAGGCGCAAACCCCCGCCACCCCCTCGACgcgcgg             >  2:324508/1‑133 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GGCTGAACAGCCCCTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACGTAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACACCCGCGCCACCGGCTCGAAGCGCTTGGCCAGGCTGC  >  NC_002947/3482860‑3483147

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