Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,780,565 C→T 6.4% A339V (GCC→GTC)  nikA → nickel ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,780,5650CT 6.4% 180.6 / 2.5 62A339V (GCC→GTC) nikAnickel ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (25/33);  new base T (2/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

AACGAGGCGGTGGACAAGCAAGCCATCATCGACAAGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGC  >  NC_002947/3780420‑3780702
                                                                                                                                                 |                                                                                                                                         
aaCGAGGCGGTGGACAAGCAAGCCATCATCGACAAGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCt                                                                                                                                        <  8:114791/149‑1 (MQ=255)
aaCGAGGCGGTGGACAAGCAAGCCATCATCGACAAGATTCTATACGCCCTTGAACCCCGGCCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCAAGCCCt                                                                                                                                        >  1:282558/1‑149 (MQ=255)
  cGAGGCGGTGGACAAGCAAGCCATCATCGACAAGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCtt                                                                                                                                       >  3:220119/1‑148 (MQ=255)
  cGAGGCGGTGGACAAGCAAGCCATCATCGACAAGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTg                                                                                                                                      <  6:23634/149‑1 (MQ=255)
  cGAGGCGGTGGACAAGCAAGCCATCATCGACAAGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTg                                                                                                                                      <  3:264918/149‑1 (MQ=255)
   gAGGCGGTGGACAAGCAAGCCATCATCGACAAGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGa                                                                                                                                     <  4:170154/149‑1 (MQ=255)
    aGGCGGTGGACAAGCAAGCCATCATCGACAAGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGAt                                                                                                                                    <  1:18558/149‑1 (MQ=255)
      tcGGTGGACAAGCAAGCCATCATCGACAAGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGAtgc                                                                                                                                  <  5:168388/148‑1 (MQ=255)
              caagcaagCCATCATCGACAAGATTCTCTACGGCCTTGAACCCCTGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGt                                                                                                                          <  6:171778/149‑1 (MQ=255)
              caagcaagCCATCATCGACAAGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGACCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATCCTGcccgct                                                                                                                          >  8:299335/1‑145 (MQ=255)
                 gcaagcCATCATCGACAAGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGa                                                                                                                       >  2:164970/1‑149 (MQ=255)
                              gACAAGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAgg                                                                                                          >  8:274991/1‑149 (MQ=255)
                              gACAAGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAg                                                                                                           <  6:298050/148‑1 (MQ=255)
                              gACAAGATTATCTACGGCCTTGAACCCCGGGGAGATACCCTGATCGCCAGCAACATGCCCTACGCCGATATCGTCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGTGTGGAGGCTGGCCCCAg                                                                                                           <  8:114047/148‑1 (MQ=255)
                               aCAAGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGgc                                                                                                         >  7:58212/1‑149 (MQ=255)
                               aCAAGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCAGCCGGGTGGACGCTGGCCCCAGgc                                                                                                         <  3:284592/149‑1 (MQ=255)
                                 aaGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGgcgc                                                                                                       <  6:28535/149‑1 (MQ=255)
                                     ttCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCt                                                                                                   >  1:226971/1‑149 (MQ=255)
                                      tctcTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGTCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCtg                                                                                                  >  7:159977/1‑149 (MQ=255)
                                      tctcTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGTCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCtg                                                                                                  >  3:3314/1‑149 (MQ=255)
                                        tctACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCtgtg                                                                                                >  7:16587/1‑149 (MQ=255)
                                         ctACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGc                                                                                                                >  7:260378/1‑132 (MQ=255)
                                         ctACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGc                                                                                                                <  8:260378/132‑1 (MQ=255)
                                              gCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAg                                                                                           >  4:145379/1‑148 (MQ=255)
                                               ccTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTa                                                                                         <  2:226971/149‑1 (MQ=255)
                                                cTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCa                                                                                            >  4:75248/1‑145 (MQ=255)
                                                cTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCa                                                                                            <  3:75248/145‑1 (MQ=255)
                                                   gAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTaaa                                                                                       <  4:175067/147‑1 (MQ=255)
                                                      ccccGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAggcgg                                                                                  <  7:274991/149‑1 (MQ=255)
                                                        ccGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGgcggcg                                                                                >  8:110984/1‑149 (MQ=255)
                                                                 tACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGCTGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTgg                                                                       <  2:238803/149‑1 (MQ=255)
                                                                           gCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCt                                                             >  8:14638/1‑149 (MQ=255)
                                                                              aGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCa                                                                                            <  5:99405/115‑1 (MQ=255)
                                                                              aGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCa                                                                                            >  6:99405/1‑115 (MQ=255)
                                                                                cAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTAt                                                        <  8:58212/149‑1 (MQ=255)
                                                                                   cATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTcc                                                     >  3:277862/1‑149 (MQ=255)
                                                                                    aTGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCt                                                    >  1:167158/1‑149 (MQ=255)
                                                                                      gCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCg                                                  >  2:48942/1‑149 (MQ=255)
                                                                                       cccTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGTCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCgg                                                 <  8:159977/149‑1 (MQ=255)
                                                                                       cccTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGTCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCgg                                                 <  4:3314/149‑1 (MQ=255)
                                                                                       cccTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCg                                                  <  8:16587/148‑1 (MQ=255)
                                                                                       cccTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCCGCGCGCGCCGCCCATGCCCTTGCTGCTGCCGGGTGGACTCTGGCCCCAGGCGCGTCTGTGCGCAGTATACGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCgg                                                 <  2:201482/149‑1 (MQ=255)
                                                                                         cTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGca                                               <  1:164970/149‑1 (MQ=255)
                                                                                          tACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGca                                               >  2:200465/1‑148 (MQ=255)
                                                                                             gCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCg                                                  <  5:12278/142‑1 (MQ=255)
                                                                                             gCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCg                                                  >  6:12278/1‑142 (MQ=255)
                                                                                                 atatCGGCCTCAAGCGCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGcc                                       <  7:299335/149‑1 (MQ=255)
                                                                                                    tCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGGTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCtt                                     <  2:282558/148‑1 (MQ=255)
                                                                                                     cGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTg                                    >  2:154193/1‑148 (MQ=255)
                                                                                                                 cTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTgg                       >  1:296109/1‑149 (MQ=255)
                                                                                                                  tACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGc                      >  6:129532/1‑149 (MQ=255)
                                                                                                                   aCCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGcc                     >  2:136296/1‑149 (MQ=255)
                                                                                                                    cccTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGc                      <  1:200465/147‑1 (MQ=255)
                                                                                                                     ccTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGa                   <  1:296062/149‑1 (MQ=255)
                                                                                                                         aCGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGt                <  2:296109/148‑1 (MQ=255)
                                                                                                                            ccccGGCGCGCGCCGCCCAGGCCCTTGATGCGGCGGGGTGGACGCTGGCCCCAGGGGCTTCTGTGCGCAGTAAAGGCGGCTAACCCTTGGCACTCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGtgctg            <  1:252381/148‑1 (MQ=255)
                                                                                                                            aCCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGtgctg            >  3:58922/1‑149 (MQ=255)
                                                                                                                             cccGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGtgctgc           <  7:14638/149‑1 (MQ=255)
                                                                                                                                ggCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAgg        <  1:154193/149‑1 (MQ=255)
                                                                                                                                 gcgcgcgcCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAggg       <  1:48942/149‑1 (MQ=255)
                                                                                                                                    cgcgcCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGa    <  4:58922/149‑1 (MQ=255)
                                                                                                                                      cgccgcCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGc  >  5:254105/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
AACGAGGCGGTGGACAAGCAAGCCATCATCGACAAGATTCTCTACGGCCTTGAACCCCGGGCAGATACCCTGTTCGCCAGCAACATGCCCTACGCCGATATCGGCCTCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGC  >  NC_002947/3780420‑3780702

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