Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 780,563 G→T 8.1% intergenic (‑90/+32) PP_0670 ← / ← glyA‑II bile acid/Na+ symporter family transporter/serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947780,5630GT 8.1% 180.1 / 5.5 62intergenic (‑90/+32)PP_0670/glyA‑IIbile acid/Na+ symporter family transporter/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (26/31);  new base T (2/3);  total (28/34)
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.22e-01

GGCCAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGCCG  >  NC_002947/780414‑780703
                                                                                                                                                     |                                                                                                                                            
ggCCAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTgg                                                                                                                                              >  1:80133/1‑150 (MQ=255)
           cggcgAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCtt                                                                                                                                   <  1:1425232/150‑1 (MQ=255)
             gcgAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTc                                                                                                                                 <  1:170463/150‑1 (MQ=255)
               gAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCAc                                                                                                                               <  1:1400797/150‑1 (MQ=255)
                aaGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCCCAGTGCCGGGCTCTTTcgtt                                                                                                                              <  1:488375/150‑4 (MQ=255)
                   tAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTa                                                                                                                           <  1:430299/150‑1 (MQ=255)
                    aGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAg                                                                                                                          <  2:284090/150‑1 (MQ=255)
                      aTGGGGATCTGCCGTGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGTTGCGCCTGTATTAAACGAAAGAGCCCCGCACTGGGCCGGGCTCTTTTACTTTAGGc                                                                                                                        <  2:128168/150‑1 (MQ=255)
                       tGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCa                                                                                                                       >  1:1056130/1‑150 (MQ=255)
                        ggggATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAgg                                                                                                                         <  2:1158786/147‑1 (MQ=255)
                        ggggATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAgg                                                                                                                         >  1:1158786/1‑147 (MQ=255)
                            aTCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAgcgc                                                                                                                  >  2:1135795/1‑150 (MQ=255)
                               tGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCaa                                                                                                               >  2:361453/1‑150 (MQ=255)
                                   gggTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAg                                                                                                           >  1:133015/1‑150 (MQ=255)
                                         cGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCg                                                                                                     <  2:80133/150‑1 (MQ=255)
                                              tGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAAc                                                                                                <  2:1381933/150‑1 (MQ=255)
                                               gCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAAcc                                                                                               >  1:190607/1‑150 (MQ=255)
                                                ctctctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCg                                                                                              >  2:1268468/1‑150 (MQ=255)
                                                  ctctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCCCAGTGCCGGGCTCtt                                                                                                                                   <  2:8451/111‑1 (MQ=39)
                                                  ctctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCCCAGTGCCGGGCTCtt                                                                                                                                   >  1:8451/1‑111 (MQ=39)
                                                   tctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCCCAGTGCCGGGCTCTTTc                                                                                                                                 >  2:1189897/1‑112 (MQ=255)
                                                   tctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCCCAGTGCCGGGCTCTTTc                                                                                                                                 <  1:1189897/112‑1 (MQ=39)
                                                   tctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGAAGGCGCCAATCAGTTGCCGTAAACCGGCa                                                                                           >  1:1071089/1‑150 (MQ=255)
                                                   tctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGGCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCa                                                                                           <  2:157480/150‑1 (MQ=255)
                                                    ctGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAg                                                                                          <  1:417625/150‑1 (MQ=255)
                                                        tCAAGGTGCGGCCTGTGGAAAGTGGGGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGcttc                                                                                      <  2:183006/150‑1 (MQ=255)
                                                        tCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGctc                                                                                                                                     >  2:1418800/1‑103 (MQ=255)
                                                        tCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGctc                                                                                                                                     <  1:1418800/103‑1 (MQ=255)
                                                                cGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGAt                                                                              <  2:190607/150‑1 (MQ=255)
                                                                   ccTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGc                                                                           <  1:215619/150‑1 (MQ=255)
                                                                     tgtgCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCt                                                                         <  1:635200/150‑1 (MQ=255)
                                                                        gCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAg                                                                                                                          <  1:278559/98‑1 (MQ=255)
                                                                        gCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAg                                                                                                                          >  2:278559/1‑98 (MQ=255)
                                                                        gCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGa                                                                      <  2:1071089/150‑1 (MQ=255)
                                                                             gTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGcc                                                                                                      <  1:225250/113‑1 (MQ=255)
                                                                             gTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGcc                                                                                                      >  2:225250/1‑113 (MQ=255)
                                                                                 aGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGAGTTGACCTTCTCac                                                             >  2:583775/1‑150 (MQ=255)
                                                                                  ggTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAAGCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCacg                                                            >  1:1379110/1‑150 (MQ=255)
                                                                                              gcAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGAt                                                <  2:169772/150‑1 (MQ=255)
                                                                                                  tCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCaa                                                                                                               <  1:343934/83‑1 (MQ=255)
                                                                                                  tCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCaa                                                                                                               >  2:343934/1‑83 (MQ=255)
                                                                                                    tttGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGACAATGAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTTTCACGGACGGCGTCGATCACCGc                                          >  1:450128/1‑150 (MQ=255)
                                                                                                        tACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCa                                                                                           <  2:1308593/97‑1 (MQ=255)
                                                                                                        tACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCa                                                                                           >  1:1308593/1‑97 (MQ=255)
                                                                                                           gAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCgttg                                   >  1:1392935/1‑150 (MQ=255)
                                                                                                             ggggCCGTGCGCCTGTCTTACACGACAGAGCCCGGCACTGGGCCGGGCCCTTCCACTTTAGGCAAGCGCCAATCAGTTCCCGCAAACCGGCAGCTCCTTGCAGATGCCCTTGCCCTTCTCACGGACGGCGTCGATCACCGCTTCgttgtt                                 >  2:424334/1‑150 (MQ=255)
                                                                                                              gggCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTc                                <  2:1392935/150‑1 (MQ=255)
                                                                                                                            tATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTAACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCAACGCTTCGTTGTTCAGGTCggatgcgat                  >  1:568182/1‑143 (MQ=255)
                                                                                                                                  aCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTtctc                                                               >  2:823549/1‑99 (MQ=255)
                                                                                                                                  aCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTtctc                                                               <  1:823549/99‑1 (MQ=255)
                                                                                                                                  aCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTtct                                                                <  2:968196/98‑1 (MQ=255)
                                                                                                                                  aCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTtct                                                                >  1:968196/1‑98 (MQ=255)
                                                                                                                                    gAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCCGCCGATGGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGa          >  1:615943/1‑150 (MQ=255)
                                                                                                                                    gAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATgg                                                                            >  2:1099051/1‑84 (MQ=255)
                                                                                                                                    gAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATgg                                                                            <  1:1099051/84‑1 (MQ=255)
                                                                                                                                      aaGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCACGAGCGCATCAGAAGCCGTCACACGGTAGTTTCTTGCATATGGCCTTTACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATc        <  2:160967/150‑1 (MQ=255)
                                                                                                                                        tagCCCGTCACTGGTCCGGTCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCCTAAACCGGCAGCTTCTTGCAGATGGCCTTTACCTTCTCATGGACGTCGTCGATCACCGCGTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCa      <  2:568182/149‑1 (MQ=255)
                                                                                                                                         agCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAg     >  2:129865/1‑150 (MQ=255)
                                                                                                                                            ccGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCa                                              >  2:1424252/1‑106 (MQ=255)
                                                                                                                                            ccGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCa                                              <  1:1424252/106‑1 (MQ=255)
                                                                                                                                            ccGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGCCg  <  2:1056130/150‑1 (MQ=255)
                                                                                                                                                  cTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGa                   <  2:193163/127‑1 (MQ=255)
                                                                                                                                                  cTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGa                   >  1:193163/1‑127 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GGCCAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTGCGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGCCG  >  NC_002947/780414‑780703

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