Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,061,400 C→T 6.6% L105L (CTC→CTT)  PP_1838 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,061,4000CT 6.6% 175.0 / 2.0 60L105L (CTC→CTT) PP_1838membrane protein
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base T (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

CCCCTTATTCGACGGTCGCCGAAGGCTTCCTCAACGGCTTGCGCTTCGACCTGCGGGTGGTGGTGTACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGACCCGAAAACCGTCATGAGCATGCTCTGGTACGGTTTCCCGGTG  >  NC_002947/2061252‑2061544
                                                                                                                                                    |                                                                                                                                                
ccccTTATTCGACGGTCGCCGAAGGCTTCCTCAACGGCTTGCGCTTCGACCTGCGGGTGGTGGTGTACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTc                                                                                                                                                  >  6:50810/1‑149 (MQ=255)
 cccTTATTCGACGGTCGCCGAAGGCTTCCTCAACGGCTTGCGCTTCGACCTGCGGGTGGTGGTGTACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCg                                                                                                                                                 >  4:78833/1‑149 (MQ=255)
            cGGTCGCCGAAGGCTTCCTCAACGGCTTGCGCTTCGACCTGCGGGTGGTGGTGTACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCgtg                                                                                                                                      <  3:78833/149‑1 (MQ=255)
                  ccGAAGGCTTCCTCAACGGCTTGCGCTTCGACCTGCGGGTGGTGGTGTACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATg                                                                                                                                >  4:2339/1‑149 (MQ=255)
                   cGAAGGCTTCCTCAACGGCTTGCGCTTCGACCTGCGGGTGGTGGTGTACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGt                                                                                                                               >  8:255464/1‑149 (MQ=255)
                                      ttGCGCTTCGACCTGCGGGTGGTGGTGTACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGaaa                                                                                                             >  7:89983/1‑148 (MQ=255)
                                      ttGCGCTTCGACCTGCGGGTGGTGGTGTACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAAt                                                                                                            >  4:333097/1‑149 (MQ=255)
                                        gcgcTTCGACCTGCGGGTTGTGGTGTACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATg                                                                                                           <  8:115337/148‑1 (MQ=255)
                                        gcgcTTCGACCTGCGGGTGGTGGTGTACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGaaa                                                                                                             <  8:89983/146‑1 (MQ=255)
                                        gcgcTTCGACCTGCGGGTGGTGGTGTACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATg                                                                                                           <  6:305474/148‑1 (MQ=255)
                                          gcTTCGACCTGCGGGTGGTGGTGTACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGAc                                                                                                        <  3:333097/149‑1 (MQ=255)
                                            ttCGACCTGCGGGTGGTGGTGTACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACt                                                                                                       >  7:331269/1‑148 (MQ=255)
                                                     cGGGTGGTGGTGTACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGa                                                                                             <  1:159809/149‑1 (MQ=255)
                                                                gtACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGAt                                                                                                                                 >  2:327360/1‑102 (MQ=255)
                                                                    aTCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGa                                                                                                                                           <  6:50468/88‑1 (MQ=255)
                                                                    aTCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGGAATGGACTTCTACCGCGAGTTCCACCAGCGCCt                                                                              >  5:271532/1‑149 (MQ=255)
                                                                    aTCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCt                                                                              >  1:129174/1‑149 (MQ=255)
                                                                     tCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCACCt                                                                              >  2:209158/1‑148 (MQ=255)
                                                                       aGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCaa                                                                           >  2:211331/1‑149 (MQ=255)
                                                                        gCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAAc                                                                          <  5:50810/149‑1 (MQ=255)
                                                                             ccccTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCt                                                                     >  6:205410/1‑149 (MQ=255)
                                                                              cccTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTg                                                                    <  3:2339/149‑1 (MQ=255)
                                                                                   ctgctgGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGATATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGgtgt                                                                <  1:253626/148‑1 (MQ=255)
                                                                                    tgctgGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCACCGGCCTGgtgtgc                                                              >  5:41139/1‑147 (MQ=255)
                                                                                    tgctgGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTc                                                              >  7:23167/1‑149 (MQ=255)
                                                                                    tgctgGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTc                                                              >  3:334845/1‑149 (MQ=255)
                                                                                     gctgGCCATGCTCAGCCCCTGGGCCATGGCCCGGCCAGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCCCCTCAACGGCCTGGTGTTc                                                              >  2:255681/1‑148 (MQ=255)
                                                                                          ccATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTa                                                         >  3:36084/1‑148 (MQ=255)
                                                                                           cATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGAc                                                                                                        >  2:254525/1‑100 (MQ=255)
                                                                                           cATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGAc                                                                                                        <  1:249888/100‑1 (MQ=255)
                                                                                           cATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGAc                                                                                                        <  1:254525/100‑1 (MQ=255)
                                                                                           cATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGAc                                                                                                        >  2:249888/1‑100 (MQ=255)
                                                                                                tCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAg                                                  >  3:148159/1‑149 (MQ=255)
                                                                                                tCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAg                                                  >  8:179297/1‑149 (MQ=255)
                                                                                                  aGCCCCTGGTCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGGTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAg                                                  <  8:331269/147‑1 (MQ=255)
                                                                                                  aGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAg                                                  <  1:211331/147‑1 (MQ=255)
                                                                                                      ccTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGGAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGAc                                            <  6:271532/149‑1 (MQ=255)
                                                                                                      ccTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGAc                                            <  2:129174/149‑1 (MQ=255)
                                                                                                      ccTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGAc                                            <  8:23167/149‑1 (MQ=255)
                                                                                                       cTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTTGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGAcc                                           >  5:70351/1‑149 (MQ=255)
                                                                                                       cTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTTGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGAcc                                           >  5:65199/1‑149 (MQ=255)
                                                                                                             cATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTc                                                              >  6:280051/1‑124 (MQ=255)
                                                                                                             cATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGa                                             >  7:82844/1‑141 (MQ=255)
                                                                                                             cATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGa                                             <  8:82844/141‑1 (MQ=255)
                                                                                                             cATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGACCCGaaaa                                     <  7:255464/149‑1 (MQ=255)
                                                                                                               tGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGACCCGAAAAcc                                   <  7:179297/149‑1 (MQ=255)
                                                                                                                  cccGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTTGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGACCCGAAAACCGTc                                <  6:65199/149‑1 (MQ=255)
                                                                                                                  cccGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTTGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGACCCGAAAACCGTc                                <  6:70351/149‑1 (MQ=255)
                                                                                                                   ccGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCa                                                            >  5:268967/1‑120 (MQ=255)
                                                                                                                   ccGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCa                                                            <  6:268967/120‑1 (MQ=255)
                                                                                                                   ccGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGACCCGAAAACCGTCa                               >  2:39672/1‑149 (MQ=255)
                                                                                                                      gcgtggcgtgTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAgg                                                 >  5:198245/1‑128 (MQ=255)
                                                                                                                      gcgtggcgtgTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAgg                                                 <  6:198245/128‑1 (MQ=255)
                                                                                                                      gcgtggcgtgTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGACCCGAAAACCGTCAt                              >  2:190273/1‑147 (MQ=255)
                                                                                                                                   cGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGACCCGAAAACCGTCATGAGCATGCTCTGGTa               >  2:81969/1‑149 (MQ=255)
                                                                                                                                    gTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGACCCGAAAACCGTCATGAGCATGCTCTGGTAc              <  1:190273/149‑1 (MQ=255)
                                                                                                                                       tCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGACCCGAAAACCGTCATGAGCATGCTCTGGTACGGt           >  3:44796/1‑149 (MQ=255)
                                                                                                                                            cTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGACCCGAAAACCGTCATGAGCATGCTCTGGTACGGTTTccc      <  4:44796/149‑1 (MQ=255)
                                                                                                                                               aCCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGACCCGAAAACCGTCATGAGCATGCTCTGGTACGGTTTCCCggt   <  4:148159/149‑1 (MQ=255)
                                                                                                                                                cccTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAATGAAGACCCGAAAACCGTCATGAGCATGCTCTGGTACGGTTTCCCggtg  <  1:255681/149‑1 (MQ=255)
                                                                                                                                                cccTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGACCCGAAAACCGTCATGAGCATGCTCTGGTACGGTTTCCCggtg  <  5:205410/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CCCCTTATTCGACGGTCGCCGAAGGCTTCCTCAACGGCTTGCGCTTCGACCTGCGGGTGGTGGTGTACATCAGCATCCCCCTGCTGCTGGCCATGCTCAGCCCCTGGGCCATGGCCCGGCGTGGCGTGTTCCGTTTCTGGCTGACCCTCGCTTCGAGCGTGGTGATGTTCCTCGGCCTGATGGAAATGGACTTCTACCGCGAGTTCCACCAGCGCCTCAACGGCCTGGTGTTCCAGTACATCAAGGAAGACCCGAAAACCGTCATGAGCATGCTCTGGTACGGTTTCCCGGTG  >  NC_002947/2061252‑2061544

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