Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,703,248 A→G 18.9% intergenic (+97/‑70) lysS → / → PP_1497 lysine‑‑tRNA ligase/TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,703,2480AG18.9% 104.1 / 18.7 58intergenic (+97/‑70)lysS/PP_1497lysine‑‑tRNA ligase/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (26/17);  new base G (6/4);  total (32/26)
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.87e-01

GCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCC  >  NC_002947/1703100‑1703393
                                                                                                                                                      |                                                                                                                                                 
gCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAGaaa                                                                                                                                                   <  4:344197/149‑1 (MQ=255)
 ccGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAg                                                                                                                                                  <  2:85470/149‑1 (MQ=255)
  cGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAgg                                                                                                                                                 >  3:47246/1‑149 (MQ=255)
       ctCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGAGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTg                                                                                                                                            <  4:234369/148‑1 (MQ=255)
       aTCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTCGTCAGTGCTGGCCCCTTCGCGGGTAAACCCGCTCCTACAAGGCCtttg                                                                                                                                           <  5:137866/149‑3 (MQ=38)
        tCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGt                                                                                                                                           <  3:110128/149‑1 (MQ=255)
        gccGCGACGTGATCCTGAGCGCGCACATGCGTCCGCAGGACTGAGTGAGCTGCACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGATGCTT‑‑TATGTTGTCAGTACCGGCCTGTTCGCGGGTATAGCCGCTCCTACAAGGCCtttgt                                                                                                                                          <  4:170807/148‑4 (MQ=255)
          cgcgACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTGCTGGCCCCTTCGCGGGTAAACCCGCTCCTACAAGGCCTTTGTAg                                                                                                                                         <  7:203881/149‑1 (MQ=38)
          cgcgACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGCCTTTGTa                                                                                                                                          <  5:86254/148‑1 (MQ=255)
          cgcgACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGCCTTTGTAg                                                                                                                                         <  6:240309/149‑1 (MQ=255)
           gcgACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAAACCGCTCCGACAAAGGGCTTGTAgg                                                                                                                                        >  7:320444/1‑149 (MQ=255)
           gcgACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAACCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAg                                                                                                                                         >  7:371843/1‑148 (MQ=255)
            cgACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCCCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGa                                                                                                                                       <  6:318630/149‑1 (MQ=255)
                gTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCgg                                                                                                                                   <  3:174884/149‑1 (MQ=255)
                         ttCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTAccc                                                                                                                           <  5:149603/148‑1 (MQ=255)
                         ttCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCcg                                                                                                                          <  7:20794/149‑1 (MQ=255)
                            ccGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCcgcg                                                                                                                        >  8:291590/1‑148 (MQ=255)
                              gCACATGCGTCCGCAGGCCTGAGTGAGCGGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGACGCTT‑‑TATGTGGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAg                                                                                                                     <  5:273197/149‑1 (MQ=255)
                                      gTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTGGTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCCCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGc                                                                                                                 >  1:314539/1‑147 (MQ=255)
                                      gTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAg                                                                                                              >  1:220018/1‑148 (MQ=255)
                                                       agCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACAcccc                                                                                             >  5:387472/1‑148 (MQ=255)
                                                        gCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCAAAAAAGGCCGTGTAGGAGCGGGTTTACCCGCGAAGGGGGCAGCACTGACGAAAcccccc                                                                                             >  8:206639/1‑147 (MQ=38)
                                                          tGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTAccc                                                                                                                           <  6:84399/115‑1 (MQ=255)
                                                          tGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTAccc                                                                                                                           >  5:84399/1‑115 (MQ=255)
                                                                 gccgccTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGGCTTGTAGGCGCGGGTTTTCCCGCGAAGGGGGCCGCACGGACGGACCACCCCCATGAGGGAtg                                                                                  >  8:232850/1‑148 (MQ=37)
                                                                    gccTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGt                                                                                                                                 <  3:280421/99‑1 (MQ=255)
                                                                    gccTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGt                                                                                                                                 >  4:280421/1‑99 (MQ=255)
                                                                         tCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGt                                                                            <  6:387472/147‑1 (MQ=255)
                                                                                ggcTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGc                                                                   <  7:291590/149‑1 (MQ=255)
                                                                                   tttttCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCAAGAGGTATACCCCGTGACAccccccct                                                                >  3:158561/1‑147 (MQ=255)
                                                                                           gCCTTATGCTT‑‑TTTGTTGTCCTTACCGGCCTTTTCGCGGGTAAACCCGCTCATAAAATGGCTTTTTTGGTGGGGTTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTg                                                           <  4:158561/146‑1 (MQ=255)
                                                                                           gCCTGAAGGTT‑‑TATGTTGTCAGTCCCGGCCTCTTCGCGGTTAAACCCGCTCCTACAATGGCCTTGTTGTTGCGGTTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGAc                                                         <  2:314539/148‑1 (MQ=255)
                                                                                               gAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTTGGAGCGGGTTTACCAGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGGATCCCCCGTGCCACCCGCCATGCCTGACCACgg                                                    >  1:48222/1‑149 (MQ=38)
                                                                                                        tATGTTGTCAGTCCCGGGGTGTTCGCGGGTAAACCCGCTCCTACAAAGGCGTTGTTGGCGCGGGTTAACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCGCCGTGACACCCGCCATGCCTGACCAAGGcgccgc                                              <  2:48222/148‑1 (MQ=37)
                                                                                                         aTGTTGTCAGTACCGGCCTCTTCGCGGTTAAACCCGCTCCTACAAAGGCCTTGTCGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCgccgcc                                             <  8:371843/148‑1 (MQ=255)
                                                                                                              gtCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATACCCCCTGACACCCGCCATGCCGGACCAAGGCGCCGCCGGCATc                                       >  6:181209/1‑149 (MQ=255)
                                                                                                               tCAGTCCCGGCCTCTTCGCGGTAAAACCCGCTCCTACAACGGCCTTGTTGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCg                                      <  2:220018/149‑1 (MQ=38)
                                                                                                                cAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGc                                     >  4:271547/1‑149 (MQ=255)
                                                                                                                      cGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCGTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCGCACTGCCGAAACACCCCCATGAGGTTTGCCCCGGGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCAccgc                               >  1:68086/1‑149 (MQ=37)
                                                                                                                        gCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAACCACCCCCATGAGGTATCCCCCGTGACACCCGCCAGGCCTGACCAAGGCGCCGCCGGCATCGCCAccgccg                             >  5:397382/1‑149 (MQ=255)
                                                                                                                          ctctTCGCGGGTAAACCCCCACCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCAAGAGGTATCACCCCTGACACCCGCCATGCCAGGCCAACGCGCCGCCGGCAACGCCAccccccgc                           >  4:54028/1‑146 (MQ=25)
                                                                                                                             ttCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAACCACCCCCATGAGGTATCCCCCGTGACCCCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGcc                        >  3:287191/1‑149 (MQ=255)
                                                                                                                             ttCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGcc                        >  6:342135/1‑149 (MQ=255)
                                                                                                                                 cgGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGa                      >  3:301583/1‑147 (MQ=255)
                                                                                                                                  gggTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGCCCAAGGCGCCGCCGGCATCGCCCCCGCCGGCGCCCAAAg                   >  3:243729/1‑149 (MQ=255)
                                                                                                                                       aaCCGGCTCCTACAATGGCCTTGTCGGAGGGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAATCGTGc              <  7:355516/149‑1 (MQ=255)
                                                                                                                                        aCCCGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAACACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGc              >  7:307643/1‑148 (MQ=255)
                                                                                                                                         cccGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGcc                        >  6:32109/1‑137 (MQ=255)
                                                                                                                                         cccGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGcc                        <  5:32109/137‑1 (MQ=255)
                                                                                                                                          ccccTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCACCGTGACACCCGCCATGCCTGACAAAGGCGCCGCCGGCATCGCCACCGCCGTCCCCGAAACCGTGCaa            >  2:345878/4‑148 (MQ=37)
                                                                                                                                              tCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATAcc        >  5:129974/1‑148 (MQ=255)
                                                                                                                                              tCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATAcc        >  3:184728/1‑148 (MQ=255)
                                                                                                                                               ccTACAAAGGCCTTGTAGGCGGGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCaa      <  2:68086/149‑1 (MQ=255)
                                                                                                                                                cTACATTGGCCTTTTAGTTGCGGTTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGGCGTCGCCGAAAGCGTGCAATACCaa      <  5:334780/148‑1 (MQ=255)
                                                                                                                                                cTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAg     >  6:114706/1‑149 (MQ=255)
                                                                                                                                                 tACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAgg    >  2:329961/1‑149 (MQ=255)
                                                                                                                                                  aCAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCACCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACAAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCCAGGc   >  3:354532/1‑149 (MQ=255)
                                                                                                                                                   cAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGcc  >  3:137052/1‑149 (MQ=255)
                                                                                                                                                   cAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCAATGCCTGAACAATGATCCTCCGGAATCGCCACCGCCGTCGCAGCAAGCGTTCAATACCAATGcc  >  5:220690/1‑149 (MQ=38)
                                                                                                                                                      |                                                                                                                                                 
GCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTT‑‑TATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCC  >  NC_002947/1703100‑1703393

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: