Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,345,552 T→G 6.2% L317R (CTG→CGG)  acsA‑II → acetyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,345,5520TG 6.2% 189.2 / 2.6 64L317R (CTG→CGG) acsA‑IIacetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base T (25/35);  new base G (2/2);  total (27/37)
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.03e-01

AAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATGCGTGAGGGCTCCGC  >  NC_002947/5345405‑5345699
                                                                                                                                                   |                                                                                                                                                       
aaGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTg                                                                                                                                                        <  6:337459/149‑1 (MQ=255)
     gTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCaa                                                                                                                                                   <  5:103549/149‑1 (MQ=255)
        ctgcATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATg                                                                                                                                                 <  7:304720/148‑1 (MQ=255)
        ctgcATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCAGCGGGCCataga                                                                                                                                                <  8:206121/149‑5 (MQ=255)
         tgcATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGgc                                                                                                                                               >  3:275642/1‑149 (MQ=255)
                  ccGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCg                                                                                                                                     <  1:352868/149‑1 (MQ=255)
                   cggcggCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGc                                                                                                                                    <  3:72593/149‑1 (MQ=255)
                   cggcggCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGc                                                                                                                                    >  2:250569/1‑149 (MQ=255)
                    ggcggcTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCt                                                                                                                                   <  4:275642/149‑1 (MQ=255)
                    ggcggcTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCt                                                                                                                                   <  2:42525/149‑1 (MQ=255)
                    ggcggcTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCt                                                                                                                                   <  2:254506/149‑1 (MQ=255)
                               gctgcAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAg                                                                                                                        >  5:354500/1‑149 (MQ=255)
                                ctgcAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGaa                                                                                                                         <  3:33914/147‑1 (MQ=255)
                                ctgcAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAg                                                                                                                        <  4:370357/148‑1 (MQ=255)
                                ctgcAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCGAGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAg                                                                                                                        <  6:358348/148‑1 (MQ=255)
                                           gATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAAc                                                                                                             >  3:99812/1‑148 (MQ=255)
                                                ccTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTAcccc                                                                                                       <  2:248180/149‑1 (MQ=255)
                                                 cTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTAcccc                                                                                                       >  5:214361/1‑148 (MQ=255)
                                                  ttCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCg                                                                                                      >  5:331820/1‑148 (MQ=255)
                                                  ttCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGa                                                                                                     >  6:13524/1‑149 (MQ=255)
                                                  ttCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATATTCGAAGGCGTGCCCAACTACCCCGa                                                                                                     >  5:303031/1‑149 (MQ=255)
                                                  ttCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCCGATGTTCGAAGGCGTGCCCAACTACCCCGa                                                                                                     >  5:69400/1‑149 (MQ=255)
                                                    gaaGGTGGTGTTCGTTTATCGTGACGGCGAGGTGTTCTGGTGCAGCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGGAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGa                                                                                                     <  5:70496/146‑1 (MQ=255)
                                                    cAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCg                                                                                                                                     >  2:18593/1‑115 (MQ=255)
                                                    cAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGa                                                                                                     <  4:159775/147‑1 (MQ=255)
                                                     aaGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑Ac                                                                                               <  5:122701/149‑1 (MQ=255)
                                                     aaGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑Ac                                                                                               <  7:63827/149‑1 (MQ=255)
                                                        gtggtgTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTc                                                                                            <  5:358153/149‑1 (MQ=255)
                                                           gtgtTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTcgcg                                                                                         >  4:140213/1‑149 (MQ=255)
                                                             gtTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCt                                                                                       >  7:143262/1‑149 (MQ=255)
                                                                 gATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCt                                                                                    <  6:344358/148‑1 (MQ=255)
                                                                   ttATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCaa                                                                                                              <  5:156417/123‑1 (MQ=255)
                                                                   ttATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCaa                                                                                                              >  6:156417/1‑123 (MQ=255)
                                                                   ttATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTgg                                                                                  >  1:55426/1‑148 (MQ=255)
                                                                        gTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAAAGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTACTGGCAggtg                                                                            >  2:307468/1‑149 (MQ=255)
                                                                         tGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAggtgg                                                                           <  5:13524/149‑1 (MQ=255)
                                                                         tGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAggtgg                                                                           >  2:241046/1‑149 (MQ=255)
                                                                         tGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCAAATGTCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAggtgg                                                                           >  7:238051/1‑149 (MQ=255)
                                                                           aCGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAggtggt                                                                          <  1:250569/148‑1 (MQ=255)
                                                                            cGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGtggtgg                                                                        <  6:217169/149‑1 (MQ=255)
                                                                             ggCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGa                                                                       >  8:182756/1‑149 (MQ=255)
                                                                               cGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACa                                                                     <  1:241046/149‑1 (MQ=255)
                                                                                  ggagtTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGc                                                                  <  3:200397/146‑1 (MQ=255)
                                                                                    tgtTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCAc                                                                <  4:99812/149‑1 (MQ=255)
                                                                                        cTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAgg                                                            >  2:164273/1‑149 (MQ=255)
                                                                                           gTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGt                                                           <  6:214361/147‑1 (MQ=255)
                                                                                              cACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGGCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACa                                                      >  1:366389/1‑149 (MQ=255)
                                                                                                    cgATGTCGGCTGGGGGACCGGCCACATCTACATCCTCTAATGCCCGCTTGCCAATTGCTGGAATCTCGCTGATTTTCGAATGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACAtcttc                                                 <  1:308002/149‑1 (MQ=255)
                                                                                                                  tGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGAcc                                  >  6:328925/1‑149 (MQ=255)
                                                                                                                   gACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCg                                 >  3:331781/1‑149 (MQ=255)
                                                                                                                           acaGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCTTGCCCAACTACCCAGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCATgcgc                         <  3:73438/149‑1 (MQ=255)
                                                                                                                            caGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGATCTCGCCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTgc                           <  1:307468/149‑1 (MQ=38)
                                                                                                                               cTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGCCCGCCCTGCGCGCa                      >  8:141192/1‑148 (MQ=255)
                                                                                                                               cTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCAc                     <  6:331820/149‑1 (MQ=255)
                                                                                                                               cTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCAc                     <  6:105363/149‑1 (MQ=255)
                                                                                                                                 aCATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGTACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACt                    >  4:284786/1‑148 (MQ=255)
                                                                                                                                     cGTCTATGGACAGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATGc               >  4:184333/1‑149 (MQ=255)
                                                                                                                                      gccattGGCCAGCGGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTATCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATg                <  1:20264/142‑1 (MQ=255)
                                                                                                                                      gccattGGCCAGCGGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTATCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATg                >  2:20264/6‑147 (MQ=255)
                                                                                                                                      gTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATGc               <  3:140213/148‑1 (MQ=255)
                                                                                                                                       tCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATGCGt             <  5:138294/149‑1 (MQ=255)
                                                                                                                                        cattGGCCAGCGGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATGCGt             >  8:73503/4‑148 (MQ=255)
                                                                                                                                           tGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATGCGTGa           <  7:182756/147‑1 (MQ=255)
                                                                                                                                                  cTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATGCGTGAGGGCTCcgc  <  8:143262/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
AAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATGCGTGAGGGCTCCGC  >  NC_002947/5345405‑5345699

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: