Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,196,441 C→T 9.3% R231R (CGG→CGA)  PP_2804 ← alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,196,4410CT 9.3% 178.2 / 5.6 64R231R (CGG→CGA) PP_2804alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base C (30/28);  new base T (3/3);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

CCGACAGGCACGTCGGCAAGTCATGCCCCATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAGGCCGGCTCGTTCCACTCCAGC  >  NC_002947/3196296‑3196582
                                                                                                                                                 |                                                                                                                                             
ccGACAGGCCCGTCGGCAAGTCATGCCCCATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCAGGAGGTCTTTGTCGCCATGAATGCCCAGCGTGGGGCAGCGAATCCGc                                                                                                                                            >  3:271495/1‑149 (MQ=255)
  gACAGGCACGTCGGCAAGTCATGCCCCATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTg                                                                                                                                          <  6:180497/149‑1 (MQ=255)
  gACAGGCACGTCGGCAAGTCATGCCCCATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTg                                                                                                                                          >  2:201080/1‑149 (MQ=255)
       gCACGTCGGCAAGTCATGCCCCATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGGATTCGCTGcccca                                                                                                                                     <  2:204107/149‑5 (MQ=255)
          cGTCGGCAAGTCATGCCCCATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCt                                                                                                                                  <  1:96410/149‑1 (MQ=255)
               gCAAGTCATGCCCCATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTa                                                                                                                             >  8:189929/1‑149 (MQ=255)
                 aaGTCATGCCCCATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACg                                                                                                                           >  5:99967/1‑149 (MQ=255)
                   gTCATGCCCCATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGt                                                                                                                         >  5:342290/1‑149 (MQ=255)
                    tCATGCCCCATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGGATTCGCGGcccca                                                                                                                                     <  1:17684/136‑5 (MQ=255)
                    tCATGCCCCATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGGATTCACTGcccca                                                                                                                                     >  2:17684/1‑132 (MQ=255)
                    tCATGCCCCATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCt                                                                                                                                           <  6:257604/130‑1 (MQ=255)
                    tCATGCCCCATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCt                                                                                                                                           >  5:257604/1‑130 (MQ=255)
                          cccATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGa                                                                                                                  <  7:220672/149‑1 (MQ=255)
                           ccATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGTCGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAg                                                                                                                 >  4:138011/1‑149 (MQ=255)
                           ccATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAg                                                                                                                 >  7:342906/1‑149 (MQ=255)
                              tCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGtt                                                                                                               >  8:96336/1‑148 (MQ=255)
                                 cGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCt                                                                                                                                           <  2:359095/117‑1 (MQ=255)
                                 cGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCt                                                                                                                                           >  1:359095/1‑117 (MQ=255)
                                    gcagcagcaCCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATc                                                                                                        >  3:48371/1‑149 (MQ=255)
                                      agcagcaCCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATcgc                                                                                                      >  5:250735/1‑149 (MQ=255)
                                       gcagcaCCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATcgcg                                                                                                     >  7:272772/1‑149 (MQ=255)
                                        cagcaCCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTtgat                                                                                                            >  5:80197/1‑141 (MQ=255)
                                        cagcaCCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTtgat                                                                                                            <  6:80197/141‑1 (MQ=255)
                                              cAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCt                                                                                              >  8:100561/1‑149 (MQ=255)
                                              cAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCt                                                                                              >  4:225109/1‑149 (MQ=255)
                                                    ggCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTg                                                                                         <  8:272772/148‑1 (MQ=255)
                                                     gcgcCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTgg                                                                                        <  6:342290/148‑1 (MQ=255)
                                                        ccgcGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACGAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGt                                                                                    <  1:114973/149‑1 (MQ=255)
                                                        ccgcGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGACGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCGACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGt                                                                                    >  1:280238/1‑149 (MQ=255)
                                                          gcgAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGc                                                                                  >  4:259223/1‑149 (MQ=255)
                                                          gcgAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGc                                                                                  <  7:189929/149‑1 (MQ=255)
                                                                 gCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGt                                                                            <  3:225109/148‑1 (MQ=255)
                                                                   cGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGc                                                                         <  7:96336/149‑1 (MQ=255)
                                                                     aCGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTg                                                                       >  3:32589/1‑149 (MQ=255)
                                                                       gCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGc                                                                     <  6:250735/149‑1 (MQ=255)
                                                                         cGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCt                                                                   <  3:259223/149‑1 (MQ=255)
                                                                               aaaCCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGc                                                             <  8:116778/149‑1 (MQ=255)
                                                                                aaCCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGc                                                                                            <  2:337933/117‑1 (MQ=255)
                                                                                aaCCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGc                                                                                            >  1:337933/1‑117 (MQ=255)
                                                                                aaCCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGcc                                                            <  1:201080/149‑1 (MQ=255)
                                                                                aaCCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGc                                                             >  6:135450/1‑148 (MQ=255)
                                                                                 aCCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCa                                                           >  1:236000/1‑149 (MQ=255)
                                                                                     gccgTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTc                                                       >  2:312705/1‑149 (MQ=255)
                                                                                         ttGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTc                                                                  >  3:103543/1‑134 (MQ=255)
                                                                                         ttGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTc                                                                  <  4:103543/134‑1 (MQ=255)
                                                                                            gtggCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGcc                                                 <  4:32589/148‑1 (MQ=255)
                                                                                                          gTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCg                                   >  5:52037/1‑148 (MQ=255)
                                                                                                                    ccATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCaa                        <  6:145040/149‑1 (MQ=255)
                                                                                                                      atgaatgaCCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAgg                      >  4:60459/1‑149 (MQ=255)
                                                                                                                       tgaatgaCAAGCGTGGTGCAGCGAATCCGCTGCAACTCCTGTGTACTGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAGGc                     <  4:48371/149‑1 (MQ=255)
                                                                                                                            gaCCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAGGCCGGCt                >  2:229347/1‑149 (MQ=255)
                                                                                                                              ccAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCg                                                                                                                    <  6:123944/47‑1 (MQ=255)
                                                                                                                              ccAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCg                                                                                                                    >  5:123944/1‑47 (MQ=255)
                                                                                                                              ccAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAGGCCGGCTCg              >  3:111319/1‑149 (MQ=255)
                                                                                                                              ccAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAGGCCGGCTCg              <  1:312705/149‑1 (MQ=255)
                                                                                                                                 gCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAGGCCGGCTCGTTc           <  6:99967/149‑1 (MQ=255)
                                                                                                                                    tGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAGGCCGGCTCGTTCCAc        <  5:135450/149‑1 (MQ=255)
                                                                                                                                    tGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAGGCCGGCGCGTTCCAc        <  8:342906/149‑1 (MQ=255)
                                                                                                                                      ttgCAGCGGATTCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAGGCCGGCTCGTTCCa         <  4:1387/144‑1 (MQ=255)
                                                                                                                                      ttgCAGCGGATTCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAGGCCGGCTCGTTCCa         >  3:1387/3‑146 (MQ=255)
                                                                                                                                      ttgCAGCGGATTCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAGGCCGGCTCGTTCCACTc      >  7:57865/3‑149 (MQ=255)
                                                                                                                                        gCAGCGAATCCGCTGCGACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAGGCCGGCTCGTTCCACTCCa    <  2:280238/149‑1 (MQ=255)
                                                                                                                                         cAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAGGCCGGCTCGTTCCACTCCAg   >  2:304340/1‑149 (MQ=255)
                                                                                                                                          aGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAGGCCGGCTCGTTCCACTCCAGc  <  7:100561/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CCGACAGGCACGTCGGCAAGTCATGCCCCATCCCGGGCAGCAGCACCAGGTTGGCGCCGCGAATCGCCGACGCCGTGGCAAACCCGCCGTTGGTGGCGACCATGAGGTCTTTGTCGCCATGAATGACCAGCGTGGGGCAGCGAATCCGCTGCAACTCCTGGGTACGGTCACCGGAGTTGATGATCGCGCCGATCTGCCTGGCGGTGCCGAGGTTGCTGGCCTCACCGCCACCTCGCGGCCAGGCTTGCATGGCGTAGTCACGCAAGGCCGGCTCGTTCCACTCCAGC  >  NC_002947/3196296‑3196582

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