Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,726,510 T→G 6.3% L163R (CTT→CGT)  nemR → transcriptional repressor for the nemRA‑gloA operon, quinone‑, glyoxal‑, and HOCl‑activated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,726,5100TG 6.3% 265.2 / 6.8 91L163R (CTT→CGT) nemRtranscriptional repressor for the nemRA‑gloA operon, quinone‑, glyoxal‑, and HOCl‑activated
Reads supporting (aligned to +/- strand):  ref base T (34/51);  new base G (3/3);  total (37/54)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

AAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGCCTTTA  >  NC_000913/1726363‑1726659
                                                                                                                                                   |                                                                                                                                                      
aaGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCttt                                                                                                                                                      <  1:1702028/150‑1 (MQ=255)
  gtgtGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTAc                                                                                                                                                    <  1:179342/150‑1 (MQ=255)
  gtgtGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTAc                                                                                                                                                    <  1:1284104/150‑1 (MQ=255)
    gtgCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGc                                                                                                                                                  >  1:506556/1‑150 (MQ=255)
       cGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑t                                                                                                                                              >  2:283202/1‑150 (MQ=255)
          tCTGCCTGAACAAATCCGCCGCGCGTGGCATACAGGTGCTCGCGGGGTGCCCGCCCTGCTCGCGCCGGCGCGGTAACATGGCCGGCACAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGTCACTAGTGCTTTACGCAC‑tgtg                                                                                                                                           <  2:1596775/150‑1 (MQ=255)
            tGTCTGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGc                                                                                                                                         <  1:90309/150‑1 (MQ=255)
            tGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGc                                                                                                                                         <  1:666093/150‑1 (MQ=255)
            tGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCAC‑AGTGCGTAAAGCAC‑ttgtgcc                                                                                                                                        <  1:117384/150‑7 (MQ=255)
            tGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCAC‑AGTGCGTAAAGCAC‑ttgtgcc                                                                                                                                        <  1:1527603/150‑7 (MQ=255)
             gTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGATCCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCt                                                                                                                                        >  1:864635/1‑150 (MQ=255)
             gTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCGTAAAGCAC‑ttgtgcc                                                                                                                                        <  2:465645/150‑7 (MQ=255)
              tCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCtt                                                                                                                                       <  2:407331/150‑1 (MQ=255)
                agaagaTATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTgg                                                                                                                                     <  2:506556/150‑1 (MQ=255)
                 gaagaTATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGgc                                                                                                                                    >  1:987426/1‑150 (MQ=255)
                 gaagaTATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGgc                                                                                                                                    <  2:1710363/150‑1 (MQ=255)
                   agaTATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGgcgc                                                                                                                                  >  2:668234/1‑150 (MQ=255)
                       atGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAAt                                                                                                                              >  2:962588/1‑150 (MQ=255)
                        tGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGGGATCGCCCTGCTCTCGCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATc                                                                                                                             <  1:291635/150‑1 (MQ=255)
                          cgcAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTg                                                                                                                           <  2:279174/150‑1 (MQ=255)
                           gcAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGc                                                                                                                          <  1:1018331/150‑1 (MQ=255)
                            cAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTTTTTAACCTTTTGTGGCGAACCTCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTCGCTTGGCGCGAATCTGCa                                                                                                                         >  2:1073495/1‑150 (MQ=255)
                            cAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCa                                                                                                                         >  1:667874/1‑150 (MQ=255)
                              gcgcGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAgg                                                                                                                       >  2:45456/1‑150 (MQ=255)
                              gcgcGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAgg                                                                                                                       <  2:791756/150‑1 (MQ=255)
                                   aTGGATAAAGGTGCTTGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGc                                                                                                                          >  2:1242579/1‑142 (MQ=255)
                                   aTGGATAAAGGTGCTTGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGc                                                                                                                          <  1:1242579/142‑1 (MQ=255)
                                     ggATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAAt                                                                                                                >  1:400317/1‑150 (MQ=255)
                                       aTAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAAttt                                                                                                              <  1:1255622/150‑1 (MQ=255)
                                            ggTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCgcgc                                                                                                         >  1:542606/1‑150 (MQ=255)
                                                cTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCGTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGtt                                                                                                     >  1:1688455/1‑150 (MQ=255)
                                                  cgcgGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTc                                                                                                   >  1:399588/1‑150 (MQ=255)
                                                   gcggcgTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCg                                                                                                  <  1:144171/150‑1 (MQ=255)
                                                     ggcgTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAg                                                                                                <  2:1566242/150‑1 (MQ=255)
                                                     ggcgTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAg                                                                                                >  2:995271/1‑150 (MQ=255)
                                                      gcgTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGc                                                                                               >  1:1583405/1‑150 (MQ=255)
                                                         tGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACATTTTGTGTCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGTCGCGAATCTGCAGGCCAAAATTTCGAGCAGTTTCGAGCCAc                                                                                            >  2:628000/1‑150 (MQ=255)
                                                                ccTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGAAGTTTCGAGCCACTGGaaaa                                                                                     <  2:400317/150‑1 (MQ=255)
                                                                  tGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAAcg                                                                                   <  1:1143479/150‑1 (MQ=255)
                                                                  tGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAAcg                                                                                   <  2:667874/150‑1 (MQ=255)
                                                                   gCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAAcgc                                                                                  >  1:1393879/1‑150 (MQ=255)
                                                                    ctctctCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAAcgcg                                                                                 <  1:529274/150‑1 (MQ=255)
                                                                     tctctcAGGCGCTGGAAACTGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGACTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACgcgc                                                                                <  1:1060576/150‑1 (MQ=255)
                                                                      ctctcAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCAGCAGGCCAAAATTTCGCGAAGTTTCTAGCCACTGGAAAAAGCGCt                                                                               >  2:767437/1‑150 (MQ=255)
                                                                      ctctcAGGCGCTGGAAAATGGCCGTGAGAAACATTGTTTAACCTGTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCt                                                                               >  1:544363/1‑150 (MQ=255)
                                                                          cAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGcc                                                                           >  1:652780/1‑150 (MQ=255)
                                                                          cAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGcc                                                                           <  2:1390419/150‑1 (MQ=255)
                                                                            ggCGCTGGAAAATGGCCGTGAGAACCATTGTTTTACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCa                                                                         <  1:403320/150‑1 (MQ=255)
                                                                                cTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTa                                                                     <  1:1157747/150‑1 (MQ=255)
                                                                                    aaaaTGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTaaaaa                                                                 <  1:1397113/150‑1 (MQ=255)
                                                                                      aaTGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAAc                                                               <  2:1074273/150‑1 (MQ=255)
                                                                                      aaTGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAAc                                                               <  2:415204/150‑1 (MQ=255)
                                                                                      aaTGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAAc                                                               <  1:624802/150‑1 (MQ=255)
                                                                                           ccGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTATAAAACattat                                                          >  2:647939/1‑150 (MQ=255)
                                                                                            cGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACattatt                                                         <  2:652780/150‑1 (MQ=255)
                                                                                              tGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGc                                                       <  2:399588/150‑1 (MQ=255)
                                                                                               gagaACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCg                                                      <  2:577454/150‑1 (MQ=255)
                                                                                                 taaCCATTGTTTAAACTTCTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGAc                                                    <  2:853156/149‑1 (MQ=255)
                                                                                                 gaACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGAc                                                    >  2:1019321/1‑150 (MQ=255)
                                                                                                 gaACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGAc                                                    >  2:1231301/1‑150 (MQ=255)
                                                                                                  aaCCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTgg                                                                                         >  1:1104673/1‑112 (MQ=255)
                                                                                                  aaCCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTgg                                                                                         <  2:1104673/112‑1 (MQ=255)
                                                                                                    ccATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCg                                                      <  2:730620/145‑1 (MQ=255)
                                                                                                    ccATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCg                                                      >  1:730620/1‑145 (MQ=255)
                                                                                                    ccATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGcc                                                 <  1:1019321/150‑1 (MQ=255)
                                                                                                    ccATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGcc                                                 <  1:1231301/150‑1 (MQ=255)
                                                                                                       ttGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTGCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGc                                              <  1:902513/150‑1 (MQ=255)
                                                                                                             aaCCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTa                                        <  1:995271/150‑1 (MQ=255)
                                                                                                                cTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCa                                     <  1:569973/150‑1 (MQ=255)
                                                                                                                cTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCa                                     >  2:589388/1‑150 (MQ=255)
                                                                                                                 ttttGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATATCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAg                                    <  1:104633/150‑1 (MQ=255)
                                                                                                                  tttGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAgg                                   <  2:848466/150‑1 (MQ=255)
                                                                                                                    tgtgGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCa                                 >  2:327896/1‑150 (MQ=255)
                                                                                                                    tgtgGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCa                                 <  2:227009/150‑1 (MQ=255)
                                                                                                                    tgtgGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCa                                 >  2:614719/1‑150 (MQ=255)
                                                                                                                     gtgGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAt                                <  1:1255317/150‑1 (MQ=255)
                                                                                                                     gtgGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAt                                >  1:1418805/1‑150 (MQ=255)
                                                                                                                      tgGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCATGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAtt                               >  1:1091404/1‑150 (MQ=255)
                                                                                                                        gCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAtttt                             >  2:928960/1‑150 (MQ=255)
                                                                                                                         cGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAttttt                            <  2:987426/150‑1 (MQ=255)
                                                                                                                           aaCCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTa                          <  2:1455308/150‑1 (MQ=255)
                                                                                                                            aCCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTAt                         >  2:1387626/1‑150 (MQ=255)
                                                                                                                                  gCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAg                   <  2:1393879/150‑1 (MQ=255)
                                                                                                                                     aCAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCACATGTAAAAAACAATATTTCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAgacg                >  2:272556/1‑150 (MQ=255)
                                                                                                                                        ggCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAACCGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGAcc             >  2:1178884/1‑150 (MQ=255)
                                                                                                                                         gccacaGTGCGTAAAGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCg            >  1:1111688/6‑150 (MQ=255)
                                                                                                                                          ccacaGTGCGTAAAGCACTTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCg            >  2:383240/5‑150 (MQ=255)
                                                                                                                                           acaAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTCCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCggg          <  2:573590/150‑1 (MQ=255)
                                                                                                                                            caAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGa         <  1:45456/150‑1 (MQ=255)
                                                                                                                                               gTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGcc      <  1:1599964/150‑1 (MQ=255)
                                                                                                                                                ggCTTTACGCAC‑TGTGGGTTGGCTCTAATCTGCAGGCCAAAATTTCTCGCAGTTTCGAGCCACTGTAAAACGCGCTGGCCCATGTAAAAACCATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCATACGACCGGGAGCCt     <  2:544363/149‑1 (MQ=255)
                                                                                                                                                  cTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGTCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAAGATTTTTGCGACGCCTGCCGTTTAGCTGGCATTTTTTATCACCAGACGACCGGGAGCCttt   <  1:628000/150‑1 (MQ=255)
                                                                                                                                                   tttACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGCCTTTa  >  2:1695121/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
AAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGCCTTTA  >  NC_000913/1726363‑1726659

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