Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,171,838 C→A 7.3% A174S (GCA→TCA)  glmS ← L‑glutamine/D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,171,8380CA 7.3% 203.5 / 3.6 64A174S (GCA→TCA) glmSL‑glutamine/D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base A (2/3);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAGAGACCTGATCACGGCGGATTTCAGCGATGTCGCCTTCTTCCAGGTACATGAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCG  >  NC_002947/6171689‑6171986
                                                                                                                                                     |                                                                                                                                                     
gATAGAGACCTGATCACGGCGGATTTCAGCGATGTCGCCTTCTTCCAGGTACATGAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGtgc                                                                                                                                                       <  1:1062318/150‑1 (MQ=255)
      gaCCTGATCACGGCGGATTTCAGCGATGTCGCCTTCTTCCAGGTACATGAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGccac                                                                                                                                                 >  2:1049533/1‑150 (MQ=255)
      gaCCTGATCACGGCGGATTTCAGCGATGTCGCCTTCTTCCAGGTACATGAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGccac                                                                                                                                                 <  1:971499/150‑1 (MQ=255)
         cTGATCACGGCGGATTTCAGCGATGTCGCCTTCTTCCAGGTACATGAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAg                                                                                                                                              >  1:1141017/1‑150 (MQ=255)
         cTGATCACGGCGGATTTCAGCGATGTCGCCTTCTTCCAGGTACATGAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAg                                                                                                                                              <  2:781236/150‑1 (MQ=255)
                  gcggATTTCAGCGATGTCGCCTTCTTCCAGGTACATGAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGATGCCACCAGGCGGTCCgg                                                                                                                                     >  2:904412/1‑150 (MQ=255)
                  gcggATTTCAGCGATGTCGCCTTCTTCCAGGTACATGAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGATGCCACCAGGCGGTCCgg                                                                                                                                     >  2:904559/1‑150 (MQ=255)
                      aTTTCAGCGATGTCGCCTTCTTCCAGGTACATGAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGgctgc                                                                                                                                 <  1:481880/150‑1 (MQ=255)
                        ttCAGCGATGTCGCCTTCTTCCAGGTACATGAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCtt                                                                                                                               <  2:690737/150‑1 (MQ=255)
                          cAGCGATGTCGCCTTCTTCCAGGTACATGAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTgg                                                                                                                             <  1:193492/150‑1 (MQ=255)
                                      cttcttCCAGGTACATGAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGAAGCCAGGAAAATCTCGCCAAATCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAgcg                                                                                                                 >  2:1136616/1‑150 (MQ=255)
                                                gTACATGAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTa                                                                                                       >  1:964337/1‑150 (MQ=255)
                                                 tACATGAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTaa                                                                                                      >  1:746487/1‑150 (MQ=255)
                                                      gAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGTCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGcc                                                                                                 >  2:1034851/1‑150 (MQ=255)
                                                      gAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGcc                                                                                                 >  1:759348/1‑150 (MQ=255)
                                                             tCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAgg                                                                                          >  1:258956/1‑150 (MQ=255)
                                                             tCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAgg                                                                                          >  1:507423/1‑150 (MQ=255)
                                                                   aCCGGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTc                                                                                    <  2:459248/150‑1 (MQ=255)
                                                                          gCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCAccgccg                                                                             >  2:539846/1‑150 (MQ=255)
                                                                                  aGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAgg                                                                    <  1:911958/150‑1 (MQ=255)
                                                                                        tCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCg                                                              >  1:382131/1‑150 (MQ=255)
                                                                                         cGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCgg                                                             >  2:962244/1‑150 (MQ=255)
                                                                                            aCGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCa                                                          >  1:48759/1‑150 (MQ=255)
                                                                                             cGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCAccgcc                                                                              <  1:933075/130‑1 (MQ=255)
                                                                                             cGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCAccgcc                                                                              >  2:933075/1‑130 (MQ=255)
                                                                                             cGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAg                                                         <  2:964337/150‑1 (MQ=255)
                                                                                              gCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAgg                                                        <  1:244972/150‑1 (MQ=255)
                                                                                                cAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAg                                                         >  1:165036/1‑147 (MQ=255)
                                                                                                cAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAg                                                         <  2:165036/147‑1 (MQ=255)
                                                                                                    aaGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCGCCTGTTGGTCCTGCTGCTTGGCACTGATCAGCGCCCGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCCGGTCTGGa                                                  <  2:455684/150‑1 (MQ=255)
                                                                                                       ttCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGtgaagc                                                       >  1:953947/1‑137 (MQ=255)
                                                                                                       ttCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGtgaagc                                                       <  2:953947/142‑6 (MQ=255)
                                                                                                       ttCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAg                                                                     <  1:1167103/129‑1 (MQ=255)
                                                                                                       ttCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAg                                                                     >  2:1167103/1‑129 (MQ=255)
                                                                                                          tcGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGctc                                            <  2:420249/150‑1 (MQ=255)
                                                                                                          ccGCCAAGGCCCAGGGCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGctc                                            <  1:372761/149‑1 (MQ=255)
                                                                                                            gCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGATCTGGAATGCTCtt                                          >  2:6247/1‑150 (MQ=255)
                                                                                                               aagtcccAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAg                                       <  1:1034851/146‑1 (MQ=255)
                                                                                                               aaGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAg                                       >  2:298070/1‑150 (MQ=255)
                                                                                                               aaGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAg                                       <  2:382131/150‑1 (MQ=255)
                                                                                                                  gCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGgt                                    >  1:614126/1‑150 (MQ=255)
                                                                                                                       ggCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtggt                               <  2:507423/150‑1 (MQ=255)
                                                                                                                        gCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtggtg                              >  1:581535/1‑150 (MQ=255)
                                                                                                                        gCCGATAACCAGCGGGCTACCACTGCGTGATGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtggtg                              <  1:904412/150‑1 (MQ=255)
                                                                                                                        gCCGATAACCAGCGGGCTACCACTGCGTGATGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtggtg                              <  1:904559/150‑1 (MQ=255)
                                                                                                                         ccGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGtggtgg                             <  2:1141017/150‑1 (MQ=255)
                                                                                                                          cGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGa                            <  2:746487/150‑1 (MQ=255)
                                                                                                                           gATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGAt                           >  1:391898/1‑150 (MQ=255)
                                                                                                                            aTAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATc                          <  1:962244/150‑1 (MQ=255)
                                                                                                                                ccAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAgg                                      <  1:55581/134‑1 (MQ=255)
                                                                                                                                ccAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAgg                                      >  2:55581/1‑134 (MQ=255)
                                                                                                                                  aGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTgg                    >  2:592980/1‑150 (MQ=255)
                                                                                                                                   gCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGa                   <  2:759348/150‑1 (MQ=255)
                                                                                                                                    cGGGCTACCACTGCGTGATGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGAc                  <  2:141623/150‑1 (MQ=255)
                                                                                                                                        cTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAggg                                                                   <  1:387362/97‑1 (MQ=255)
                                                                                                                                        cTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAggg                                                                   >  2:387362/1‑97 (MQ=255)
                                                                                                                                           ccACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACt           >  2:773867/1‑150 (MQ=255)
                                                                                                                                           ccACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACt           >  2:951304/1‑150 (MQ=255)
                                                                                                                                             aCTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAg                                                                                           <  2:570059/69‑1 (MQ=255)
                                                                                                                                             aCTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAg                                                                                           >  1:570059/1‑69 (MQ=255)
                                                                                                                                              cTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCg        >  2:781052/1‑150 (MQ=255)
                                                                                                                                              cTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCg        <  1:592980/150‑1 (MQ=255)
                                                                                                                                              cTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCg        <  1:1136616/150‑1 (MQ=255)
                                                                                                                                                  gTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGtgtg                                 >  1:407459/1‑121 (MQ=255)
                                                                                                                                                  gTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGtgtg                                 <  2:407459/121‑1 (MQ=255)
                                                                                                                                                    gctgcCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCg  <  2:391898/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATAGAGACCTGATCACGGCGGATTTCAGCGATGTCGCCTTCTTCCAGGTACATGAAGCGGTCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCG  >  NC_002947/6171689‑6171986

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