Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,083,053 G→C 6.4% A18P (GCC→CCC)  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,0530GC 6.4% 199.8 / 3.7 62A18P (GCC→CCC) PP_3593amino acid ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (30/28);  new base C (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

TTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGC  >  NC_002947/4082908‑4083194
                                                                                                                                                 |                                                                                                                                             
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                                                                                                                                             <  2:44505/112‑1 (MQ=255)
                                    ggCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTgcc                                                                                                                                             >  1:44505/1‑112 (MQ=255)
                                    ggCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAg                                                                                                                                        >  5:206837/1‑117 (MQ=255)
                                    ggCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAg                                                                                                                                        <  6:206837/117‑1 (MQ=255)
                                            cTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCg                                                                                                >  4:389129/1‑149 (MQ=255)
                                                tGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCt                                                                                            <  5:227460/149‑1 (MQ=255)
                                                     gACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGGGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAAc                                                                                       >  1:236920/1‑149 (MQ=255)
                                                     gACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAAc                                                                                       >  5:309132/1‑149 (MQ=255)
                                                     gACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAAc                                                                                       <  7:195059/149‑1 (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                                                                    <  1:388784/127‑1 (MQ=255)
                                                                                              aTGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACa                                                                    >  2:388784/1‑127 (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       >  5:42910/3‑149 (MQ=255)
                                                                                                                                     gcgctgGCGGCACCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATg       >  2:149325/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)
                                                                                                                                                 |                                                                                                                                             
TTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGC  >  NC_002947/4082908‑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: ‑