Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,753,353 G→A 6.6% R108H (CGC→CAC)  PP_1555 → putative Phage integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,753,3530GA 6.6% 178.2 / 2.1 60R108H (CGC→CAC) PP_1555putative Phage integrase
Reads supporting (aligned to +/- strand):  ref base G (29/27);  new base A (2/2);  total (31/29)
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.59e-01

GCTGCCAAGGTGGCCGCTGCGCTGCCGGATATCCGCCTGCGCCTGGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGCCGCTGCAAGAGTCGCTGTCGATCGACTATGTGCGCCTGGTC  >  NC_002947/1753205‑1753501
                                                                                                                                                    |                                                                                                                                                    
gcTGCCAAGGTGGCCGCTGCGCTGCCGGATATCCGCCTGCGCCTGGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAgcg                                                                                                                                                      >  5:47403/1‑149 (MQ=255)
    ccAAGGTGGCCGCTGCGCTGCCGGATATCCGCCTGCGCCTGGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCa                                                                                                                                                    <  6:47403/147‑1 (MQ=255)
     cAAGGTGGCCGCTGCGCTGCCGGATATCCGCCTGCGCCTGGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAAc                                                                                                                                                 >  7:27968/1‑149 (MQ=255)
         gTGGCCGCTGCGCTGCCGGATATCCGCCTGCGCCTGGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCAc                                                                                                                                             <  1:69792/149‑1 (MQ=255)
             ccGCTGCGCTGCCGGATATCCGCCTGCGCCTGGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCAcggcg                                                                                                                                         <  6:217748/149‑1 (MQ=255)
                 ggcgcaggCGGATATCCGGCAGCGCCTGGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGt                                                                                                                                     >  4:159919/2‑149 (MQ=255)
                   cgctgcCGGATATCCGCCTGCGCCTGGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCg                                                                                                                                   <  3:262362/149‑1 (MQ=255)
                    gctgcCGGATATCCGCCTGCGCCTGGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCgg                                                                                                                                  <  7:287670/149‑1 (MQ=255)
                         cGGATATCCGCCTGCGCCTGGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATc                                                                                                                             <  8:257720/149‑1 (MQ=255)
                               tCCGCCTGCGCCTGGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAagca                                                                                                                        <  8:281242/148‑1 (MQ=255)
                               tCCGCCTGCGCCTGGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAag                                                                                                                          >  7:281242/1‑146 (MQ=255)
                                 cgcctgcgcctgGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAgcagca                                                                                                                     <  5:169330/149‑1 (MQ=255)
                                  gcctgcgcctgGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCACAAACGCACGGCGGCGTCGGCGATCAAgcagca                                                                                                                     >  4:209858/1‑148 (MQ=255)
                                     tgcgcctgGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCACAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTg                                                                                                                 <  3:209858/149‑1 (MQ=255)
                                       cgcctgGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGgc                                                                                                               <  1:176354/149‑1 (MQ=255)
                                        gcctgGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCAc                                                                                                                                             >  2:158117/1‑118 (MQ=255)
                                            ggAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGgcgccgc                                                                                                          >  7:170979/1‑149 (MQ=255)
                                                    gCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTgg                                                                                                    >  8:121915/1‑147 (MQ=255)
                                                         ggATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCACAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAAc                                                                                             >  5:180842/1‑149 (MQ=255)
                                                           aTCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGc                                                                                           >  1:167112/1‑149 (MQ=255)
                                                               aGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCAtt                                                                                       >  5:238118/1‑149 (MQ=255)
                                                                gCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTg                                                                                      >  6:158223/1‑149 (MQ=255)
                                                                gCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTg                                                                                      <  8:27968/149‑1 (MQ=255)
                                                                  cTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCACAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGc                                                                                    <  6:180842/149‑1 (MQ=255)
                                                                       gAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGcc                                                                                   >  6:292679/1‑145 (MQ=255)
                                                                         gCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGcc                                                                                   >  2:263732/1‑143 (MQ=255)
                                                                             ggTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGa                                                                           <  3:159919/147‑1 (MQ=255)
                                                                              gTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGCCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCAGGGAGAAACGCCATTGGCCCAGATCGACa                                                                         >  4:270028/1‑148 (MQ=255)
                                                                                    aCTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAgcgc                                                                   <  2:167112/148‑1 (MQ=255)
                                                                                        gagAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGTCGATc                                                                                                                             >  7:199107/3‑86 (MQ=255)
                                                                                        gCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATc                                                                                                                             <  8:199107/86‑1 (MQ=255)
                                                                                        gCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCgctgc                                                                <  1:263732/147‑1 (MQ=255)
                                                                                           aGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGAc                                                           <  5:292679/149‑1 (MQ=255)
                                                                                           aGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGAc                                                           >  8:115309/1‑149 (MQ=255)
                                                                                            gTTGCTGGACTGGTTCGCTGAGCGCATGTGCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGAcc                                                          <  3:155964/149‑1 (MQ=255)
                                                                                                cTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGcc                                                                                   <  1:196561/120‑1 (MQ=255)
                                                                                                cTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGcc                                                                                   >  2:196561/1‑120 (MQ=255)
                                                                                                cTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGa                                                      <  7:115309/149‑1 (MQ=255)
                                                                                                    aCTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTg                                                  <  6:238118/149‑1 (MQ=255)
                                                                                                     cTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGa                                                 <  1:198786/149‑1 (MQ=255)
                                                                                                      tGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGAt                                                >  1:29202/1‑149 (MQ=255)
                                                                                                        gTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGccc                                                                                  >  2:116389/1‑113 (MQ=255)
                                                                                                        gTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGccc                                                                                  <  1:116389/113‑1 (MQ=255)
                                                                                                        gTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACcgc                                                        >  7:97445/1‑139 (MQ=255)
                                                                                                        gTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACcgc                                                        <  8:97445/139‑1 (MQ=255)
                                                                                                        gTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGAtgt                                              <  8:170979/149‑1 (MQ=255)
                                                                                                            gCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGcc                                          >  3:68942/1‑149 (MQ=255)
                                                                                                              tGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCTCGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCGGGGAGAACCGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGCCGc                                        >  3:47370/1‑149 (MQ=255)
                                                                                                                       gTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGCCGCTGCAAGAGt                               >  8:242818/1‑149 (MQ=255)
                                                                                                                          ccGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGCGCTGATGTGGCCGCTGCAAGAGTCGc                            >  3:25892/1‑149 (MQ=255)
                                                                                                                          ccGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGCCGCTGCAAGAGTCGc                            >  5:178345/1‑149 (MQ=255)
                                                                                                                            gTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGCCGCTGCAAGAGTCGCTg                          <  3:270028/149‑1 (MQ=255)
                                                                                                                                   cAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGCCGCTGCAaga                                 >  2:137888/1‑135 (MQ=255)
                                                                                                                                   cAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGCCGCTGCAaga                                 <  1:137888/135‑1 (MQ=255)
                                                                                                                                        tctcGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGCCGCTGCAAGAGTCGCTGTCGATCGACTAt              >  1:199107/1‑149 (MQ=255)
                                                                                                                                        tctcGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGCCGCTGCAAGAGTCGCTGTCGATCGACTAt              >  7:11428/1‑149 (MQ=255)
                                                                                                                                        tctcGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGCCGCTGCAAGAGTCGCTGTCGATCGACTAt              >  2:68735/1‑149 (MQ=255)
                                                                                                                                         ctcGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGCCGCTGCAAGAGTCGCTGTCGATCGACTAtg             >  1:82201/1‑149 (MQ=255)
                                                                                                                                            ggCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGCCGCTGCAAGAGTCGCTGTCGATCGACTATGTgc          <  2:82201/149‑1 (MQ=255)
                                                                                                                                                aaGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGAcc                                                          >  8:99714/1‑97 (MQ=255)
                                                                                                                                                aaGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGAcc                                                          <  7:99714/97‑1 (MQ=255)
                                                                                                                                                  gcgcAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGCCGCTGCAAGAGTCGCTGTCGATCGACTATGTGCGCCTgg    >  1:12940/1‑149 (MQ=255)
                                                                                                                                                    gcAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGCCGCTGCAAGAGTCGCTGTCGATCGACTATGTGCGCCTGGTc  >  3:188227/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTGCCAAGGTGGCCGCTGCGCTGCCGGATATCCGCCTGCGCCTGGAAGCCGGCGAGGGATCGAGCCTGTCGAGCTGGGTGCTGACTGGCGAGTTGCTGGACTGGTTCGCTGAGCGCATGTCCCGTGACCGCAGCCTCTCGGCCAAGCGCAAACGCACGGCGGCGTCGGCGATCAAGCAGCATCTGGCGCCGCGCCTGGGAGAAACGCCATTGGCCCAGATCGACAAAGCGCTGCTCGACCGCGAGCTGATGTGGCCGCTGCAAGAGTCGCTGTCGATCGACTATGTGCGCCTGGTC  >  NC_002947/1753205‑1753501

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