Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,618,708 G→A A41T (GCG→ACG)  nikR → DNA‑binding transcriptional repressor NikR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,618,7080GA97.2% 227.2 / ‑4.8 72A41T (GCG→ACG) nikRDNA‑binding transcriptional repressor NikR
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (33/36);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCGTCGAAACCCAGGTGGTGGGAGAGAAATTAACCTTTTCCTCTGACGCCGGACGTGTGCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGGCGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGC  >  NC_000913/3618461‑3618954
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
aaaTCGTCGAAACCCAGGTGGTGGGAGAGAAATTAACCTTTTCCTCTGACGCCGGACGTGTGCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGacgc                                                                                                                                                                                                                                                     >  1:23480/1‑251 (MQ=255)
     gtcgAAACCCAGGTGGTGGGAGAGAAATTAACCTTTTCCTCTGACGCCGGACGTGTGCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAagag                                                                                                                                                                                                                                                <  1:224106/251‑1 (MQ=255)
                ggtggtGGGAGAGAAATTAACCTTTTCCTCTGACGCCGGACGTGTGCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCcagca                                                                                                                                                                                                                                     <  1:313978/251‑1 (MQ=255)
                    gtggGAGAGAAATTAACCTTTTCCTCTGACGCCGGACGTGTGCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAgcac                                                                                                                                                                                                                                 <  1:79713/251‑1 (MQ=255)
                    gtggGAGAGAAATTAACCTTTTCCTCTGACGCCGGACGTGTGCTACAAAACGCGGTAGTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAACGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAgcac                                                                                                                                                                                                                                 <  1:416339/251‑1 (MQ=255)
                     ggggAGAGAAATTAACCGTTCTCTCTGACGCCGGACGTGGTCTACAAAACGCGGTATTACCCGCATTCCCCGTGCACCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGGGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTGTCGTCGACATTCTGCGTAGCGCCCGGACGCAAGAGGCCACTcagcatcatg                                                                                                                                                                                                                                <  2:174814/250‑6 (MQ=255)
                               aTTAACCTTTTCCTCTGACGCCGGACGTGTGCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAgg                                                                                                                                                                                                                      >  2:233500/1‑251 (MQ=255)
                               aTTAACCTTTTCCTCTGACGCCGGACATGTGCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAgg                                                                                                                                                                                                                      >  1:320539/1‑251 (MQ=255)
                                 tAACCTTTTCCTCTGACGCCGGACGTGTGCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGtt                                                                                                                                                                                                                    >  1:149392/1‑251 (MQ=255)
                                    ccTTTTCCTCGGACGCCGGACGTGTGCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATGTACTGGAGACGCTGGACAGCCTGAGCCGGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCTCGCAAGGTTTcg                                                                                                                                                                                                                 <  2:193865/251‑1 (MQ=255)
                                          ccTCTGACGCCGGACGTGTGCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCAACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGc                                                                                                                                                                                                           >  2:82722/1‑251 (MQ=255)
                                                  gCCGGACGTGTGCTACAAAACGCGGTATTACCCGCATTGCCCGTGCGCCGGCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGATGCTGGACAGCCTGGGCCAGCGTTGTGGTTATAACAACCTTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCGTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTAt                                                                                                                                                                                                   <  1:268328/251‑1 (MQ=255)
                                                    cGGACGTGTGCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGACGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCCGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTAtgt                                                                                                                                                                                                 >  2:16864/1‑251 (MQ=255)
                                                     ggACGTGTGCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAAAGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTAtgtg                                                                                                                                                                                                >  1:279831/1‑251 (MQ=255)
                                                      gACGTGTGCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATgtgt                                                                                                                                                                                               >  2:344697/1‑251 (MQ=255)
                                                      gACGAGTGCTACAAGCAGCGGTATTACCCGCGTTCCCCGTGCGCCGTCGCACCACAGATAAGGTTTAACGCCGATGCAACGAGTCACCATCGCGCTTGCTGACGATTTACTGGAGACGCTGGACAGCCTGAGACAGCGTCGTGGTTATAACAACCGTTCCAAAGCTATCCGCGCCATTCTGCGTAGCACCCTGACGCAAGAGGCCACGCAGCAGCACGGCACGCAAGGTTTCGCGGTGCGGTCGTATgtgt                                                                                                                                                                                               >  2:352685/1‑251 (MQ=255)
                                                         gtgcgCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCTCACAACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTCCTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGCCATTCTGCGTAGCGCCCTGCAGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATg                                                                                                                                                                                            >  1:292099/5‑251 (MQ=255)
                                                             gCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAac                                                                                                                                                                                         >  2:127318/1‑250 (MQ=255)
                                                             gCTACAAAACGCGGTATGACCTGCATTCCCCGTGCGCCGTCGAACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCTCCAGCACGCTGGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACGGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGATGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAaca                                                                                                                                                                                        <  1:252895/251‑1 (MQ=255)
                                                                  aaaaCGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGaa                                                                                                                                                                                    >  2:395277/1‑250 (MQ=255)
                                                                          gTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGGTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGAc                                                                                                                                                                           >  1:114968/1‑251 (MQ=255)
                                                                              tACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACtt                                                                                                                                                                         >  1:300291/1‑249 (MQ=255)
                                                                              tACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTa                                                                                                                                                                        <  1:334191/250‑1 (MQ=255)
                                                                                cccGCAGTCCACGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTagcc                                                                                                                                                                     <  1:5705/251‑1 (MQ=255)
                                                                                  cGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATTACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTagccag                                                                                                                                                                   <  1:457543/251‑1 (MQ=255)
                                                                                    cATTCCCCGTGCGCCGTAGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTagccagcc                                                                                                                                                                 <  1:258866/251‑1 (MQ=255)
                                                                                        ccccGTGCGCCGTCGCACCACAGAAAATGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCAt                                                                                                                                                             <  2:313604/251‑1 (MQ=255)
                                                                                                 ccGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTAGGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCAc                                                                                                                                                    <  1:43542/251‑1 (MQ=255)
                                                                                                    tCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCAccc                                                                                                                                                  >  1:192873/1‑250 (MQ=255)
                                                                                                       caccacAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAag                                                                                                                                                                                                                                                  <  2:252843/151‑1 (MQ=255)
                                                                                                       caccacAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAag                                                                                                                                                                                                                                                  >  1:252843/1‑151 (MQ=255)
                                                                                                             aGAAAAGGTTCAACGCAAATCCAACGAGTCGCCCTCACGCTTGATGACGATTTACTGGGGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTACCCGCGACATTCTGCGTAGCGCCCTGGCGCAAGAGGCCGCCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACGCGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCGCCCAGCAtcatca                                                                                                                                        <  2:251942/251‑1 (MQ=255)
                                                                                                               aaaaGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATcacc                                                                                                                                      <  2:360039/251‑1 (MQ=255)
                                                                                                                 aaGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACACCCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGGCGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATcaccac                                                                                                                                    <  2:120420/251‑1 (MQ=255)
                                                                                                                        aaCGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACctc                                                                                                                              <  2:145558/250‑1 (MQ=255)
                                                                                                                         aCGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCtctc                                                                                                                            <  1:123893/251‑1 (MQ=255)
                                                                                                                         aCGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCtctc                                                                                                                            >  1:194828/1‑251 (MQ=255)
                                                                                                                            cAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACGTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGt                                                                                                                         <  2:26889/251‑1 (MQ=255)
                                                                                                                                tGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGcc                                                                                                                     >  1:240338/1‑251 (MQ=255)
                                                                                                                                 gCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGcc                                                                                                                     <  1:347496/250‑1 (MQ=255)
                                                                                                                                   aaCGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGAATCATCACCACGACCTCTCCGTCGCCACg                                                                                                                  >  2:222009/1‑251 (MQ=255)
                                                                                                                                          cACCATCACGCTTGATGACGATTTACTTGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCAtg                                                                                                           >  2:211790/1‑251 (MQ=255)
                                                                                                                                               tCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAAAACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCAACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGca                                                                                                       >  2:262110/1‑250 (MQ=255)
                                                                                                                                                   gCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTCTAACACCCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATc                                                                                                   >  2:191221/1‑250 (MQ=255)
                                                                                                                                                   gCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCAAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTGAGCCAGCCGCATTGTCTCCACCCAGAATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATc                                                                                                   >  1:99798/1‑250 (MQ=255)
                                                                                                                                                          gaCGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCacga                                                                                           <  1:200085/251‑1 (MQ=255)
                                                                                                                                                          gaCGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCacg                                                                                            <  1:127318/250‑1 (MQ=255)
                                                                                                                                                           aCGATTTACTGGAGACGCTGGACAGCCTGAGCCAGAGTCGTGGTTATAACAACCGTTCCGCAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGCGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAACAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCacga                                                                                           >  1:153982/1‑250 (MQ=255)
                                                                                                                                                            cGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCacgacg                                                                                         >  1:177696/1‑251 (MQ=255)
                                                                                                                                                             gATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCacgacg                                                                                         >  1:73635/1‑250 (MQ=255)
                                                                                                                                                                ttACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACTTCAACCACGACGACt                                                                                      <  1:122043/250‑1 (MQ=255)
                                                                                                                                                                  aCTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGt                                                                                    <  2:23480/250‑1 (MQ=255)
                                                                                                                                                                              ggACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGtctc                                                                                                                                                       <  2:152020/171‑1 (MQ=255)
                                                                                                                                                                              ggACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGtctc                                                                                                                                                       >  1:152020/1‑171 (MQ=255)
                                                                                                                                                                                   gCCTGAGCCAGCGTCGTGGTTATAACAACAGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCACCGTGTTg                                                                  >  1:68842/1‑251 (MQ=255)
                                                                                                                                                                                           cAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGACGCATTGTCTCCACCCAGCATCATCACCAAGACCTCTCAGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTg                                                           >  2:69186/1‑250 (MQ=255)
                                                                                                                                                                                                     gTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGAc                                                <  1:99277/251‑1 (MQ=255)
                                                                                                                                                                                                       tataACAACCGTTCCGAAACTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCAAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGt                                              >  1:401313/1‑251 (MQ=255)
                                                                                                                                                                                                        ataACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTTTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTg                                             <  1:401396/251‑1 (MQ=255)
                                                                                                                                                                                                        ataACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGTAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTg                                             <  1:160839/251‑1 (MQ=255)
                                                                                                                                                                                                          aacaacCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGTCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGACCACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTTCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTgc                                            <  1:60504/250‑1 (MQ=255)
                                                                                                                                                                                                          aacaacCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTgc                                            >  1:90821/1‑250 (MQ=255)
                                                                                                                                                                                                           acaacCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTgcag                                          <  2:73635/251‑1 (MQ=255)
                                                                                                                                                                                                           acaacCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTgcag                                          <  2:320539/251‑1 (MQ=255)
                                                                                                                                                                                                              acCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCAt                                       >  1:184091/1‑251 (MQ=255)
                                                                                                                                                                                                                     cGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCg                                <  1:298776/251‑1 (MQ=255)
                                                                                                                                                                                                                      gAAGCTATCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCg                                >  1:286563/1‑250 (MQ=255)
                                                                                                                                                                                                                            aTCCGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGAc                           <  1:222009/249‑1 (MQ=255)
                                                                                                                                                                                                                              ccGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCTATGACGTTa                       <  2:435171/251‑1 (MQ=255)
                                                                                                                                                                                                                              ccGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTa                       <  1:395277/251‑1 (MQ=255)
                                                                                                                                                                                                                              ccGCGACATTCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGGTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCGCCACCCAGCAGCATCACCACGACCTCTCCGTCGCCACGCTGCAGGTGCACATCAACCACGACGACTGTCGCGAAGACGCCGTGGTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTa                       <  1:290974/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tCTGCGTAGCGCCCTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAgc              >  2:4015/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ccTGACGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACCTCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGGGCAGCATTTTGCCGATTACGTTATCTCCCAGCGCGGCGTGCGGc  >  2:11775/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
AAATCGTCGAAACCCAGGTGGTGGGAGAGAAATTAACCTTTTCCTCTGACGCCGGACGTGTGCTACAAAACGCGGTATTACCCGCATTCCCCGTGCGCCGTCGCACCACAGAAAAGGTTTAACGCAAATGCAACGAGTCACCATCACGCTTGATGACGATTTACTGGAGACGCTGGACAGCCTGAGCCAGCGTCGTGGTTATAACAACCGTTCCGAAGCTATCCGCGACATTCTGCGTAGCGCCCTGGCGCAAGAGGCCACCCAGCAGCACGGCACGCAAGGTTTCGCGGTGCTGTCGTATGTGTATGAACACGAAAAACGCGACTTAGCCAGCCGCATTGTCTCCACCCAGCATCATCACCACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGC  >  NC_000913/3618461‑3618954

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: