Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,160 C→T 7.1% E408K (GAG→AAG)  trmF ← tRNA uridine 5‑carboxymethylaminomethyl modification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,1600CT 7.1% 163.0 / 2.4 56E408K (GAG→AAG) trmFtRNA uridine 5‑carboxymethylaminomethyl modification protein
Reads supporting (aligned to +/- strand):  ref base C (26/26);  new base T (2/2);  total (28/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.74e-01

CCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGCCCGGCGAAGAACAGGCCGCCGATGACTTTC  >  NC_002947/3025‑3302
                                                                                                                                         |                                                                                                                                              
ccTTTTCGGTCAGGCTCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCggg                                                                                                                                   <  4:95787/149‑1 (MQ=255)
ccTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCggg                                                                                                                                   <  8:110342/149‑1 (MQ=255)
 cTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCaggg                                                                                                                                  >  7:157021/1‑148 (MQ=255)
        gTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCtcgt                                                                                                                                              >  6:166082/1‑130 (MQ=255)
        gTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCtcgt                                                                                                                                              <  5:166082/130‑1 (MQ=255)
            ggCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCa                                                                                                                       <  2:311048/149‑1 (MQ=255)
             gcgcAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTTGTCGCGGCGCGGGCACCAGCTGTCAc                                                                                                                      >  6:92309/1‑149 (MQ=255)
             gcgcAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTTGTCGCGGCGCGGGCACCAGCTGTCAc                                                                                                                      >  6:96924/1‑149 (MQ=255)
             gcgcAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCAc                                                                                                                      <  1:34746/149‑1 (MQ=255)
             gcgcAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCAc                                                                                                                      >  2:217276/1‑149 (MQ=255)
              cgcAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAATCCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCAc                                                                                                                      <  1:271780/149‑1 (MQ=255)
                 aGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGcc                                                                                                                  <  4:279468/149‑1 (MQ=255)
                       gCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTgcgc                                                                                                            <  7:164309/149‑1 (MQ=255)
                        cGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCTAAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCtcgt                                                                                                                                              <  1:67150/115‑1 (MQ=255)
                        cGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCtcgt                                                                                                                                              >  2:67150/1‑114 (MQ=255)
                            gTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTTGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGa                                                                                                        <  5:92309/148‑1 (MQ=255)
                            gTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTTGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGa                                                                                                        <  5:96924/148‑1 (MQ=255)
                              cTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGt                                                                                                     >  2:11368/1‑149 (MQ=255)
                                tCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTgc                                                                                                   >  8:81206/1‑149 (MQ=255)
                                 cgcgCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGt                                                                                                     <  1:11368/146‑1 (MQ=255)
                                     cAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGt                                                                                              <  3:250095/149‑1 (MQ=255)
                                     cAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGt                                                                                              <  4:151791/149‑1 (MQ=255)
                                               cGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGc                                                                                    >  7:280716/1‑149 (MQ=255)
                                                            gcgAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTgg                                                                       >  3:10088/1‑149 (MQ=255)
                                                              gAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGt                                                                                              <  7:308392/124‑1 (MQ=255)
                                                              gAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGt                                                                                              >  8:308392/1‑124 (MQ=255)
                                                               aGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGgcg                                                                    >  1:207513/1‑149 (MQ=255)
                                                                 gTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCcagcag                                                                              <  4:183996/137‑1 (MQ=255)
                                                                 gTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCcagcag                                                                              >  3:183996/1‑137 (MQ=255)
                                                                    aaCATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTgc                                                                                                   >  7:197967/1‑113 (MQ=255)
                                                                    aaCATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTgc                                                                                                   <  8:197967/113‑1 (MQ=255)
                                                                          ctgttcGGCTCCTGTCGACTC‑AAGGTAAT‑GAGGTCGTCGACCAGTACGCGGATGTATGGCTCGGCGCGGCGCGGGGACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCttcttc                                                         <  8:324974/144‑1 (MQ=255)
                                                                                 gCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCg                                                  >  8:180941/1‑149 (MQ=255)
                                                                                  cTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCgg                                                 >  8:120143/1‑149 (MQ=255)
                                                                                     cTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCggtg                                               <  8:280716/148‑1 (MQ=255)
                                                                                      tGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCg                                                  >  5:100382/1‑144 (MQ=255)
                                                                                      tGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCg                                                  <  6:100382/144‑1 (MQ=255)
                                                                                       gAGTACCC‑AAGGTAAT‑CAGGTCGTCCCCCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGgtggtg                                            >  4:308366/1‑149 (MQ=255)
                                                                                            ccc‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGtt                                       >  8:89892/1‑149 (MQ=255)
                                                                                            ccc‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGtt                                       >  4:294259/1‑149 (MQ=255)
                                                                                            ccc‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGt                                        <  7:81206/148‑1 (MQ=255)
                                                                                             cc‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGtt                                       >  6:226344/1‑148 (MQ=255)
                                                                                             cc‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGtt                                       >  8:210751/1‑148 (MQ=255)
                                                                                             cc‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTg                                      <  1:217276/149‑1 (MQ=255)
                                                                                                aaGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGcttct                                                           <  1:87699/126‑1 (MQ=255)
                                                                                                aaGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGcttct                                                           >  2:87699/1‑126 (MQ=255)
                                                                                                aaGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGAt                                    >  8:166830/1‑149 (MQ=255)
                                                                                                  ggTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCt                                  <  7:210751/149‑1 (MQ=255)
                                                                                                     aaT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGcc                               >  4:8670/1‑149 (MQ=255)
                                                                                                     aaT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGcc                               >  2:100861/1‑149 (MQ=255)
                                                                                                           ggTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTCTCACGGCCCTGCGCCAGAAGAGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCCGCTTCTTCGTACCCGGAGGTGCCGTTCATCTGCCCGGc                           >  2:95210/1‑148 (MQ=255)
                                                                                                             tcgtcCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGCCCGGCgaa                        <  3:294259/149‑1 (MQ=255)
                                                                                                                  caccaGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGCCCGGCGAAGAAc                    <  1:100861/148‑1 (MQ=255)
                                                                                                                    ccaGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGCCCGGCg                          <  4:189559/140‑1 (MQ=255)
                                                                                                                    ccaGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGCCCGGCg                          >  3:189559/1‑140 (MQ=255)
                                                                                                                    ccaGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGCCCGGCGAAGAACAgg                 <  7:120143/149‑1 (MQ=255)
                                                                                                                           gCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGgc                                                                     >  5:110600/1‑90 (MQ=255)
                                                                                                                           gCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGgc                                                                     <  6:110600/90‑1 (MQ=255)
                                                                                                                                   aTGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGCCCGGCGAAGAACAGGCCGCCGATGACTTTc  >  2:288567/1‑149 (MQ=255)
                                                                                                                                    tGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGCCCGGCGAAGAACAGGCCGCCGATGACTTTc  <  3:8670/148‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                              
CCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCC‑AAGGTAAT‑CAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGCCCGGCGAAGAACAGGCCGCCGATGACTTTC  >  NC_002947/3025‑3302

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