Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,267,264 G→C 6.2% P242A (CCA→GCA)  quiP ← acyl‑homoserine lactone acylase QuiP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,267,2640GC 6.2% 189.0 / 2.3 60P242A (CCA→GCA) quiPacyl‑homoserine lactone acylase QuiP
Reads supporting (aligned to +/- strand):  ref base G (31/24);  new base C (2/2);  total (33/27)
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

GGGCTCAAGGCCCAGGCGGCGCGGCTGTCGCTGGCCAGCAGGCTCTTGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGCTACCGGCTTTCTGCGCCAGGGTCAGGGCGCTGAGTTCCTCT  >  NC_002947/1267116‑1267412
                                                                                                                                                    |                                                                                                                                                    
gggCTCAAGGCCCAGGCGGCGCGGCTGTCGCTGGCCAGCAGGCTCTTGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTgg                                                                                                                                                      >  6:46375/1‑149 (MQ=255)
    tCAAGGCCCAGGCGGCGCGGCTGTCGCTGGCCAGCAGGCTCTTGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGc                                                                                                                                                  >  4:277589/1‑149 (MQ=255)
      aaGGCCCAGGCGGCGCGGCTGTCGCTGGCCAGCAGGCTCTTGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTgg                                                                                                                                                      <  7:147090/143‑1 (MQ=255)
      aaGGCCCAGGCGGCGCGGCTGTCGCTGGCCAGAAGGCTCTTGCCGCTGCGGCTGCGCTGCGGACCCAGGCCCAGGATGCCCGACCCCGTGCTACCGAGCAGGTACAGGTCGGCCAGCTTCAGGCTGGCAGCCGCCAGAGCTCGCAGctc                                                                                                                                                >  8:147090/1‑147 (MQ=255)
      aaGGCCCAGGCGGCGCGGCTGTCGCTAGCCAGCAGGCTCTTGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGccc                                                                                                                                                >  7:270386/1‑149 (MQ=255)
             aGGCGGCGCGGCTGTCGCTGGCCAGCAGGCTCTTGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAAc                                                                                                                                         <  8:239491/149‑1 (MQ=255)
                    gcggcTGTCGCTGGCCAGCAGGCTCTTGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAActggctgg                                                                                                                                  >  1:280133/1‑149 (MQ=255)
                       gcTGTCGCTGGCCAGCAGGCTCTTGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACtggctggc                                                                                                                                 >  7:214314/1‑147 (MQ=255)
                        ctgaaccTGGCCAGCAGGGTCTTGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGTCAGCCCCGGCAACTGGCTGGCCAg                                                                                                                              <  3:244739/143‑1 (MQ=255)
                              cTGGCCAGCAGGCTCTTGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAg                                                                                                                        <  2:220822/149‑1 (MQ=255)
                                   cagcagGCTCTTGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCt                                                                                                                   <  4:259397/149‑1 (MQ=255)
                                   cagcagGCTCTTGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCt                                                                                                                   >  5:149556/1‑149 (MQ=255)
                                    agcagGCTCTTGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCtt                                                                                                                  <  5:46375/149‑1 (MQ=255)
                                          ctctTGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCtt                                                                                                            <  8:214314/149‑1 (MQ=255)
                                           tctTGCCGCTGCGGCGGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTg                                                                                                           <  2:176605/149‑1 (MQ=255)
                                             ttGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTc                                                                                                         <  8:270386/149‑1 (MQ=255)
                                               gCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGt                                                                                                          <  6:279216/146‑1 (MQ=255)
                                                 cGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTc                                                                                                         >  6:291821/1‑145 (MQ=255)
                                                 cGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTc                                                                                                         <  5:291821/145‑1 (MQ=255)
                                                       ggctgcgCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGACAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCAtt                                                                                               >  1:193396/1‑149 (MQ=255)
                                                           gcgctgcgGACCCCGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCGTGAGCTTGTCCACTTCCATTTCgg                                                                                           >  5:302474/1‑149 (MQ=255)
                                                               tgcgGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCaa                                                                                       >  5:372190/1‑149 (MQ=255)
                                                                gcgGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAAc                                                                                      >  2:366471/1‑149 (MQ=255)
                                                                     cccAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGGCCACTTCCATTTCGGCCAACGGCTc                                                                                 >  6:262335/1‑149 (MQ=255)
                                                                      ccAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCg                                                                                >  8:29214/1‑149 (MQ=255)
                                                                            ccAGGTTGGCCGACCCCGGGCTACCGCGCAGGTCGAGGTCGGCCAGCTTCTGGCGGGCAGCCGTCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCggg                                                                          <  4:206934/149‑1 (MQ=255)
                                                                             cAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCgggg                                                                         >  5:373787/1‑149 (MQ=255)
                                                                                gTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCa                                                                                                                                                    >  8:261005/1‑71 (MQ=255)
                                                                                gTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCa                                                                                                                                                    <  7:261005/71‑1 (MQ=255)
                                                                                gTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGCCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGccc                                                                      <  2:280133/149‑1 (MQ=255)
                                                                                   ggCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCg                                                                    <  6:374703/148‑1 (MQ=255)
                                                                                    gCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCgg                                                                   <  6:149556/148‑1 (MQ=255)
                                                                                      cGACCCCGGGCTACCGAGCAGGTCGCGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGgca                                                                <  7:267923/149‑1 (MQ=255)
                                                                                      cGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGgca                                                                <  6:71420/149‑1 (MQ=255)
                                                                                          cccGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGgca                                                                <  8:153791/145‑1 (MQ=255)
                                                                                          cccGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGgca                                                                >  7:153791/1‑145 (MQ=255)
                                                                                            cGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAAc                                                           >  4:51334/1‑148 (MQ=255)
                                                                                             gggCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCa                                                         <  6:373787/149‑1 (MQ=255)
                                                                                               gCTACCTAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCAAGGTTCAGGCCCTTGAGCTTGTCAAATTCCATTTCGGACAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAgg                                                       >  2:200735/1‑149 (MQ=255)
                                                                                                 tACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGc                                                      >  3:360853/1‑148 (MQ=255)
                                                                                                 tACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGc                                                      >  2:308715/1‑148 (MQ=255)
                                                                                                    cGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCTCTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGc                                                  >  1:302248/1‑149 (MQ=255)
                                                                                                    cGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGc                                                  >  7:139539/1‑149 (MQ=255)
                                                                                                     gAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCt                                                 >  5:166393/1‑149 (MQ=255)
                                                                                                        cAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGt                                              <  6:372190/149‑1 (MQ=255)
                                                                                                            tcgaggtcgGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGc                                           >  4:97082/1‑148 (MQ=255)
                                                                                                            tcgaggtcgGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGc                                           >  4:102185/1‑148 (MQ=255)
                                                                                                               aggtcgGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTc                                                                                 >  6:367830/1‑107 (MQ=255)
                                                                                                               aggtcgGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTc                                                                                 <  5:367830/107‑1 (MQ=255)
                                                                                                               aggtcgGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGCTAcc                                       <  2:193396/149‑1 (MQ=255)
                                                                                                                      ccAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGCGACCGGCTTTc                                >  2:53836/1‑149 (MQ=255)
                                                                                                                       cAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGCTACCGGCTTTCt                               <  1:179642/149‑1 (MQ=255)
                                                                                                                          cTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGCTACCGGCTTTCTgcg                            >  4:140321/1‑149 (MQ=255)
                                                                                                                           ttCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGCTACCGGCTTTCTgcgc                           >  2:300459/1‑149 (MQ=255)
                                                                                                                                gctggcAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCACGGGCAGCAACCAGGCCAGCTTGTCGCTACCGGCTTTCTGCGCCAggg                      >  5:199107/1‑149 (MQ=255)
                                                                                                                                  tggcAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGCTACCGGCTTTATGCGCCAGGGt                     >  5:246890/1‑148 (MQ=255)
                                                                                                                                       gttgccGGGGCTGCCAGCCCTGGCAACTGGCTTGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGgca                                                                <  6:171274/97‑1 (MQ=255)
                                                                                                                                       gttgccGGGGCTGCCAGCCCTGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGCTACCGGCTTTCTGCGCCAGGGTCAGGgc               >  4:354422/4‑149 (MQ=255)
                                                                                                                                       gttgccGGGGCTGCCAGCCCTGGCAACTGGCTGGCCAGGGTCAGGACCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCAGGgca                                                                >  5:171274/4‑100 (MQ=255)
                                                                                                                                          gccAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGCTACCGGCTTTCTGCGCCAGGGTCAGGGCGCt            <  3:51334/149‑1 (MQ=255)
                                                                                                                                            aaGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTGGTCCACTTCCACTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGCTACCGGCTTTCTGCGCCAGGGTCAGGGCGCTGa          <  7:194844/148‑1 (MQ=255)
                                                                                                                                                 cTGCCAGCCCTGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGCTACCGGCTTTCTGCGCCAGGGTCAGGGCGCTGAGTTcc     <  3:354422/149‑1 (MQ=255)
                                                                                                                                                   ggCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGCTACCGGCTTTCTGCGCCAGGGTCAGGGCGCTGAGTTCctc   >  2:318823/1‑149 (MQ=255)
                                                                                                                                                    gCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCAGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGCTACCGGCTTTCTGCGCCAGGGTCAGGGCGCTGAGTTCctct  >  7:337986/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGGCTCAAGGCCCAGGCGGCGCGGCTGTCGCTGGCCAGCAGGCTCTTGCCGCTGCGGCTGCGCTGCGGACCCAGGGCCAGGTTGGCCGACCCCGGGCTACCGAGCAGGTCGAGGTCGGCCAGCTTCTGGCTGGCAGCCGCCAGGGCTGGCAGCCCCGGCAACTGGCTGGCCAGGTTCAGGCCCTTGAGCTTGTCCACTTCCATTTCGGCCAACGGCTCGTCCGGGGCCCCGGGCAGCAACCAGGCCAGCTTGTCGCTACCGGCTTTCTGCGCCAGGGTCAGGGCGCTGAGTTCCTCT  >  NC_002947/1267116‑1267412

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