Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,155,230 T→G 10.1% intergenic (‑12/‑125) yqhC ← / → yqhD DNA‑binding transcriptional activator YqhC/NADPH‑dependent aldehyde reductase YqhD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,155,2300TG10.1% 256.6 / 20.0 89intergenic (‑12/‑125)yqhC/yqhDDNA‑binding transcriptional activator YqhC/NADPH‑dependent aldehyde reductase YqhD
Reads supporting (aligned to +/- strand):  ref base T (29/51);  new base G (2/7);  total (31/58)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AGTCCTGACAAACTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATGAACAACTTTAATCTGCACACC  >  NC_000913/3155082‑3155378
                                                                                                                                                    |                                                                                                                                                    
aGTCCTGACAAACTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCActc                                                                                                                                                     <  1:667170/150‑1 (MQ=255)
  tCCTGACAAACTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTg                                                                                                                                                   <  2:762223/150‑1 (MQ=255)
   ccTGACAATCTATTTTCTTTATTTTTCAGTTTATTAACTTTATCGGCCACCAAGCGGGAAATCTCTTCACTTTTGATATCGCGTAATTTCTTAGCAATAATGCGGCACTTTGCTTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTgg                                                                                                                                                  <  2:824920/150‑1 (MQ=255)
   ccTGACAAACTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTgg                                                                                                                                                  <  2:1391172/150‑1 (MQ=255)
   ccTGACAAACTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTgg                                                                                                                                                  <  1:644990/150‑1 (MQ=255)
     tGACAAACTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGgag                                                                                                                                                <  1:1541368/150‑1 (MQ=255)
     tGACAAACTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGgag                                                                                                                                                <  2:234042/150‑1 (MQ=255)
        cAAACTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGaaa                                                                                                                                             <  2:1391079/150‑1 (MQ=255)
         aaaCTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAAt                                                                                                                                            >  2:734387/1‑150 (MQ=255)
          aaCTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATa                                                                                                                                           <  1:43339/150‑1 (MQ=255)
          aaCTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATa                                                                                                                                           >  1:479980/1‑150 (MQ=255)
          aaCTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGAGTGCTGGAGAAAtg                                                                                                                                           <  2:1454249/150‑2 (MQ=255)
           aCTATTTTCTTTATTTTTCAGTTTATTTACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGAGTGCTGGAGAAAtgc                                                                                                                                          <  1:1488445/150‑3 (MQ=255)
           aCTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGAGTGCTGGAGAAAtgc                                                                                                                                          <  1:1650020/150‑3 (MQ=255)
           aCTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGAGTGCTGGAGAAAtgc                                                                                                                                          <  2:1733354/150‑3 (MQ=255)
            cTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGAGTGCTGGAGAAAtgct                                                                                                                                         <  1:1471312/150‑4 (MQ=255)
            cTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGAGTGCTGGAGAAAtgct                                                                                                                                         <  1:1574118/150‑4 (MQ=255)
              aTTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGAGTGCTGGAGAAAtgctac                                                                                                                                       <  2:736149/150‑6 (MQ=255)
                tttCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCTCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCATCACTCTGGAGAAATAGGCAAg                                                                                                                                     <  2:1055768/150‑1 (MQ=255)
                       aTTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTgg                                                                                                                              >  1:786781/1‑150 (MQ=255)
                        tttttCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGc                                                                                                                             <  2:538879/150‑1 (MQ=255)
                          tttCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAg                                                                                                                           >  1:1691573/1‑150 (MQ=255)
                          tttCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAg                                                                                                                           >  1:1337237/1‑150 (MQ=255)
                           ttCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGa                                                                                                                          <  1:1814071/150‑1 (MQ=255)
                               gTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATg                                                                                                                      <  1:698509/150‑1 (MQ=255)
                                tttATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGa                                                                                                                     >  1:1247119/1‑150 (MQ=255)
                                  tattaACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGc                                                                                                                   >  2:946897/1‑150 (MQ=255)
                                  tattaACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGc                                                                                                                   >  1:334111/1‑150 (MQ=255)
                                   attaACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCa                                                                                                                  >  2:1299263/1‑150 (MQ=255)
                                     taACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCAtt                                                                                                                <  1:493470/150‑1 (MQ=255)
                                        cTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTgag                                                                                                             <  1:81643/150‑1 (MQ=255)
                                         tttATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTgaga                                                                                                            >  1:1974735/1‑150 (MQ=255)
                                           tATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGc                                                                                                          >  1:149248/1‑150 (MQ=255)
                                              ccGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAg                                                                                                       >  2:1953555/1‑150 (MQ=255)
                                                gCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAggg                                                                                                     <  1:946897/150‑1 (MQ=255)
                                                   aGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGgcgc                                                                                                  >  1:1595944/1‑150 (MQ=255)
                                                         cGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGa                                                                                            >  2:1374578/1‑150 (MQ=255)
                                                          ggCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGAt                                                                                           <  2:555266/150‑1 (MQ=255)
                                                           gCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATc                                                                                          <  2:1369707/150‑1 (MQ=255)
                                                           gCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATc                                                                                          <  2:98665/150‑1 (MQ=255)
                                                           gCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATc                                                                                          <  2:786781/150‑1 (MQ=255)
                                                           gCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATc                                                                                          <  1:1012673/150‑1 (MQ=255)
                                                             aaaTCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATcac                                                                                        <  1:223602/150‑1 (MQ=255)
                                                             aaaTCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATcac                                                                                        <  1:1495391/150‑1 (MQ=255)
                                                               aTCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACaa                                                                                      <  2:28362/150‑1 (MQ=255)
                                                                      cACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCAATCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGaaaaa                                                                               >  1:1902721/1‑150 (MQ=255)
                                                                       aCGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAAc                                                                              <  1:1149030/150‑1 (MQ=255)
                                                                           ttCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCa                                                                          <  1:1221430/150‑1 (MQ=255)
                                                                              atatCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGc                                                                       <  1:1303482/150‑1 (MQ=255)
                                                                              atatCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGc                                                                       <  1:1964765/150‑1 (MQ=255)
                                                                               tatCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCa                                                                      >  1:1818369/1‑150 (MQ=255)
                                                                               tatCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCa                                                                      <  1:1881456/150‑1 (MQ=255)
                                                                                  cgcgTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGAt                                                                   >  2:1763105/1‑150 (MQ=255)
                                                                                   gcgTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATc                                                                  <  1:734387/150‑1 (MQ=255)
                                                                                   gcgTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATc                                                                  >  1:1086083/1‑150 (MQ=255)
                                                                                     gTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGt                                                                >  1:354361/1‑150 (MQ=255)
                                                                                       aaTTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCGGGCAGATCGTtc                                                              <  1:1570465/150‑1 (MQ=255)
                                                                                       aaTTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTtc                                                              <  2:457277/150‑1 (MQ=255)
                                                                                         tttCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTtctc                                                            <  1:1093769/150‑1 (MQ=255)
                                                                                         tttCTTAGGAATAATGAGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAAAAGGCAAGACATTGCCAGAAATGAGCAATGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGAACGTtctc                                                            >  2:1190293/1‑150 (MQ=255)
                                                                                          ttCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTctct                                                           <  2:1462826/150‑1 (MQ=255)
                                                                                             ttAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGcc                                                        >  2:374998/1‑150 (MQ=255)
                                                                                                ggAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTc                                                     <  1:1374578/150‑1 (MQ=255)
                                                                                                ggAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTc                                                     >  1:759262/1‑150 (MQ=255)
                                                                                                   ataatGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCata                                                  <  2:334111/150‑1 (MQ=255)
                                                                                                      atGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTg                                               <  2:864055/150‑1 (MQ=255)
                                                                                                        gCGGCAATTTGATTGTTCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAAGATCAGGCAGATCGTTCTCTGCCCTCATATTGGc                                             <  1:1193498/150‑1 (MQ=255)
                                                                                                         cGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGcc                                            <  1:1299263/150‑1 (MQ=255)
                                                                                                            cAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAg                                         <  1:1741094/150‑1 (MQ=255)
                                                                                                                  gATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAgg                                   <  1:1717601/150‑1 (MQ=255)
                                                                                                                      gtgCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGc                               <  2:1086083/150‑1 (MQ=255)
                                                                                                                         cgcAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAg                            <  2:588095/150‑1 (MQ=255)
                                                                                                                         cgcAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAg                            <  2:149248/150‑1 (MQ=255)
                                                                                                                          gcAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGt                           <  2:639018/150‑1 (MQ=255)
                                                                                                                             aTTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAAt                        <  1:99856/150‑1 (MQ=255)
                                                                                                                             aTTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAAt                        <  2:117691/150‑1 (MQ=255)
                                                                                                                              ttttGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATg                       <  1:937903/150‑1 (MQ=255)
                                                                                                                              ttttGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATg                       <  1:1239741/150‑1 (MQ=255)
                                                                                                                                    cctaTTTCTCCAGAGTGCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATGaacaac                 >  2:362632/4‑150 (MQ=255)
                                                                                                                                     ctaTTTCTCCAGAGTGCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATGAACAACt                >  2:328521/3‑150 (MQ=255)
                                                                                                                                      cATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCGCAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATGAACAACtt               >  2:1854129/1‑150 (MQ=255)
                                                                                                                                          tctcCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATGAACAACTTTAAt           >  1:136433/1‑150 (MQ=255)
                                                                                                                                           ctcCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATGAACAACTTTAATc          >  2:648433/1‑150 (MQ=255)
                                                                                                                                           ctcCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATGAACAACTTTAATc          <  2:1818369/150‑1 (MQ=255)
                                                                                                                                           ctcCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGCCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGCCCCAGCAAACGGAGGAAGTAATGAACAACTTTAATc          >  2:655027/1‑150 (MQ=255)
                                                                                                                                            tcCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATGaac                    >  1:1434611/1‑139 (MQ=255)
                                                                                                                                            tcCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATGaac                    <  2:1434611/139‑1 (MQ=255)
                                                                                                                                             ccAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATGAACAACTTTAATCTg        >  1:633039/1‑150 (MQ=255)
                                                                                                                                              cAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATGAACAACTTTAATCTGc       >  1:374466/1‑150 (MQ=255)
                                                                                                                                               aGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATGAACAACTTTAATCTGca      <  2:1974735/150‑1 (MQ=255)
                                                                                                                                                   ctctGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATGAACAACTTTAATCTGCACAcc  >  2:68757/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGTCCTGACAAACTATTTTCTTTATTTTTCAGTTTATTAACTTTATCCGCCAGCAAGCGGCAAATCTCTTCACGTTTCATATCGCGTAATTTCTTAGGAATAATGCGGCAATTTGATTGTGCGCAATTTTGTAGCATTTCTCCAGCACTCTGGAGAAATAGGCAAGACATTGGCAGAAATGAGCATTGAGAGCCAGGGCGCTGGCGATCACAATGAAAAACATCAGGCAGATCGTTCTCTGCCCTCATATTGGCCCAGCAAAGGGAGCAAGTAATGAACAACTTTAATCTGCACACC  >  NC_000913/3155082‑3155378

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