Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,726,514 C→A 5.4% Y164* (TAC→TAA) ‡ 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,5140CA 5.4% 287.2 / 4.3 91Y164* (TAC→TAA) ‡nemRtranscriptional repressor for the nemRA‑gloA operon, quinone‑, glyoxal‑, and HOCl‑activated
Reads supporting (aligned to +/- strand):  ref base C (36/50);  new base A (2/3);  total (38/53)
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.85e-01

GTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGCCTTTATGTC  >  NC_000913/1726365‑1726663
                                                                                                                                                     |                                                                                                                                                      
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                                                                                       <  2:667874/150‑1 (MQ=255)
                                                                tGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAAcg                                                                                       <  1:1143479/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                                                                               <  2:1390419/150‑1 (MQ=255)
                                                                        cAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGcc                                                                               >  1:652780/1‑150 (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                                                                   <  1:624802/150‑1 (MQ=255)
                                                                                    aaTGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAAc                                                                   <  2:415204/150‑1 (MQ=255)
                                                                                    aaTGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAAc                                                                   <  2:1074273/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:1231301/1‑150 (MQ=255)
                                                                                               gaACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGAc                                                        >  2:1019321/1‑150 (MQ=255)
                                                                                                aaCCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTgg                                                                                             <  2:1104673/112‑1 (MQ=255)
                                                                                                aaCCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTgg                                                                                             >  1:1104673/1‑112 (MQ=255)
                                                                                                  ccATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCg                                                          <  2:730620/145‑1 (MQ=255)
                                                                                                  ccATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCg                                                          >  1:730620/1‑145 (MQ=255)
                                                                                                  ccATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGcc                                                     <  1:1231301/150‑1 (MQ=255)
                                                                                                  ccATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGcc                                                     <  1:1019321/150‑1 (MQ=255)
                                                                                                     ttGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTGCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGc                                                  <  1:902513/150‑1 (MQ=255)
                                                                                                           aaCCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTa                                            <  1:995271/150‑1 (MQ=255)
                                                                                                              cTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCa                                         >  2:589388/1‑150 (MQ=255)
                                                                                                              cTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCa                                         <  1:569973/150‑1 (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:614719/1‑150 (MQ=255)
                                                                                                                  tgtgGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCa                                     <  2:227009/150‑1 (MQ=255)
                                                                                                                   gtgGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAt                                    >  1:1418805/1‑150 (MQ=255)
                                                                                                                   gtgGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAt                                    <  1:1255317/150‑1 (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)
                                                                                                                                                   tACGCAC‑TGGGGCTTGGCGCGAATCTGCAGGCCAAATTTTCGCGCATTTTCTTGCCACTGGAAAACGCGCTGGACCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGCCTTTATg    <  1:767437/150‑1 (MQ=255)
                                                                                                                                                     cGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGCCTTTATGTc  <  2:542606/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCAC‑TGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGCCTTTATGTC  >  NC_000913/1726365‑1726663

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: ‑