Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,153,595 A→C 7.7% intergenic (‑497/+101) ppc ← / ← argE phosphoenolpyruvate carboxylase/acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,153,5950AC 7.7% 189.6 / 9.5 65intergenic (‑497/+101)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base C (3/2);  total (33/32)
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.26e-01

AAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTT  >  NC_000913/4153311‑4153879
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                             
aaaTGTTAAAAAATCGCCCCAAGGAACACAAAAGGTGTGGGTCGGATAAGATGCGAAAGTAGCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGGTGGGCATCGTTTGAGACACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGt                                                                                                                                                                                                                                                                                            <  2:400254/288‑1 (MQ=255)
    gTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGAt                                                                                                                                                                                                                                                                                       >  1:388860/1‑289 (MQ=255)
       aaaaaaTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAg                                                                                                                                                                                                                                                                                    <  2:388860/289‑1 (MQ=255)
                  caaaGTAAAACCACAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGCCATGATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgca                                                                                                                                                                                                                                                                         <  2:522112/287‑1 (MQ=255)
                        aaCACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatcgca                                                                                                                                                                                                                                                                    <  1:276943/288‑1 (MQ=255)
                           acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCGTCGTTCGAGTCACTGTAACTCGGATACGATGCGCAAGCATCTCAACCg                                                                                                                                                                                                                                                                >  2:416273/1‑289 (MQ=255)
                                             aTAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACT‑GTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTTAGTCACTGTAAGTCGGCTAAGATGCGCAAGCATCGCATCCGACCTTATTGCGGCACTgg                                                                                                                                                                                                                                              >  2:388843/1‑289 (MQ=255)
                                                       cAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCAGTGTAGGTCCGATAAGATTCGTAAGCATCGCATCCGGCATTATTGGGGCACTGGAATTTggta                                                                                                                                                                                                                                     >  2:290469/1‑286 (MQ=255)
                                                                          gattGCGGCACGGGAGTTTGGCAACGGTGCCGGAGGGGGCGCGTGCGCCTTATCCGGCGTAAAGTTGGGAATCGGTTGAGTCGCT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGcgcg                                                                                                                                                                                                                 <  2:136804/288‑1 (MQ=17)
                                                                                 gCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGCTGCGgcgccgcggcct                                                                                                                                                                                                          >  2:153272/1‑281 (MQ=255)
                                                                                 gCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGACTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGAGAAGATGCGCCAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCCACAGGGCCGGATGCGGCTCGAGCGCCt                                                                                                                                                                                                          >  2:211905/1‑289 (MQ=255)
                                                                                     gggAGTTTGGCAACAGTGCCGAATGGGGCGCGGGCGCCTTATCCGGCCTAAAGTTGGGCATCGTTTGGGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                                                                                                      <  2:408801/288‑1 (MQ=255)
                                                                                            tGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCCTCGTTTGAGTCACTGTCGGCCGGATAAGATGCGCAAGTATCGCATCCGTCATTATTGCCGCACTGGAGTTTGGCGACAGTGCCGGATGCGGTGCGAGCGCCTTCTCCGGCCt                                                                                                                                                                                                >  2:290329/1‑288 (MQ=17)
                                                                                                    gTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGGTGCGGCGCGAGCGCCTTAGCCGGCCTTCGCTTTaa                                                                                                                                                                                       >  2:41991/1‑289 (MQ=255)
                                                                                                                ggcggAGGGCCTTATCCGGCCTACGTTTGGGAATCGTTTGAGT‑ACTGGTAGGTCGGATAAGATGGGCAAGTATCGATTCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCGGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCaaaa                                                                                                                                                                            <  1:388843/284‑1 (MQ=17)
                                                                                                                 cgcgAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACCGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGACTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCCAGCATCGCATCCGACATTATTGCCGCACTGGCGCTTGGCAACAGCCCCGGACGCGGCGCGAGCCCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAAt                                                                                                                                                                          >  1:96677/1‑289 (MQ=255)
                                                                                                                      gcgcCTGTTCCGGCCAACAGTTGGCCATGTTTTGAGTCACT‑GTCGGTCGGTTAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTg                                                                                                                                                                      <  1:66950/288‑1 (MQ=255)
                                                                                                                       ctccTTATCCGGCTTACAGTTGGGCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGa                                                                                                                                                                     <  1:60205/286‑1 (MQ=255)
                                                                                                                            tatcagagttaaaGTTGGGCATCGTTTGAGTCATT‑GTCGGTCGGTTAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTAcc                                                                                                                                                               <  2:500708/277‑1 (MQ=255)
                                                                                                                            tATCCGGCCTACAGTTGGGCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGGCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGAAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTAcc                                                                                                                                                               <  1:153272/289‑1 (MQ=255)
                                                                                                                                    ttACGGTTGGGCATCGTTTGAGTCACT‑GTGGGTCGGATAAGATGCGCAAGTATCGCATCCGAAATTATTGCGGCACTGGGGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTAt                                                                                                                                                       <  1:290329/288‑1 (MQ=255)
                                                                                                                                             ggCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCCCTGGAGTTTGGCAACAGCGCCGGCTGCGGGGCGAGCCCCTTATCCGGCCTACGTTTTACTGCCACCAAAAATGGTGAATTAACTGGGTTATCAGttaggg                                                                                                                                              >  1:399955/1‑285 (MQ=255)
                                                                                                                                                 tCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatcg                                                                                                                                                                                                                                                                      <  2:453330/165‑1 (MQ=255)
                                                                                                                                                 tCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatcg                                                                                                                                                                                                                                                                      >  1:453330/1‑165 (MQ=255)
                                                                                                                                                    tttGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGCGTCCCTGTAGGTCGGATAAGATGCGCCAGCATCGCATCCGACATTATTGCGGCCCTGGAGATAGGCCACCGCGCCGCATGGGGCGCGAGCCCCTTATCCGGCATACCTTTTACGGCCAGCAAAAATGGTGAATTACCTGGGTCTTCAGGTCGCCGGGTg                                                                                                                                         >  2:415790/1‑287 (MQ=17)
                                                                                                                                                       gAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGCTGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTACTGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑TCGGTTGGCTTg                                                                                                                                   >  1:167810/1‑289 (MQ=255)
                                                                                                                                                          aaggt‑GTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATa                                                                                                                                >  2:111795/5‑289 (MQ=255)
                                                                                                                                                              t‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATgcgc                                                                                                                                                                                                                                                                              <  2:206910/144‑1 (MQ=37)
                                                                                                                                                              t‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCCACCGCGCCGGATGCGGCGCGCGCGCCCTATCCGGGCGACGTTTTAATTCCAGCACAAATGGTTAATTATCCGGGTTATCAGTTTG‑CGCGTGGGCTTgttcaat                                                                                                                             >  2:429295/1‑282 (MQ=17)
                                                                                                                                                                    gtcggAAAAGATGCGCATGTATCGCATCCGACAGTATTGCGGAACTGGAGTGTGGCAACAGTGCCGGAGGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCgtgt                                                                                                                        <  1:186084/288‑1 (MQ=255)
                                                                                                                                                                       ggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTcc                                                                                                                    <  1:111795/289‑1 (MQ=255)
                                                                                                                                                                       ggATAAGATGCGAAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTcc                                                                                                                    <  1:290469/289‑1 (MQ=255)
                                                                                                                                                                         aTAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatcgca                                                                                                                                                                                                                                                                    >  1:206910/1‑144 (MQ=37)
                                                                                                                                                                           aaGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGat                                                                                                                <  1:68157/289‑1 (MQ=255)
                                                                                                                                                                                   cAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGGCTTTATTTCGGCCCCGGAGTTTGGCAACAGCGCTGGATGCGGCGTCAGCGCCTTATCCGGCCTACGTTTTAATGCCCTCAACAATGGATGCTTACCTGGGTTATCCATTTG‑CGGGTTAGCTTGTCACACCGTGCTTCCAAATATtcatca                                                                                                        >  2:297357/1‑289 (MQ=14)
                                                                                                                                                                                   cAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTTAATTACCTGGTTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATtcatcc                                                                                                        >  2:81144/1‑288 (MQ=255)
                                                                                                                                                                                   aaagcagcGCATCCGACTATATGCCGTCACTGGAGTTTGGCAACAATGCCGGATGCGGCGCGAGCGCCATAGCCGGCCTACGGTTGGGCATCTGTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATtcatc                                                                                                         <  2:302675/281‑1 (MQ=21)
                                                                                                                                                                                       caTCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                                                                                                      <  2:270842/191‑1 (MQ=37)
                                                                                                                                                                                              tcggaaattgttGCGGCACGGGAGTTTGGCAACAGGGCCGGAGGCGGCGCGAGTGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCAACGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGc                                                                                             <  2:96677/283‑1 (MQ=255)
                                                                                                                                                                                              tCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGGTGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAg                                                                                              >  1:460659/1‑288 (MQ=255)
                                                                                                                                                                                              tCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTAcc                                                                                                                                                               >  1:190252/1‑224 (MQ=255)
                                                                                                                                                                                              tCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTAcc                                                                                                                                                               <  2:190252/224‑1 (MQ=255)
                                                                                                                                                                                                 gACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGt                                                                                                                                                                                            >  1:270842/1‑192 (MQ=37)
                                                                                                                                                                                                   cATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGa                                                                                         <  1:41991/288‑1 (MQ=255)
                                                                                                                                                                                                       atggCGGAACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCAAGCGCCTTATCCGGCTTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGAATAAGATGCGCAAGCATCGCATCCGAAATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAAt                                                                                    <  1:464772/286‑1 (MQ=255)
                                                                                                                                                                                                           cGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGACGCGGCGAGAGAGCCTTATCCGGCCTACGTTTTATTGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAGACCGTGCTTCCCGATATTCATTAGGTTGATGATCCTGATGTATTGa                                                                                 >  1:416855/1‑288 (MQ=21)
                                                                                                                                                                                                             gCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTAt                                                                                                                                                       >  1:15327/1‑217 (MQ=255)
                                                                                                                                                                                                             gCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTAt                                                                                                                                                       <  2:15327/217‑1 (MQ=255)
                                                                                                                                                                                                             gCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGc                                                                               >  2:101517/1‑288 (MQ=255)
                                                                                                                                                                                                                     gTTTGGAAACAGGGCCGGATGCGGGGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCAGCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGcccc                                                                       <  1:416273/288‑1 (MQ=255)
                                                                                                                                                                                                                             aCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGc                                                                                             >  2:378421/1‑258 (MQ=255)
                                                                                                                                                                                                                             aCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGc                                                                                             <  1:378421/258‑1 (MQ=255)
                                                                                                                                                                                                                               aaggCCGGGTGCGGAGCGAGCGCCGTAGCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATGATTGCGGAACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGt                                                            <  2:169875/286‑1 (MQ=255)
                                                                                                                                                                                                                                    cGGATGGGGCGCGAGCGCCTTAGCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGGGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCggg                                                        <  1:211905/288‑1 (MQ=255)
                                                                                                                                                                                                                                                       ttATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACg                                     >  2:427793/1‑288 (MQ=255)
                                                                                                                                                                                                                                                            cgtccTACGGTTAGGATGCGTTTGAGTCACGGTGGGTCGGATAAGATGCGCAAGCATTGGATCCGACATTCTTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGAGCCTTATCCGGCCTACGTTCTAATGCCAGCAAAAATGGTGAATTACCTGGTTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCtt                               <  2:115986/286‑1 (MQ=255)
                                                                                                                                                                                                                                                            cGGCCTACGGTTGGGCATCGTTTGATGCACTGTAGGTCGGATAAGATGCGAAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCtt                               <  1:414359/289‑1 (MQ=255)
                                                                                                                                                                                                                                                              gCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGAATAACGTTTGAATAAACGGCGCTTCg                             >  1:16851/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                        gggCATCGTTTGAGTCACTGTACGACGGATAATATGCGCAAGCATCGCATCCGACATCATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGACTTATCCGGCCTACGCTTTAATTTCCGCTCAACTGGCGAATTATCTGGGCTATCAGTTCG‑CGGTTGGTCTTGATAAACCGTGTTTTCAGATATTCATCACGTTGATTAGGCTGATTATTTGAGACAGGCCCCAACACCATCGACGCACATAACCTTTGAATATCAGGCGCTTCAGTACagtcgt                   >  2:74026/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                         ggCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCCCTTCGGTACagtagt                   >  2:237510/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                             tCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTcacc              >  1:288983/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                              cGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTcacc              >  2:468938/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                   gAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGt                                                            >  2:355032/1‑237 (MQ=255)
                                                                                                                                                                                                                                                                                   gAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGt                                                            <  1:355032/237‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          tGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGttttt  >  1:465890/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                             
AAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACT‑GTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCG‑CGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTT  >  NC_000913/4153311‑4153879

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑