Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,083,051 G→C 6.2% G17A (GGT→GCT)  PP_3593 → amino acid ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,083,0510GC 6.2% 203.7 / 3.7 64G17A (GGT→GCT) PP_3593amino acid ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base C (2/2);  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 = 4.35e-01

ACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGC  >  NC_002947/4082904‑4083194
                                                                                                                                                   |                                                                                                                                               
aCGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGt                                                                                                                                                <  8:55821/149‑1 (MQ=255)
    ttGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTgccg                                                                                                                                            >  6:325054/1‑149 (MQ=255)
            aaCTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCt                                                                                                                                    >  5:373029/1‑149 (MQ=255)
              cTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGTGCAAACAACAATAGCAAGAGTTCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGC‑GCAGGCGGCACCCGCCAgcagcagg                                                                                                                                 <  7:254233/149‑7 (MQ=255)
                 ccGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAgcgtcgggtg                                                                                                                               >  6:167715/1‑142 (MQ=255)
                 ccGCTCTACATAACGAATTCGTCGGCTACACCTCCTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGATGTGCGAGATGAACAAGTCCATGGCCTTGCTGTGTGCGTGCGCCCTGCTGCCGGCGGCACCCGCCAgcagcagg                                                                                                                                 <  7:59692/147‑7 (MQ=255)
                    ctctACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCGATGGCCTTGATTGGTGCGTGCGCGCTGCTGCTGGCGGGTGCCGGCAGCGCGGAACCCCt                                                                                                                            >  8:387297/1‑149 (MQ=255)
                    ctctACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCt                                                                                                                            >  8:195059/1‑149 (MQ=255)
                         cATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGt                                                                                                                       >  4:227261/1‑149 (MQ=255)
                         cATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGGGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCg                                                                                                                         <  6:24675/147‑1 (MQ=255)
                          aTAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGtt                                                                                                                      >  4:45736/1‑149 (MQ=255)
                          aTAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTCCAATGAGGTGCGAGATGAAAAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCGGAAACCCTGCGGtt                                                                                                                      >  5:343381/1‑149 (MQ=255)
                               gAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCa                                                                                                                 >  1:124558/1‑149 (MQ=255)
                                 aTTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCAcc                                                                                                               >  5:361048/1‑149 (MQ=255)
                                    cgtcgGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAg                                                                                                            >  4:92325/1‑149 (MQ=255)
                                     gtcgGCTACACCTCGTGGTCGACGTGCAAACAACAATAGCAAGATTGCAATTAGGTGCGATATGAACAAGTCAATGGCCTTGCTTGGTGCGTGCGCCCTGGTGCTGGCGGGTGCCGGCAGCGCTGAAACCCTTCGTTTCGCCACCGatg                                                                                                           <  6:343381/149‑3 (MQ=255)
                                     gtcgGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAg                                                                                                            <  3:227261/148‑1 (MQ=255)
                                        ggCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTgcc                                                                                                                                             >  1:44505/1‑112 (MQ=255)
                                        ggCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTgcc                                                                                                                                             <  2:44505/112‑1 (MQ=255)
                                        ggCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAg                                                                                                                                        <  6:206837/117‑1 (MQ=255)
                                        ggCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAg                                                                                                                                        >  5:206837/1‑117 (MQ=255)
                                                cTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCg                                                                                                >  4:389129/1‑149 (MQ=255)
                                                    tGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCt                                                                                            <  5:227460/149‑1 (MQ=255)
                                                         gACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGGGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAAc                                                                                       >  1:236920/1‑149 (MQ=255)
                                                         gACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAAc                                                                                       <  7:195059/149‑1 (MQ=255)
                                                         gACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAAc                                                                                       >  5:309132/1‑149 (MQ=255)
                                                             ggCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATg                                                                                   <  3:389129/149‑1 (MQ=255)
                                                                aaaCAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATg                                                                                   >  2:324589/1‑146 (MQ=255)
                                                                    aacaaTAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGc                                                                            >  2:215153/1‑149 (MQ=255)
                                                                      caaTAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCg                                                                          <  1:215153/149‑1 (MQ=255)
                                                                            cAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCATCTATGTAGACGCTGATAAc                                                                     >  2:382592/1‑148 (MQ=255)
                                                                              agagTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCATCTATGTAGACGCTGATAACAAg                                                                  <  1:382592/149‑1 (MQ=255)
                                                                                    cAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCAc                                                            <  3:45736/149‑1 (MQ=255)
                                                                                     aaTGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACg                                                           >  2:59951/1‑149 (MQ=255)
                                                                                          ggTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCttt                                                      <  2:236920/149‑1 (MQ=255)
                                                                                               gagaTGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATg                                                  >  4:240304/1‑148 (MQ=255)
                                                                                               gagaTGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGt                                                 >  7:25033/1‑149 (MQ=255)
                                                                                                agaTGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCGTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTc                                                <  2:130724/149‑1 (MQ=255)
                                                                                                agaTGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTc                                                <  1:31914/149‑1 (MQ=255)
                                                                                                agaTGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATATc                                                <  2:124558/149‑1 (MQ=255)
                                                                                                  aTGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACa                                                                    >  2:388784/1‑127 (MQ=255)
                                                                                                  aTGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACa                                                                    <  1:388784/127‑1 (MQ=255)
                                                                                                   ccagcAAGGCCATGGACTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGa                                              >  6:351273/5‑148 (MQ=255)
                                                                                                    gAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACa                                            >  3:65851/1‑149 (MQ=255)
                                                                                                     agcAAGGCCATGGACTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACAt                                           >  7:356105/3‑149 (MQ=255)
                                                                                                          gTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCa                                         <  4:93559/146‑1 (MQ=255)
                                                                                                          gTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCa                                         >  3:93559/1‑146 (MQ=255)
                                                                                                             cATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATg                                   >  5:215366/1‑149 (MQ=255)
                                                                                                                  ccTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCtg                              >  4:211932/1‑149 (MQ=255)
                                                                                                                   cTTGCTGGGTGCGTGCGCCCTGCTGCTGGCTGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACGCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCtgt                             <  5:325054/149‑1 (MQ=255)
                                                                                                                   cTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCtgt                             <  3:92325/149‑1 (MQ=255)
                                                                                                                    ttGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCtgtg                            >  7:356280/1‑149 (MQ=255)
                                                                                                                         gggTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAAc                       <  6:373029/149‑1 (MQ=255)
                                                                                                                          ggTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAAc                       <  3:240304/148‑1 (MQ=255)
                                                                                                                          ggTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTCAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGaa                        <  5:351273/147‑1 (MQ=255)
                                                                                                                             gcgtgcgCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGAt                   <  2:325322/149‑1 (MQ=255)
                                                                                                                              cgtgcgCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATg                  <  6:215366/149‑1 (MQ=255)
                                                                                                                                    cccTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTc            <  3:211932/149‑1 (MQ=255)
                                                                                                                                         gcgctgGCGGCACCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATg       >  2:149325/3‑149 (MQ=255)
                                                                                                                                         gcgctgGCGGCACCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATg       >  5:42910/3‑149 (MQ=255)
                                                                                                                                         ctgctgGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAAt        <  8:25033/148‑1 (MQ=255)
                                                                                                                                            ctgGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGc                                  <  7:311491/119‑1 (MQ=255)
                                                                                                                                            ctgGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGc                                  >  8:311491/1‑119 (MQ=255)
                                                                                                                                               gCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGc  <  8:356280/148‑1 (MQ=255)
                                                                                                                                                 gggTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGAc                                             <  5:308721/103‑1 (MQ=255)
                                                                                                                                                 gggTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGAc                                             >  6:308721/1‑103 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGC  >  NC_002947/4082904‑4083194

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