Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,707 C→A intergenic (+160/+148) nadR → / ← ettA DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,628,7070CA97.1% 196.7 / NA 71intergenic (+160/+148)nadR/ettADNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base C (2/0);  major base A (35/32);  minor base G (0/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGT‑AGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAA  >  NC_000913/4628468‑4628944
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
gggTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCAcgcc                                                                                                                                                                                                                                   >  1:156055/1‑251 (MQ=255)
    tGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCAcgccgtat                                                                                                                                                                                                                               >  1:238556/1‑248 (MQ=255)
     gAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCAcgccgtat                                                                                                                                                                                                                               <  2:103129/250‑4 (MQ=255)
             cgATTAGGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGCCCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTc                                                                                                                                                                                                                      <  1:14936/250‑1 (MQ=255)
                 taggACAGCCGTTTCCGGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTAAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCAGGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGt                                                                                                                                                                                                                   <  2:238556/247‑1 (MQ=255)
                              tCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGCTGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGTTGCTCAATGCCTg                                                                                                                                                                                                       >  2:53551/1‑249 (MQ=255)
                               ccggcgcggcgTGGAGCGGGGGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGGGAGGTATAAAATAAGTTTTTG‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATg                                                                                                                                                                                                    <  2:199824/248‑1 (MQ=255)
                                    gctgcgTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTTAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGTATACTAAGGCTGTTGT‑AGGCCGGAAAAGGTGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACg                                                                                                                                                                                               <  1:53551/251‑1 (MQ=255)
                                         gTGGAGCTGGTGCGGGAGCTGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGACTGATGCGACGCTTAc                                                                                                                                                                                          >  2:442066/1‑251 (MQ=255)
                                          tggagctcgtcagggAGATGATGGGGGAGCAGAGATAACCGCGCTGCAACGGCTCAAAGGCGAGGTATCAAATAAGTTTTTT‑GATGAGTTGACAACCTCTCTGGAAGAGGCTGTCGAGATTAACCACGGTTTGAAAAAACCTGCACGGGGGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGCAAAGGCGTTCCCGCCGCATCCGGCATTCCGTGCTCAATGCATGATGCGCCGCttccc                                                                                                                                                                                         >  2:58016/13‑248 (MQ=255)
                                                 ggTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCtt                                                                                                                                                                                  <  1:238162/251‑1 (MQ=255)
                                                          gatgatGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAACACCTGCACGGGTGACCCGCTACGAACTTGAGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGGGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCt                                                                                                                                                                         >  2:357324/1‑251 (MQ=255)
                                                                gggggAACAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATAAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCGGGAAAAGGCTTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCtggaga                                                                                                                                                                    <  2:341776/250‑6 (MQ=255)
                                                                 ggggAGCAGAGATAACCGCGATGAAACGGCTCGAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCtggagatt                                                                                                                                                                  <  1:460295/251‑8 (MQ=255)
                                                                   ggAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCtggagattc                                                                                                                                                                 >  1:477268/1‑242 (MQ=255)
                                                                    gagacgagaTAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCAGCTACGAAATTGCGGATACTAAGGCTGTTGA‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCtggcgattcac                                                                                                                                                               >  2:352596/6‑241 (MQ=255)
                                                                                   ggatgaAACGGCTAAAAGGCGAGGTATAAAAT‑AATTTTTTTGAAGAGTTGAAAACCTCTCGGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAAttt                                                                                                                                                <  1:186740/250‑2 (MQ=255)
                                                                                       gaaacgtcTCAAAGGCGATGTATAAAATAAGTTTTTG‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGTTTTGAAAAAGCCTGCACGGGTGACCCGCTACGAAATTGCGTATACTAAGGCTGTGGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAAtttctt                                                                                                                                             <  1:352596/243‑5 (MQ=255)
                                                                                       gaAACGGCTCAAAGGCGAGGTAGAAAATAAGTTTTTT‑GAGGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAAtttcttg                                                                                                                                            <  1:308118/251‑6 (MQ=255)
                                                                                         aacggctaaaaGGTGAGGTATAAAA‑GAAGTTTTTTGATGAGTGGAAAACCTCTCTGGAAGAGGCTGTAGAGATGAAACAAGGTTTTAAAAAACCTGCCCGGGTGACCCGCTACGAAATTTCGGATGCTAAGGCTGTTGT‑AGTCCGTAAAAGGCGTTCACGCCTTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAAtttcttggt                                                                                                                                          <  2:531537/243‑8 (MQ=255)
                                                                                         aacggCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGCAAAACTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGCAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAAtttcttggt                                                                                                                                          >  1:392602/5‑244 (MQ=255)
                                                                                             gCTCAAAGGCAGGGTAAAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGTCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAATTTGCGGATACTAAGTCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAAtttcttggtttgt                                                                                                                                      <  2:171481/251‑12 (MQ=255)
                                                                                              ctaaaaGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGACAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAACGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGCAAtttcttggtttgt                                                                                                                                      >  2:185724/4‑239 (MQ=255)
                                                                                              cTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAAtttcttggtttgt                                                                                                                                      >  1:8248/1‑239 (MQ=255)
                                                                                                   aGGCGTGGTATAAAATAAGTTTTTT‑GAGGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATCCTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCg                                                                                                                                <  2:474966/251‑1 (MQ=255)
                                                                                                       gAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGAt                                                                                                                             >  1:228279/1‑250 (MQ=255)
                                                                                                                aaataagttttttGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGATATTAAACAAGTGTTGAAAAAACCTGCACGGGGGACCCGCTACGAAATTGCTGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGAGTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTc                                                                                                                    <  1:378969/244‑1 (MQ=255)
                                                                                                                          tttGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCa                                                                                                                                                                                                                                       <  2:258674/126‑1 (MQ=255)
                                                                                                                          tttGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCa                                                                                                                                                                                                                                       >  1:258674/1‑126 (MQ=255)
                                                                                                                           ttGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGTCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                         >  1:554936/1‑251 (MQ=255)
                                                                                                                           tggaggaGTTAAAACCCTCTCGGGAAGGGGCTGTCGAAAGTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAGGGCTGTTGT‑AGGCGGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                          <  2:394139/245‑1 (MQ=255)
                                                                                                                            tgatgaGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAACCCTGCACGGGTGCCCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                         <  1:158455/250‑1 (MQ=255)
                                                                                                                             gatgaGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACACGGTTTGAAAAAACCTGCACGGGTGACCCGCTAAGAAATTGCGGATCCTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGACTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGGGATTCCCAACGTATAGAAATTTCTTGGTTTGTCGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                       >  1:581536/1‑251 (MQ=255)
                                                                                                                             gatgaGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                       >  2:169732/1‑251 (MQ=255)
                                                                                                                             gatgaGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGGTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGCGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                       >  2:247283/1‑251 (MQ=255)
                                                                                                                                gaGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTTTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                    <  1:274349/251‑1 (MQ=255)
                                                                                                                                gaGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                    >  2:433219/1‑251 (MQ=255)
                                                                                                                                 cgTTGAAAACCTCTCTGGAAGAGGCTGTCGCGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATACGGCGTTCACGCCGCATCCGGCAtt                                                                                                   >  1:12283/2‑251 (MQ=255)
                                                                                                                                 aGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCATGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtt                                                                                                   >  1:12297/1‑251 (MQ=255)
                                                                                                                                     gAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCAcgccgtat                                                                                                                                                                                                                               <  2:9387/123‑4 (MQ=255)
                                                                                                                                     gAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCAcgccgtat                                                                                                                                                                                                                               >  1:9387/1‑120 (MQ=255)
                                                                                                                                     gAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTAc                                                                                               >  1:359738/1‑251 (MQ=255)
                                                                                                                                      aaaaCCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGCAATTGCGGAAACTAAGGCTGTTGT‑AGGCCGGTAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAAGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGGATTTTACg                                                                                              >  2:586396/1‑251 (MQ=255)
                                                                                                                                                tGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCAAGCCGCATCCGGGATTTTACGCCTTACTTCg                                                                                    >  2:335357/1‑251 (MQ=255)
                                                                                                                                                   aaGAGGTTGTCAAGATTAAACAAGGTTTGAAAAAACCGGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCaa                                                                                 <  1:335357/251‑1 (MQ=255)
                                                                                                                                                           ttCGATATTAAACAAGTTTTAAAAAAGCCTGCACGGGTGACCCGCTCCGAAATTGCGGATACTAAGGCTTTTGG‑AGTCCGAAAAAGCCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCtt                                                                          <  1:49990/249‑1 (MQ=255)
                                                                                                                                                           gacGAGATTAAACAAGGTTGGAAAAGACCTGCACGGGTGCCCCGCTACGAAATCGCGGATACTAAGGCTGTGGT‑AGGCCGAAAAAGGCGTCCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACTTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCt                                                                           <  1:564753/247‑1 (MQ=255)
                                                                                                                                                           gTCGAGATAAAACAAGGTTTGAAAAACCCTGCACGGGTGACCCGAAACGAAATGGCGGATACTGAGGGTGTTGT‑AGGCCGAAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTg                                                                         <  2:544944/251‑1 (MQ=255)
                                                                                                                                                              gagaTTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAACTTGCGGATACTAAGGCTGTAGT‑AGGCCGGAAAAGGCGTTCACGACGTCTCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTAc                                                                      >  1:48423/1‑251 (MQ=255)
                                                                                                                                                              gagaTTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTAc                                                                      >  1:236384/1‑251 (MQ=255)
                                                                                                                                                                 atcaaccAAGTTTTGAAACACCCTGCCCGGTTGCCCCGCTACGAAATTGCGGATGCTAAGG‑TTTTGTTGGGCCGAAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTg                                                                   <  1:497924/245‑1 (MQ=255)
                                                                                                                                                                          ggTTTGAAAAACCCGTCACGGGTGCCCCGCTAGGAAATTGGGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTc                                                          <  1:468156/251‑1 (MQ=255)
                                                                                                                                                                           gTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGCAAAGGCGTTCAAGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGATTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCg                                                         >  1:518260/1‑251 (MQ=255)
                                                                                                                                                                             ttGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGc                                                       >  1:194317/1‑251 (MQ=255)
                                                                                                                                                                                aaaaacCCTGCCCGGGTGACCCGCTACGAAATTGCTGATACTAATGCTGTTGT‑AGGCCGAAAAGGGCGTTCACGCCGTATCCGGCATTCGGTTCTCAATGCCTGATGCGACGCTTACCGCGGCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATATGTCTTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTcc                                                    <  1:254115/251‑1 (MQ=255)
                                                                                                                                                                                   aaaCCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAgt                                                 >  2:85377/1‑250 (MQ=255)
                                                                                                                                                                                    aaCCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAgtg                                                >  1:219658/1‑249 (MQ=255)
                                                                                                                                                                                           aCGGTTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGA‑AGGCCGAAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATTCGCTTCGGCTCCAGTGCGTCTgc                                         <  2:399047/251‑1 (MQ=255)
                                                                                                                                                                                              cgatccccGCTACAAAATTGCGGAGACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTTCGTCTGCGcc                                      <  1:187064/246‑1 (MQ=255)
                                                                                                                                                                                                    ccGCTACGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGcc                                      >  2:187064/1‑245 (MQ=255)
                                                                                                                                                                                                     cgctaCGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAgcgtg                               >  1:580402/5‑251 (MQ=255)
                                                                                                                                                                                                      gctaCGAACTTGCGGATACTAAGGCTGTTGT‑AGGCCGGCAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGCTAAGGCGTTCACGCCGCCTCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTAGGCTCCAGTTCTTCTGCGCCCAgcgtgc                              >  2:232791/4‑251 (MQ=255)
                                                                                                                                                                                                        taCGAAATTGCGGATACTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTAATTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAgcgtgcgt                            >  2:48748/2‑251 (MQ=255)
                                                                                                                                                                                                         aCGAAATGGCGGATTCTAAGGCTGTTGG‑AGGCCGAAAAGGGCGTCCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCTACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAAGACGCTTGTACTTGATACGCTTCGGCGCCAGTGCGTCTGCGCCCAgtgtgcg                             <  2:307218/249‑1 (MQ=255)
                                                                                                                                                                                                                 gcggATACTAAGGCGGTGTT‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGGTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTGTGTAGGCCGGAGAAGGCTTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGTTACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTAc                   <  1:70430/249‑1 (MQ=255)
                                                                                                                                                                                                                     atacTAAGGCTGTTGT‑AGGCCGGAAAAGGCGTTCACGCCGAATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTCCGCATTACTTCGCAATCCGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCtt               >  2:444608/4‑251 (MQ=255)
                                                                                                                                                                                                                         tAAGGCTGTTGT‑AGGCCGGAAACGGCGTCCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGCGATTCACAACGTATTGAAATTTCTTGCTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGTTACGCCTCCGCTCACGATCGTCTGAGCCCAGCGTGCAGTTCTTGTACTCTTCGt            >  2:547331/1‑250 (MQ=255)
                                                                                                                                                                                                                               tgttgt‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTcag      <  1:431006/250‑3 (MQ=255)
                                                                                                                                                                                                                                 tttgagggccGGAAAAGGCGTTCGCGCCGTATCCGGCATTCGGTGCACAATGCCTGATGGGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTcag      <  1:586396/241‑3 (MQ=255)
                                                                                                                                                                                                                                   gt‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGttcgcagtaaa  <  2:52530/250‑11 (MQ=255)
                                                                                                                                                                                                                                   gt‑AGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTaa   >  2:210169/1‑249 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
GGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGT‑AGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAA  >  NC_000913/4628468‑4628944

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