Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,719,255 Δ1 bp 8.3% coding (382/585 nt) nfuA → Fe/S biogenesis protein NfuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,719,2550A. 8.3% 184.3 / 11.0 60coding (382/585 nt)nfuAFe/S biogenesis protein NfuA
Reads supporting (aligned to +/- strand):  ref base A (29/26);  new base . (3/2);  total (32/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 = 8.24e-01

CGCGCTGGTCGACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATT  >  NC_002947/2719107‑2719395
                                                                                                                                                      |                                                                                                                                             
cgcgCTGGTCGACGACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGgcca                                                                                                                                               <  5:210868/149‑1 (MQ=255)
           aCTACCCCACCCTCCGGCTTGGTGGCCATCTGACCATCATTGCGCCAAACGGCAATGTGCCGATTGGCAACGTATTCAGCCCGTTCTCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATggcg                                                                                                                                     <  1:295354/149‑1 (MQ=255)
                  cACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCGCCATGGCTGGCCAggccg                                                                                                                            >  6:124205/1‑146 (MQ=255)
                     cgaccgCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCGCCATGGCTGGCCAggccgggg                                                                                                                         <  5:124205/149‑7 (MQ=255)
                            aTGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGa                                                                                                                  <  8:209044/149‑1 (MQ=255)
                               ggTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGt                                                                                                               <  3:22931/149‑1 (MQ=255)
                                   gCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGt                                                                                                            <  8:6733/148‑1 (MQ=255)
                                      aGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGAt                                                                                                        >  1:18564/1‑149 (MQ=255)
                                           aCCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACg                                                                                                    <  2:253705/148‑1 (MQ=255)
                                               tCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAACACAGCCCG‑‑AACAACGAACGAATAAACTACTACCTACAGACCGAGATCAAACCCGGCCAGGCCAGCCATCGCGGCCAGGAGAGCCTGGTGCCCGTGGTCGATGCCGTCAtt                                                                                               >  8:108597/1‑149 (MQ=255)
                                                cAATGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTg                                                                                              <  3:181131/149‑1 (MQ=255)
                                                cAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTg                                                                                              >  4:148019/1‑149 (MQ=255)
                                                     cgcCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAg                                                                                                                           <  7:327324/115‑1 (MQ=255)
                                                     cgcCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAg                                                                                                                           >  8:327324/1‑115 (MQ=255)
                                                      gcCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTAc                                                                                        >  6:245383/1‑149 (MQ=255)
                                                        cAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTg                                                                                      <  2:18564/149‑1 (MQ=255)
                                                        cAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTg                                                                                      >  4:330045/1‑149 (MQ=255)
                                                         aaaCGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTg                                                                                      >  5:89932/1‑148 (MQ=255)
                                                          aaCGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCa                                                                                    <  3:330045/149‑1 (MQ=255)
                                                                 aGGTGCCGTTGCTCAACTACTACCGCCCG‑‑AGCAACGAACGTAGCAACTACTACCTCCAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGGCGATGACGGCATTGCCGTACTGCAGTTCggt                                                                             <  6:101073/149‑1 (MQ=255)
                                                                          tGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGc                                                                    >  8:62702/1‑149 (MQ=255)
                                                                            gTCAACGATGACAGCCCG‑‑ATCATCTAACGTATCAACTAGTACCTTCAGACCTAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCa                                                                  <  4:26723/149‑1 (MQ=255)
                                                                             tCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAg                                                                 >  6:332750/1‑149 (MQ=255)
                                                                               aaCGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAggg                                                               >  4:121638/1‑149 (MQ=255)
                                                                                aCGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGt                                                              >  8:155477/1‑149 (MQ=255)
                                                                                aCGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGt                                                              <  3:148019/149‑1 (MQ=255)
                                                                                 cGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGtt                                                             >  8:241400/1‑149 (MQ=255)
                                                                                    aGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCg                                                          >  3:114375/1‑149 (MQ=255)
                                                                                     gACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCgg                                                         <  7:346174/149‑1 (MQ=255)
                                                                                      aCAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGc                                                        >  3:15347/1‑149 (MQ=255)
                                                                                        aGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCa                                                      <  7:62702/149‑1 (MQ=255)
                                                                                            cG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGcc                                                  >  6:95732/1‑149 (MQ=255)
                                                                                                 tCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTTGTGGACGTGGTCGATGCCGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGAt                                               <  7:108597/149‑1 (MQ=255)
                                                                                                   aacgaacgTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCg                                                          <  3:217772/136‑1 (MQ=255)
                                                                                                   aacgaacgTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCg                                                          >  4:217772/1‑136 (MQ=255)
                                                                                                     cgaacgTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTTGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCa                                           <  4:15347/149‑1 (MQ=255)
                                                                                                     cgaacgTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTCCCAGGGTTGCGGCCAGGCCGATGTCa                                           <  7:155477/149‑1 (MQ=255)
                                                                                                          gTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATCTCACCCTg                                      <  6:56611/149‑1 (MQ=255)
                                                                                                               aaCTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCtg                                                                                                                       >  4:45046/1‑63 (MQ=255)
                                                                                                               aaCTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCtg                                                                                                                       <  3:45046/63‑1 (MQ=255)
                                                                                                                actactACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaagg                                  >  3:342517/1‑147 (MQ=255)
                                                                                                                actactACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaagg                                  >  3:337224/1‑147 (MQ=255)
                                                                                                                 ctactaCCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaaggaag                               <  5:245383/149‑1 (MQ=255)
                                                                                                                 ctactaCCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaaggaag                               >  6:281773/1‑149 (MQ=255)
                                                                                                                 ctactaCCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaaggaag                               >  6:361343/1‑149 (MQ=255)
                                                                                                                    ctacctacAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGc                             <  2:207964/148‑1 (MQ=255)
                                                                                                                     tacctacAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCAt                           >  5:70357/1‑149 (MQ=255)
                                                                                                                            aGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGcc                                                       <  2:350026/114‑1 (MQ=255)
                                                                                                                            aGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGcc                                                       >  1:350026/1‑114 (MQ=255)
                                                                                                                               ccGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCAcc                 <  2:41613/149‑1 (MQ=255)
                                                                                                                                cGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaa                                    >  6:372798/1‑129 (MQ=255)
                                                                                                                                cGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTgaa                                    <  5:372798/129‑1 (MQ=255)
                                                                                                                                cGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAg                       >  7:44080/1‑142 (MQ=255)
                                                                                                                                cGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAg                       <  8:44080/142‑1 (MQ=255)
                                                                                                                                cGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCAccc                >  7:296015/1‑149 (MQ=255)
                                                                                                                                           ccGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCg      >  8:262239/1‑148 (MQ=255)
                                                                                                                                            tcaccTGGCCGCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGc       >  5:288428/4‑146 (MQ=255)
                                                                                                                                            tcaccTGGCCGCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGc       <  6:288428/143‑1 (MQ=255)
                                                                                                                                            cGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTa    >  1:194645/1‑149 (MQ=255)
                                                                                                                                              gCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGCCGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCCAGCGCACACTGCTCGAGCGTAtt  >  7:253425/1‑149 (MQ=255)
                                                                                                                                              accTGGCCGCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTAtt  >  1:292013/2‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                             
CGCGCTGGTCGACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCG‑‑ATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATT  >  NC_002947/2719107‑2719395

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: