Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,907,691 T→A 13.0% Q335L (CAG→CTG)  hasE ← secretion metalloprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,907,6910TA13.0% 154.3 / 11.7 60Q335L (CAG→CTG) hasEsecretion metalloprotease
Reads supporting (aligned to +/- strand):  ref base T (23/29);  new base A (3/5);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGTCGAGCAGCCGGTCGGCGGAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCATAGTCCAGTGCGCGC  >  NC_002947/2907543‑2907837
                                                                                                                                                    |                                                                                                                                                  
ttcGTCGAGCAGCCGGTCGGCGGAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCt                                                                                                                                                    >  3:295692/1‑149 (MQ=255)
ttcGTCGAGCAGCCGGTCGGCGGAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCt                                                                                                                                                    >  2:93355/1‑149 (MQ=255)
        gcagcCGGTCGGCGGAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGc                                                                                                                                            <  2:239333/149‑1 (MQ=255)
          agcCGGTCGGCGTAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTg                                                                                                                                          <  1:353264/149‑1 (MQ=255)
             cGGTCGGCGGAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTc                                                                                                                                       <  4:295692/149‑1 (MQ=255)
             cGGTCGGCGGAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCAGCCGCTGCAGGTc                                                                                                                                       <  4:113303/149‑1 (MQ=255)
               gTCGGCGGAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGa                                                                                                                                     <  6:339924/149‑1 (MQ=255)
                  ggcggAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGc                                                                                                                                   <  8:284644/148‑1 (MQ=255)
                  ggcggAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCAGCCGCTGCAGGTCGatg                                                                                                                                   <  7:128457/148‑3 (MQ=255)
                   gcggAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTg                                                                                                                                          >  6:273877/1‑140 (MQ=255)
                   gcggAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTg                                                                                                                                          <  5:273877/140‑1 (MQ=255)
                    cggAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGt                                                                                                                                <  1:161873/149‑1 (MQ=255)
                     ggAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCAGCCGCTGCAGGTCGat                                                                                                                                    >  8:145093/1‑143 (MQ=255)
                     ggAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCAGCCGCTGCAGGTCGat                                                                                                                                    <  7:145093/144‑2 (MQ=255)
                        aaTGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCa                                                                                                                            >  8:228118/1‑149 (MQ=255)
                         aTGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCAc                                                                                                                           <  7:228118/149‑1 (MQ=255)
                          tGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCAGCCGCTGCAGGTCGat                                                                                                                                    <  5:138762/139‑2 (MQ=255)
                          tGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCAGCCGCTGCAGGTCGat                                                                                                                                    >  6:138762/1‑138 (MQ=255)
                             tttttACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTcc                                                                                                                       <  1:176657/149‑1 (MQ=255)
                                 tACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATg                                                                                                                    <  4:102103/148‑1 (MQ=255)
                                  aCCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGa                                                                                                                   >  5:100007/1‑148 (MQ=255)
                                   cccGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATc                                                                                                                 <  1:173738/149‑1 (MQ=255)
                                       gccggGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTgg                                                                                                             >  7:124705/1‑149 (MQ=255)
                                       gccggGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGGCGAGCGGTACCACTTCCATGATCTTgg                                                                                                             <  5:207955/149‑1 (MQ=255)
                                          gggAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCACATCGACCAGCCGCTGCAGGTCGat                                                                                                                                    <  6:12052/123‑2 (MQ=255)
                                          gggAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCACATCGACCAGCCGCTGCAGGTCGat                                                                                                                                    >  5:12052/1‑122 (MQ=255)
                                                       gTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGgatgat                                                                                             >  4:89662/1‑149 (MQ=255)
                                                              cATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATgccgcc                                                                                       >  1:225465/1‑148 (MQ=255)
                                                               aTGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGa                                                                                     >  1:352740/1‑149 (MQ=255)
                                                                  gTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGAcc                                                                                   <  8:124705/148‑1 (MQ=255)
                                                                   ttGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATc                                                                                                                 >  6:2620/1‑117 (MQ=255)
                                                                   ttGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATc                                                                                                                 <  5:2620/117‑1 (MQ=255)
                                                                      aaGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTgg                                                                               <  2:225465/148‑1 (MQ=255)
                                                                      aaGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTgg                                                                               <  6:810/148‑1 (MQ=255)
                                                                      aaGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGc                                                                              <  6:173537/149‑1 (MQ=255)
                                                                         gCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGAt                                                                           >  2:191794/1‑149 (MQ=255)
                                                                                    gATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCt                                                                <  2:182443/149‑1 (MQ=255)
                                                                                          cACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCACGGGCTGGATGATGCCGCCGACCGTGGCGATGCACAGGCCCTGGACCa                                                          >  3:329988/1‑149 (MQ=255)
                                                                                            cGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATg                                                        <  1:93355/149‑1 (MQ=255)
                                                                                             ggCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGc                                                       >  4:363920/1‑149 (MQ=255)
                                                                                                        ggACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGc                                                                         >  5:136684/1‑120 (MQ=255)
                                                                                                        ggACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGc                                                                         <  6:136684/120‑1 (MQ=255)
                                                                                                         gACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTgg                                                                                                  >  8:106198/1‑94 (MQ=255)
                                                                                                         gACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTgg                                                                                                  <  7:106198/94‑1 (MQ=255)
                                                                                                                ttGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTACATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGTCGTCGATCGGCGAGCGGAt                                    >  1:170703/1‑149 (MQ=255)
                                                                                                                      aTCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCgt                              <  3:123386/149‑1 (MQ=255)
                                                                                                                       tCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCgtg                             >  7:291970/1‑149 (MQ=255)
                                                                                                                        cGCCTGCACCGGGTTCATCGATTAGACGTTTAGCGGCTGGTCGTTCTGTACCACTTCTATGTTCTTGGGCGCGGGGTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCGTGGACCATGCCGTCGATCGGCGAGCGGATCACCgtg                             <  4:329988/148‑1 (MQ=255)
                                                                                                                          ccTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGTCGTCGATCGGCGAGCGGATCACCGTGTgg                          <  2:170703/149‑1 (MQ=255)
                                                                                                                          ccTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTgg                          >  5:135948/1‑149 (MQ=255)
                                                                                                                               aCCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTgg                          <  4:189861/144‑1 (MQ=255)
                                                                                                                               aCCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTgg                          >  3:189861/1‑144 (MQ=255)
                                                                                                                                ccGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACt                    <  3:89662/149‑1 (MQ=255)
                                                                                                                                      tcatcGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCAATCGGCGAGCGGATCACCGTGTGGGTCACt                    <  6:314804/143‑1 (MQ=255)
                                                                                                                                      tcatcGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCAATCGGCGAGCGGATCACCGTGTGGGTCACt                    >  5:314804/1‑143 (MQ=255)
                                                                                                                                           gcatcgACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCATAGTCCAg         >  8:12856/1‑149 (MQ=255)
                                                                                                                                           gcatcgACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCATAGTCCAg         >  8:18421/1‑149 (MQ=255)
                                                                                                                                            catcgACCTGCAGCGGCTGGTCGAGCGGTACCACTTCGATGATCTTGGCCCCGGGCTGGATGATGCCGCTGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGGGTGGGTCACTTCATAGTCCAGt        <  6:100007/149‑1 (MQ=255)
                                                                                                                                             atcgACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCATAGTCCAGTg       >  8:194191/1‑149 (MQ=255)
                                                                                                                                              tcgACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCATAGTCCAGTgc      <  7:194191/149‑1 (MQ=255)
                                                                                                                                                  ccTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCATAGTCCAGTgcgcg   >  8:254162/1‑148 (MQ=255)
                                                                                                                                                   cTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCATAGTCCAGTgcgcgc  >  8:335283/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTCGTCGAGCAGCCGGTCGGCGGAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGAGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCATAGTCCAGTGCGCGC  >  NC_002947/2907543‑2907837

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: