Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,496,815 T→G 6.2% intergenic (+47/+25) mdeA → / ← PP_1309 methionine gamma‑lyase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,496,8150TG 6.2% 183.0 / 2.0 59intergenic (+47/+25)mdeA/PP_1309methionine gamma‑lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base G (2/2);  total (30/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 = 1.00e+00

ATCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTT  >  NC_002947/1496668‑1496958
                                                                                                                                                      |                                                                                                                                               
aTCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCtt                                                                                                                                                <  2:354925/149‑1 (MQ=255)
  cACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGACCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGc                                                                                                                                              >  7:92229/1‑149 (MQ=255)
   aCGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCg                                                                                                                                             <  3:321922/149‑1 (MQ=255)
              gAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGa                                                                                                                                            >  1:270097/1‑139 (MQ=255)
              gAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGa                                                                                                                                            <  2:270097/139‑1 (MQ=255)
                 gggTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCGCAAGAGGCCttttttgtg                                                                                                                               <  3:56510/149‑7 (MQ=255)
                  ggTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACAt                                                                                                                              <  6:72432/149‑1 (MQ=255)
                   gTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCGATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCGTTGTcacagt                                                                                                                             >  6:191338/1‑147 (MQ=255)
                    ttGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTc                                                                                                                            <  6:87262/149‑1 (MQ=255)
                          cGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCa                                                                                                                      <  5:348939/149‑1 (MQ=255)
                               gtaggTGGGGCTGGAGGATGTGGTTGTCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGACCCACAAAAAAGGCC‑‑TCTTGCGAg                                                                                                                                           <  7:223458/121‑1 (MQ=255)
                               gTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAg                                                                                                                                           >  8:223458/1‑123 (MQ=255)
                                tCGGTGGGACTGAAGAATGTGAAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTGGCGCGGCGTTTGTCACGTTGAGATCAACGcgg                                                                                                                >  8:53897/1‑147 (MQ=255)
                                            gaggaTGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCt                                                                                                    >  7:130952/1‑149 (MQ=255)
                                              ggaTGTGGAGGTCCTGCTGGCCGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCttt                                                                                                                                     <  2:60732/114‑1 (MQ=255)
                                              ggaTGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCttt                                                                                                  >  8:299230/1‑149 (MQ=255)
                                                 tgtgGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGAt                                                                                                                        <  6:359166/124‑1 (MQ=255)
                                                 tgtgGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGAt                                                                                                                        >  5:359166/1‑124 (MQ=255)
                                                   tgGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGACCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTtcgtc                                                                                              <  8:92229/148‑1 (MQ=255)
                                                     gaggaCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGa                                                                                            <  1:164557/148‑1 (MQ=255)
                                                     gaggaCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAg                                                                                           >  7:219499/1‑149 (MQ=255)
                                                            tgctgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGc                                                                                      >  7:156126/1‑147 (MQ=255)
                                                                ggCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGc                                                                                 >  2:352982/1‑148 (MQ=255)
                                                                 gCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGc                                                                                 >  8:41547/1‑147 (MQ=255)
                                                                   aGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGt                                                                             >  4:237269/1‑149 (MQ=255)
                                                                      tGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAg                                                                          <  8:156126/149‑1 (MQ=255)
                                                                          gAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGgcg                                                                       <  4:110146/148‑1 (MQ=255)
                                                                              aGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCAGGCTAGCGAGCCCTTTGTCACATACAGATAAACGCTAGAGCCCGTAGCTTTCGACGAGCATGCCCGGCGAGTTAGGCGCtt                                                                    >  1:111973/1‑149 (MQ=255)
                                                                                gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTa                                                                                                       <  6:233806/110‑1 (MQ=255)
                                                                                gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTa                                                                                                       >  5:233806/1‑110 (MQ=255)
                                                                                gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTATGCGCTTTctg                                                                >  4:252410/1‑147 (MQ=255)
                                                                                  aTTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTg                                                              >  4:128624/1‑149 (MQ=255)
                                                                                             gcttTTCGGTTGCCCGCGAACCCGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTTCGATGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATTCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgcgc                                                   <  7:53897/149‑1 (MQ=255)
                                                                                              ctgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgcgc                                                   <  7:299230/148‑1 (MQ=255)
                                                                                                 ttAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCa                                               <  8:219499/149‑1 (MQ=255)
                                                                                                    ggTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACtt                                            >  3:146491/1‑149 (MQ=255)
                                                                                                                    gCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAgg                            <  7:41547/149‑1 (MQ=255)
                                                                                                                     ccATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGt                           >  2:344025/1‑149 (MQ=255)
                                                                                                                      cATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTg                          >  5:7658/1‑149 (MQ=255)
                                                                                                                       aTTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgg                         <  6:7658/149‑1 (MQ=255)
                                                                                                                       aTTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgg                         >  1:46304/1‑149 (MQ=255)
                                                                                                                       aTTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgg                         <  1:352982/149‑1 (MQ=255)
                                                                                                                        ttGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTggg                        <  7:73874/149‑1 (MQ=255)
                                                                                                                           aaTCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGt                      >  4:26668/1‑148 (MQ=255)
                                                                                                                            aTCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTAAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCa                    >  6:217315/1‑149 (MQ=255)
                                                                                                                               cAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAgg                  <  5:217315/148‑1 (MQ=255)
                                                                                                                                 gCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCgg               <  3:237269/149‑1 (MQ=255)
                                                                                                                                    ccACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAg                   <  1:83638/143‑1 (MQ=255)
                                                                                                                                    ccACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAg                   >  2:83638/1‑143 (MQ=255)
                                                                                                                                     cacaAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGa                                                                                                              >  7:230441/1‑51 (MQ=255)
                                                                                                                                     cacaAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGa                                                                                                              <  8:230441/51‑1 (MQ=255)
                                                                                                                                     cacaAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCa            >  4:117735/1‑149 (MQ=255)
                                                                                                                                        gacaaaGGCC‑‑TCGCAAGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTcg                                          <  2:313505/113‑1 (MQ=255)
                                                                                                                                        gacaaaGGCC‑‑TCGCAAGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTcg                                          >  1:313505/4‑116 (MQ=255)
                                                                                                                                        gacaaaGGCC‑‑TCGCAAGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGAGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGGGTTTAGGTGGGGTCAGGCGGCCCcgc         >  8:258622/4‑149 (MQ=255)
                                                                                                                                          aaaaGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGc                                                <  4:269868/108‑1 (MQ=255)
                                                                                                                                          aaaaGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGc                                                >  3:269868/1‑108 (MQ=255)
                                                                                                                                            aaGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGGGGGGTCAGGCGGTCACGCGGGc     <  3:175035/149‑1 (MQ=255)
                                                                                                                                             aGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGcc    <  5:191338/149‑1 (MQ=255)
                                                                                                                                               gCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCtt  >  4:227591/1‑149 (MQ=255)
                                                                                                                                                 c‑‑tctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTtcg                                                                                                >  1:72806/1‑53 (MQ=255)
                                                                                                                                                 c‑‑tctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTtcg                                                                                                <  2:72806/53‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
ATCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCC‑‑TCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTT  >  NC_002947/1496668‑1496958

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