Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,726,512 T→A 5.5% Y164N (TAC→AAC) ‡ 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,5120TA 5.5% 269.0 / 4.3 91Y164N (TAC→AAC) ‡nemRtranscriptional repressor for the nemRA‑gloA operon, quinone‑, glyoxal‑, and HOCl‑activated
Reads supporting (aligned to +/- strand):  ref base T (35/51);  new base A (2/3);  total (37/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

AAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGCCTTTATG  >  NC_000913/1726363‑1726661
                                                                                                                                                     |                                                                                                                                                      
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)
       cGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACt                                                                                                                                                >  2:283202/1‑150 (MQ=255)
          tCTGCCTGAACAAATCCGCCGCGCGTGGCATACAGGTGCTCGCGGGGTGCCCGCCCTGCTCGCGCCGGCGCGGTAACATGGCCGGCACAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGTCACTAGTGCTTTACGCACtgtg                                                                                                                                             <  2:1596775/150‑1 (MQ=255)
            tGTCTGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGc                                                                                                                                           <  1:90309/150‑1 (MQ=255)
            tGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGc                                                                                                                                           <  1:666093/150‑1 (MQ=255)
            tGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAGTGCGTAAAGCACttgtgcc                                                                                                                                          <  1:117384/150‑7 (MQ=255)
            tGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAGTGCGTAAAGCACttgtgcc                                                                                                                                          <  1:1527603/150‑7 (MQ=255)
             gTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGATCCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCt                                                                                                                                          >  1:864635/1‑150 (MQ=255)
             gTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCGTAAAGCACttgtgcc                                                                                                                                          <  2:465645/150‑7 (MQ=255)
              tCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCtt                                                                                                                                         <  2:407331/150‑1 (MQ=255)
                agaagaTATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTgg                                                                                                                                       <  2:506556/150‑1 (MQ=255)
                 gaagaTATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGgc                                                                                                                                      >  1:987426/1‑150 (MQ=255)
                 gaagaTATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGgc                                                                                                                                      <  2:1710363/150‑1 (MQ=255)
                   agaTATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGgcgc                                                                                                                                    >  2:668234/1‑150 (MQ=255)
                       atGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAAt                                                                                                                                >  2:962588/1‑150 (MQ=255)
                        tGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGGGATCGCCCTGCTCTCGCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATc                                                                                                                               <  1:291635/150‑1 (MQ=255)
                          cgcAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTg                                                                                                                             <  2:279174/150‑1 (MQ=255)
                           gcAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGc                                                                                                                            <  1:1018331/150‑1 (MQ=255)
                            cAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTTTTTAACCTTTTGTGGCGAACCTCTGCAACAGGCACAAGTGCTTTACGCACTGTCGCTTGGCGCGAATCTGCa                                                                                                                           >  2:1073495/1‑150 (MQ=255)
                            cAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCa                                                                                                                           >  1:667874/1‑150 (MQ=255)
                              gcgcGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAgg                                                                                                                         <  2:791756/150‑1 (MQ=255)
                              gcgcGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAgg                                                                                                                         >  2:45456/1‑150 (MQ=255)
                                   aTGGATAAAGGTGCTTGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGc                                                                                                                            <  1:1242579/142‑1 (MQ=255)
                                   aTGGATAAAGGTGCTTGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGc                                                                                                                            >  2:1242579/1‑142 (MQ=255)
                                     ggATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAAt                                                                                                                  >  1:400317/1‑150 (MQ=255)
                                       aTAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAAttt                                                                                                                <  1:1255622/150‑1 (MQ=255)
                                            ggTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCgcgc                                                                                                           >  1:542606/1‑150 (MQ=255)
                                                cTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCGTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGtt                                                                                                       >  1:1688455/1‑150 (MQ=255)
                                                  cgcgGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTc                                                                                                     >  1:399588/1‑150 (MQ=255)
                                                   gcggcgTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCg                                                                                                    <  1:144171/150‑1 (MQ=255)
                                                     ggcgTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAg                                                                                                  >  2:995271/1‑150 (MQ=255)
                                                     ggcgTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAg                                                                                                  <  2:1566242/150‑1 (MQ=255)
                                                      gcgTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGc                                                                                                 >  1:1583405/1‑150 (MQ=255)
                                                         tGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACATTTTGTGTCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGTCGCGAATCTGCAGGCCAAAATTTCGAGCAGTTTCGAGCCAc                                                                                              >  2:628000/1‑150 (MQ=255)
                                                                ccTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGAAGTTTCGAGCCACTGGaaaa                                                                                       <  2:400317/150‑1 (MQ=255)
                                                                  tGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAAcg                                                                                     <  1:1143479/150‑1 (MQ=255)
                                                                  tGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAAcg                                                                                     <  2:667874/150‑1 (MQ=255)
                                                                   gCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAAcgc                                                                                    >  1:1393879/1‑150 (MQ=255)
                                                                    ctctctCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAAcgcg                                                                                   <  1:529274/150‑1 (MQ=255)
                                                                     tctctcAGGCGCTGGAAACTGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGACTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACgcgc                                                                                  <  1:1060576/150‑1 (MQ=255)
                                                                      ctctcAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCAGCAGGCCAAAATTTCGCGAAGTTTCTAGCCACTGGAAAAAGCGCt                                                                                 >  2:767437/1‑150 (MQ=255)
                                                                      ctctcAGGCGCTGGAAAATGGCCGTGAGAAACATTGTTTAACCTGTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCt                                                                                 >  1:544363/1‑150 (MQ=255)
                                                                          cAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGcc                                                                             >  1:652780/1‑150 (MQ=255)
                                                                          cAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGcc                                                                             <  2:1390419/150‑1 (MQ=255)
                                                                            ggCGCTGGAAAATGGCCGTGAGAACCATTGTTTTACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCa                                                                           <  1:403320/150‑1 (MQ=255)
                                                                                cTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTa                                                                       <  1:1157747/150‑1 (MQ=255)
                                                                                    aaaaTGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTaaaaa                                                                   <  1:1397113/150‑1 (MQ=255)
                                                                                      aaTGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAAc                                                                 <  2:1074273/150‑1 (MQ=255)
                                                                                      aaTGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAAc                                                                 <  2:415204/150‑1 (MQ=255)
                                                                                      aaTGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAAc                                                                 <  1:624802/150‑1 (MQ=255)
                                                                                           ccGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTATAAAACattat                                                            >  2:647939/1‑150 (MQ=255)
                                                                                            cGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACattatt                                                           <  2:652780/150‑1 (MQ=255)
                                                                                              tGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGc                                                         <  2:399588/150‑1 (MQ=255)
                                                                                               gagaACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCg                                                        <  2:577454/150‑1 (MQ=255)
                                                                                                 taaCCATTGTTTAAACTTCTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGAc                                                      <  2:853156/149‑1 (MQ=255)
                                                                                                 gaACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGAc                                                      >  2:1019321/1‑150 (MQ=255)
                                                                                                 gaACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGAc                                                      >  2:1231301/1‑150 (MQ=255)
                                                                                                  aaCCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTgg                                                                                           >  1:1104673/1‑112 (MQ=255)
                                                                                                  aaCCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTgg                                                                                           <  2:1104673/112‑1 (MQ=255)
                                                                                                    ccATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCg                                                        >  1:730620/1‑145 (MQ=255)
                                                                                                    ccATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCg                                                        <  2:730620/145‑1 (MQ=255)
                                                                                                    ccATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGcc                                                   <  1:1019321/150‑1 (MQ=255)
                                                                                                    ccATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGcc                                                   <  1:1231301/150‑1 (MQ=255)
                                                                                                       ttGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTGCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGc                                                <  1:902513/150‑1 (MQ=255)
                                                                                                             aaCCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTa                                          <  1:995271/150‑1 (MQ=255)
                                                                                                                cTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCa                                       >  2:589388/1‑150 (MQ=255)
                                                                                                                cTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCa                                       <  1:569973/150‑1 (MQ=255)
                                                                                                                 ttttGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATATCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAg                                      <  1:104633/150‑1 (MQ=255)
                                                                                                                  tttGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAgg                                     <  2:848466/150‑1 (MQ=255)
                                                                                                                    tgtgGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCa                                   >  2:327896/1‑150 (MQ=255)
                                                                                                                    tgtgGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCa                                   <  2:227009/150‑1 (MQ=255)
                                                                                                                    tgtgGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCa                                   >  2:614719/1‑150 (MQ=255)
                                                                                                                     gtgGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAt                                  <  1:1255317/150‑1 (MQ=255)
                                                                                                                     gtgGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAt                                  >  1:1418805/1‑150 (MQ=255)
                                                                                                                      tgGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCATGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAtt                                 >  1:1091404/1‑150 (MQ=255)
                                                                                                                        gCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAtttt                               >  2:928960/1‑150 (MQ=255)
                                                                                                                         cGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAttttt                              <  2:987426/150‑1 (MQ=255)
                                                                                                                           aaCCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTa                            <  2:1455308/150‑1 (MQ=255)
                                                                                                                            aCCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTAt                           >  2:1387626/1‑150 (MQ=255)
                                                                                                                                  gCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAg                     <  2:1393879/150‑1 (MQ=255)
                                                                                                                                     aCAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCACATGTAAAAAACAATATTTCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAgacg                  >  2:272556/1‑150 (MQ=255)
                                                                                                                                        ggCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAACCGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGAcc               >  2:1178884/1‑150 (MQ=255)
                                                                                                                                         gccacaGTGCGTAAAGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCg              >  1:1111688/6‑150 (MQ=255)
                                                                                                                                          ccacaGTGCGTAAAGCACTTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCg              >  2:383240/5‑150 (MQ=255)
                                                                                                                                           acaAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTCCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCggg            <  2:573590/150‑1 (MQ=255)
                                                                                                                                            caAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGa           <  1:45456/150‑1 (MQ=255)
                                                                                                                                               gTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGcc        <  1:1599964/150‑1 (MQ=255)
                                                                                                                                                ggCTTTACGCACTGTGGGTTGGCTCTAATCTGCAGGCCAAAATTTCTCGCAGTTTCGAGCCACTGTAAAACGCGCTGGCCCATGTAAAAACCATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCATACGACCGGGAGCCt       <  2:544363/149‑1 (MQ=255)
                                                                                                                                                  cTTTACGCACTGTGGCTTGGCGCGAATCTGCAGTCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAAGATTTTTGCGACGCCTGCCGTTTAGCTGGCATTTTTTATCACCAGACGACCGGGAGCCttt     <  1:628000/150‑1 (MQ=255)
                                                                                                                                                   tttACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGCCTTTa    >  2:1695121/1‑150 (MQ=255)
                                                                                                                                                     tACGCACTGGGGCTTGGCGCGAATCTGCAGGCCAAATTTTCGCGCATTTTCTTGCCACTGGAAAACGCGCTGGACCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGCCTTTATg  <  1:767437/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AAGTGTGCGATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGACGACCGGGAGCCTTTATG  >  NC_000913/1726363‑1726661

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: