Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,043,722 T→A 7.8% intergenic (+18/+16) gmk → / ← PP_5297 deoxyguanylate kinase/guanylate kinase/amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,043,7220TA 7.8% 396.4 / 18.4 129intergenic (+18/+16)gmk/PP_5297deoxyguanylate kinase/guanylate kinase/amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (55/63);  new base A (5/5);  total (61/68)
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.26e-01

Only 100 of 132 total aligned reads displayed.

GACTACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAACACCAGC  >  NC_002947/6043573‑6043870
                                                                                                                                                     |                                                                                                                                                    
gACTACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAt                                                                                                                                                      >  1:1669404/1‑150 (MQ=255)
    aCGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCa                                                                                                                                                  <  1:971034/150‑1 (MQ=255)
            gtcaATGATGATCTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAAt                                                                                                                                          >  2:1652524/2‑150 (MQ=255)
            atcaATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAAt                                                                                                                                          <  1:1305850/150‑1 (MQ=255)
            agcaATGATGATTTCGACACGGCGCTAGCAGACTTGAGAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGATGGAGATCGATGGGGCTGCATTGCAGCCCCAAt                                                                                                                                          <  1:429506/148‑1 (MQ=255)
             tcaATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATg                                                                                                                                         >  2:1375848/1‑150 (MQ=255)
              caATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATgg                                                                                                                                        >  2:2253822/1‑150 (MQ=255)
               aaTGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGt                                                                                                                                       <  2:1628098/150‑1 (MQ=255)
                atgatgATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCCTTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                      <  2:1898663/150‑1 (MQ=255)
                 tgatgaTTTCGACACGGCGCTGGCAGCCTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCAGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtat                                                                                                                                     >  2:2053722/1‑148 (MQ=255)
                    tgatTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGt                                                                                                                                  >  2:41827/1‑150 (MQ=255)
                       tttCGAGACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCTGCTGCTGGTGAAGAAGCAGGAGCAGCGTCATGCGTTCTTGCTCAAGCAAATGATCGGCAGACGCTCGATGGGGCCGTATTGCAGTCCCACGGGTGTAgggcg                                                                                                                               >  2:394679/1‑146 (MQ=255)
                        ttCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTg                                                                                                                              <  1:2211958/150‑1 (MQ=255)
                              acGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCgg                                                                                                                        >  1:1758594/1‑150 (MQ=255)
                              acGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCgg                                                                                                                        >  1:48480/1‑150 (MQ=255)
                                 gcgcTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCaa                                                                                                                     <  2:312222/150‑1 (MQ=255)
                                  cgcTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAg                                                                                                                    <  1:24907/150‑1 (MQ=255)
                                     tGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTa                                                                                                                 <  1:1771791/150‑1 (MQ=255)
                                      ggCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAg                                                                                                                >  2:1660672/1‑150 (MQ=255)
                                       gCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATg                                                                                                                                         <  2:1147332/124‑1 (MQ=255)
                                       gCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATg                                                                                                                                         >  1:1147332/1‑124 (MQ=255)
                                          gACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTa                                                                                                                                    >  2:772042/1‑126 (MQ=255)
                                          gACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTa                                                                                                                                    <  1:772042/126‑1 (MQ=255)
                                             ttGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAAt                                                                                                                                          >  1:1630873/1‑117 (MQ=255)
                                             ttGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAg                                                                                                                           <  2:1596450/132‑1 (MQ=255)
                                             ttGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAg                                                                                                                           >  1:1596450/1‑132 (MQ=255)
                                             ttGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAATGCAGCCCCa                                                                                                                                            <  1:322892/115‑1 (MQ=255)
                                             ttGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAATGCAGCCCCa                                                                                                                                            >  2:322892/1‑115 (MQ=255)
                                                 aaGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGc                                                                                                     <  2:48480/150‑1 (MQ=255)
                                                  aGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAg                                                                                                                    >  1:2116332/1‑134 (MQ=255)
                                                  aGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCctgcgtcagca                                                                                                    >  2:1860085/1‑143 (MQ=255)
                                                  aGCAGTGTTACGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAATGCAGCCCCat                                                                                                                                           <  2:1308760/111‑2 (MQ=255)
                                                   gCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                      <  2:403729/115‑1 (MQ=255)
                                                   gCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                      >  1:403729/1‑115 (MQ=255)
                                                    cAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAgg                                                                                                  <  2:1758594/150‑1 (MQ=255)
                                                    cAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTCACGATCGATGGGGCTGCATTGCAGCCACAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAgg                                                                                                  <  1:1375848/150‑1 (MQ=255)
                                                       tgtTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTa                                                                                               <  1:1652524/150‑1 (MQ=255)
                                                        gtTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCgg                                                                                                                        >  1:1746868/1‑124 (MQ=255)
                                                         ttCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCAATGGGGCTGCATTGCAGccc                                                                                                                                              >  2:626970/1‑101 (MQ=255)
                                                         ttCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCAATGGGGCTGCATTGCAGccc                                                                                                                                              <  1:626970/101‑1 (MQ=255)
                                                            cGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAgg                                                                                                  <  2:460307/142‑1 (MQ=255)
                                                            cGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAgg                                                                                                  >  1:460307/1‑142 (MQ=255)
                                                             gCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCg                                                                                         >  2:164465/1‑150 (MQ=255)
                                                             gCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCg                                                                                         >  2:1042163/1‑150 (MQ=255)
                                                                  aaTCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGaa                                                                                    <  2:1669404/150‑1 (MQ=255)
                                                                   aTCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCAGCAAGGTAGCCTGCTGACGCAGGGTAACCCCGAGGAAg                                                                                   <  2:637150/150‑1 (MQ=255)
                                                                     cGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTa                                                                                 <  1:1628070/150‑1 (MQ=255)
                                                                       gctgctgctGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAt                                                                                                                                                      >  1:115452/1‑79 (MQ=255)
                                                                       gctgctgctGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAt                                                                                                                                                      <  2:115452/79‑1 (MQ=255)
                                                                         tgctggtgAATAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGcc                                                                             >  2:2211782/1‑150 (MQ=255)
                                                                         tgctgctgAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGcc                                                                             >  2:1644686/1‑150 (MQ=255)
                                                                         tgctgctgAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGcc                                                                             >  1:523906/1‑150 (MQ=255)
                                                                           ctgctgAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCgg                                                                           <  1:1042163/150‑1 (MQ=255)
                                                                                gaagaaGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGcgcac                                                                      >  1:408824/1‑150 (MQ=255)
                                                                                gaagaaGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGcgcac                                                                      <  1:1822449/150‑1 (MQ=255)
                                                                                      gcagcagcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAg                                                                                                                    >  1:1020625/1‑98 (MQ=255)
                                                                                          cagcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGa                                                            <  1:2211782/150‑1 (MQ=255)
                                                                                            gcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                      >  2:2187926/1‑74 (MQ=255)
                                                                                            gcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                      <  1:2187926/74‑1 (MQ=255)
                                                                                            gcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAg                                                                                       >  2:512623/1‑121 (MQ=255)
                                                                                            gcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAg                                                                                       <  1:512623/121‑1 (MQ=255)
                                                                                                     tGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCAtcatca                                                 <  1:33905/150‑1 (MQ=255)
                                                                                                     tGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCAtcatca                                                 >  2:212335/1‑150 (MQ=255)
                                                                                                        ggCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGc                                              <  2:523906/150‑1 (MQ=255)
                                                                                                         gCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGcc                                             <  2:453313/150‑1 (MQ=255)
                                                                                                          cATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCg                                            <  1:910779/150‑1 (MQ=255)
                                                                                                           aTTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGc                                                                                                               >  2:1302268/1‑82 (MQ=255)
                                                                                                           aTTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGc                                                                                                               <  1:1302268/82‑1 (MQ=255)
                                                                                                           aTTGCTCAAGCAAATCAGAAGCTGAAGAGCGATGGGGCTGCATGTCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTAACGCACGGTAGCCCCGACGTAGTACGCCTGCGCACTCATGCGCGACAGCATCATCAGGCCGa                                           >  2:1728108/1‑150 (MQ=255)
                                                                                                              gCTCAAGCAAATGATCGGCTGAAGATCGATGGAGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGc                                        <  2:2238521/150‑1 (MQ=255)
                                                                                                               cTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAg                                                                                   <  1:464620/106‑1 (MQ=255)
                                                                                                                tCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCa                                      >  1:2397627/1‑150 (MQ=255)
                                                                                                                tCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCa                                      >  2:588291/1‑150 (MQ=255)
                                                                                                                tCAAGCAAATGATCGGCCGAAGCTCGATGGGGCTGCATTGCAGCCCCAATGGTGTAGTCCTGCAGCAGCAAGGTAGCCTGCTGACGCAGGGTAGCCGCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCa                                      >  2:735327/1‑150 (MQ=255)
                                                                                                                  aaGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAg                                    >  1:1293077/1‑150 (MQ=255)
                                                                                                                     cAAATGATCGGCTGCAGATTGATGCGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGTCAAGGTAACCTGCTCACGCAGGGTAGCCCCTAGGAAGTTCGCCTCCTCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGGAc                                 <  2:762400/150‑1 (MQ=255)
                                                                                                                     cAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAAt                                                                                                                                          >  2:2147597/1‑45 (MQ=255)
                                                                                                                     cAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAAt                                                                                                                                          <  1:2147597/45‑1 (MQ=255)
                                                                                                                       aaTGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCaa                                     >  1:2011892/1‑144 (MQ=255)
                                                                                                                       aaTGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCaa                                     <  2:2011892/144‑1 (MQ=255)
                                                                                                                        aTGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCga                              <  2:1799513/150‑1 (MQ=255)
                                                                                                                        aTGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCga                              <  1:2253822/150‑1 (MQ=255)
                                                                                                                          gATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGAt                                          <  2:2243582/136‑1 (MQ=255)
                                                                                                                           agcGGCTGAAAATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGGCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAgat                           <  1:394679/148‑1 (MQ=255)
                                                                                                                              ggCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTACAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAgatgat                        >  2:378473/1‑150 (MQ=255)
                                                                                                                               gCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATc                       >  1:1654929/1‑150 (MQ=255)
                                                                                                                                  gAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCgag                             <  2:1197600/141‑1 (MQ=255)
                                                                                                                                  gAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGAGTCAGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATcgcg                    >  1:150133/1‑150 (MQ=255)
                                                                                                                                     gATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGa                                                                                        <  2:285856/79‑1 (MQ=255)
                                                                                                                                     gATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGa                                                                                        >  1:285856/1‑79 (MQ=255)
                                                                                                                                         attGGGGCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCAtcatca                                                 <  1:1948688/112‑1 (MQ=255)
                                                                                                                                         attGGGGCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCAtcatca                                                 >  2:1948688/3‑114 (MQ=255)
                                                                                                                                         attGGGGCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCa             >  1:1491108/3‑150 (MQ=255)
                                                                                                                                              ggCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAac        >  1:2095257/1‑150 (MQ=255)
                                                                                                                                              ggCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACCCAGGGGAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAac        >  1:903196/1‑150 (MQ=255)
                                                                                                                                              ggCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCTTGCGCGACAGCATCATCAGGCCGATGCCAAGCACCga                              >  2:2060595/1‑128 (MQ=255)
                                                                                                                                              ggCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCTTGCGCGACAGCATCATCAGGCCGATGCCAAGCACCga                              <  1:2060595/128‑1 (MQ=255)
                                                                                                                                               gCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAaca       >  1:831049/1‑150 (MQ=255)
                                                                                                                                               ccTGCATTGCAGCCCCAATGGTTGAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAaca       >  2:151361/2‑150 (MQ=255)
                                                                                                                                                    aTTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAACAccagc  <  2:2095257/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACTACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAACACCAGC  >  NC_002947/6043573‑6043870

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