Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,507,016 T→G 5.9% W93G (TGG→GGG)  cycA → D‑alanine/beta‑alanine/D‑serine/glycine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,507,0160TG 5.9% 214.9 / 4.1 68W93G (TGG→GGG) cycAD‑alanine/beta‑alanine/D‑serine/glycine permease
Reads supporting (aligned to +/- strand):  ref base T (37/26);  new base G (2/2);  total (40/28)
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.70e-01

CTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAGGGCTGCTG  >  NC_002947/5506875‑5507162
                                                                                                                                              |                                                                                                                                                  
cTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCgg                                                                                                                                             <  3:190159/149‑1 (MQ=255)
 gggCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCAGGGGCCgag                                                                                                                                            <  7:51977/148‑3 (MQ=255)
        gCCATCGATCATCTTCGTCTACATGATTATCGGCTTCCTGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTa                                                                                                                                     <  1:331437/149‑1 (MQ=255)
          cATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTTTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTAct                                                                                                                                   <  1:126517/149‑1 (MQ=255)
          cATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTAct                                                                                                                                   <  5:220628/149‑1 (MQ=255)
          cATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCAGGGGCCgagc                                                                                                                                           <  8:147906/141‑4 (MQ=255)
          cATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCAGGGGCCgagc                                                                                                                                           >  7:147906/1‑138 (MQ=255)
           aTCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACt                                                                                                                                        <  4:117620/143‑1 (MQ=255)
           aTCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACt                                                                                                                                        >  3:117620/1‑143 (MQ=255)
               atcatcTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTActggctg                                                                                                                              <  5:104489/149‑1 (MQ=255)
                tcatcTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTActggctg                                                                                                                              >  1:296237/1‑148 (MQ=255)
                   tcTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCAGACGACCTGATCGGCCCCTGGGCCGGCTACTTTACTGGCTGGAcg                                                                                                                          >  4:299616/1‑148 (MQ=255)
                       cGTCTACATGATTCTCGGCTTCCTGCTGTTCATCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTa                                                                                                                        <  4:101736/147‑1 (MQ=255)
                       cGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTAGAACCTGAACTACAAGTCGTTCATCGTCTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGAcc                                                                                                                          <  4:119458/145‑1 (MQ=255)
                        gTCTACATGATTATCGGCTGCATGCTGGTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGGCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTg                                                                                                                     <  8:208288/149‑1 (MQ=255)
                            aCATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTc                                                                                                                 >  4:341794/1‑149 (MQ=255)
                                   tATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCt                                                                                                                                                    <  4:218405/107‑1 (MQ=255)
                                   tATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCt                                                                                                                                                    >  3:213722/1‑107 (MQ=255)
                                   tATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCt                                                                                                                                                    >  3:218405/1‑107 (MQ=255)
                                   tATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCt                                                                                                                                                    <  4:213722/107‑1 (MQ=255)
                                     tCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTAGTGGTTCTGCTGggtg                                                                                                        >  5:296949/1‑149 (MQ=255)
                                     tCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGggtg                                                                                                        >  2:171671/1‑149 (MQ=255)
                                       ggCTTCATGCTGTTCTTCGTCATGCGCGCCAATGGGCGAACTGCTGCTGTCGAACCTGAAATACAAGTCGTACATCGACTTTTCCGCCGACATGATCGGCCCATGGGCCGGCTACTTTACTGCATGGACGTACTGGTTATGATGggtgg                                                                                                       >  1:39489/1‑149 (MQ=255)
                                         cTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTctgctg                                                                                                            <  7:19137/141‑1 (MQ=255)
                                         cTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTctgctg                                                                                                            >  8:19137/1‑141 (MQ=255)
                                           tCATCCTCTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGATGTCGAACCTGAACTACAAGTAGTACATCGACATTTACGCAGACATGCTAGGCCCCACGGCCCACTAAATTAATTGCTGGACCTACAGCTTCTGCTTGGTGGTAAcc                                                                                                  >  4:104686/1‑149 (MQ=255)
                                               gCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTActggctg                                                                                                                              <  7:321177/117‑1 (MQ=255)
                                               gCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTActggctg                                                                                                                              >  8:321177/1‑117 (MQ=255)
                                               gCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTctg                                                                                                               <  4:90008/132‑1 (MQ=255)
                                               gCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTctg                                                                                                               >  3:90008/1‑132 (MQ=255)
                                               gCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCa                                                                                              >  6:99000/1‑149 (MQ=255)
                                                cTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCa                                                                                              <  2:236518/148‑1 (MQ=255)
                                                 tGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTCAACTACAAGTCGTTCATCGACTTTTACGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCAtt                                                                                            >  8:37511/1‑149 (MQ=255)
                                                  gTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTg                                                                                           >  5:336315/1‑149 (MQ=255)
                                                   ttcttcGGCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCAtt                                                                                            <  6:332408/147‑1 (MQ=255)
                                                     cttcGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGcc                                                                                         >  8:8200/1‑148 (MQ=255)
                                                            aTGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCAGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGt                                                                                 <  6:200165/149‑1 (MQ=255)
                                                                                   tgctgTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAg                                                           >  6:167178/1‑149 (MQ=255)
                                                                                   tgctgTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAg                                                           >  6:338619/1‑149 (MQ=255)
                                                                                     ctgTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTa                                                                                                                                     <  5:375893/73‑1 (MQ=255)
                                                                                     ctgTCGAACCTGAACTACAAGTCGATCATCGACTTTTCCGCCGACCTGCTCGGCCCCTCGGCCCGCTACTTTa                                                                                                                                     >  6:375893/1‑73 (MQ=255)
                                                                                        tCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTg                                                      <  1:7778/149‑1 (MQ=255)
                                                                                        tCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTg                                                      <  5:338619/149‑1 (MQ=255)
                                                                                          gAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTgg                                                     <  2:10427/148‑1 (MQ=255)
                                                                                          gAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGt                                                    >  2:160799/1‑149 (MQ=255)
                                                                                          gAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGt                                                    >  2:283228/1‑149 (MQ=255)
                                                                                             ccTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTcc                                                 <  2:296237/149‑1 (MQ=255)
                                                                                               tGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGTTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTcccc                                               >  7:115628/1‑149 (MQ=255)
                                                                                               tGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTcccc                                               >  3:362576/1‑149 (MQ=255)
                                                                                               tGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTcccc                                               >  1:339157/1‑149 (MQ=255)
                                                                                                gAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTcccc                                               >  6:251976/1‑148 (MQ=255)
                                                                                                  cggaaaaGTCGATGAACGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGAc                                            >  7:5424/6‑149 (MQ=255)
                                                                                                  cggaaaaGTCGATGAACGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGAc                                            >  7:196555/6‑149 (MQ=255)
                                                                                                   ggaaaaGTCGATGAACGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGAcc                                           >  6:346268/5‑149 (MQ=255)
                                                                                                    tACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCt                                          >  3:113764/1‑149 (MQ=255)
                                                                                                            gTTCATCGACTTTTCCGCCGACCTGTTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGt                                  >  7:101476/1‑149 (MQ=255)
                                                                                                            gTTCATCGACTTTTCCGCCGACCTGTTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGt                                  <  8:101476/149‑1 (MQ=255)
                                                                                                            gTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAg                                   <  3:104686/148‑1 (MQ=255)
                                                                                                            gTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGt                                  >  5:188834/1‑149 (MQ=255)
                                                                                                               cATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTgg                                <  3:341794/148‑1 (MQ=255)
                                                                                                                      ttttCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGgc                        >  8:215726/1‑149 (MQ=255)
                                                                                                                        ttCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGgcgc                      >  3:68029/1‑149 (MQ=255)
                                                                                                                           cgccgACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGa                   >  2:36420/1‑149 (MQ=255)
                                                                                                                           cgccgACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTCGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACCTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGa                   >  8:249480/1‑149 (MQ=255)
                                                                                                                             ccgACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATcgccgc                                                                     <  2:159677/97‑1 (MQ=255)
                                                                                                                                     gccGGCCCAGGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAgg         >  2:248892/3‑149 (MQ=255)
                                                                                                                                      cgatctccTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATCGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAggg        <  7:136343/143‑1 (MQ=255)
                                                                                                                                        ggCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAggg        <  5:41812/147‑1 (MQ=255)
                                                                                                                                         gCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAGGgctg     >  1:335451/1‑149 (MQ=255)
                                                                                                                                           cccTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAGGgctgct   >  6:368155/1‑149 (MQ=255)
                                                                                                                                            ccTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCTTTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAGGgctgctg  >  1:103457/1‑149 (MQ=255)
                                                                                                                                             cTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATCGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGCTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCCTAggg        >  8:136343/1‑142 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCC‑ATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAGGGCTGCTG  >  NC_002947/5506875‑5507162

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