Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,303,243 A→G 14.5% S85S (TCT→TCC napF ← ferredoxin‑type protein, role in electron transfer to periplasmic nitrate reductase NapA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,303,2430AG14.5% 162.7 / 2.4 70S85S (TCT→TCCnapFferredoxin‑type protein, role in electron transfer to periplasmic nitrate reductase NapA
Reads supporting (aligned to +/- strand):  ref base A (20/37);  new base G (9/3);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 2.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGAGATACTCCGCAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGG  >  NC_000913/2303013‑2303480
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 
tGGAGATCCTCCGGAGTTATGGCTGATACCGGACAGCTGGCTGAACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAg                                                                                                                                                                                                                                          <  1:110317/239‑1 (MQ=255)
tGGAGATACTCCGCAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACCCCAGGCCCCCGTCCCCCATGGCGAACTGCCAATCCCAGGCCCTGGGGTGGCGCGGGGgaaaat                                                                                                                                                                                                                                           >  1:195014/1‑236 (MQ=255)
    gaTACTCCGCAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACCGACTGATACCCCAGGGACCGGTCCCCGATGGTGAACTGCCAACCCCAGGCCCCTGGGTGGCGCGGCGCAAATAACGAt                                                                                                                                                                                                                                      >  1:166785/1‑239 (MQ=255)
      tACTCCCCAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAGAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGAtt                                                                                                                                                                                                                                     <  1:143460/238‑1 (MQ=255)
      tACCCCGCAGTTAGGGTGGTTCCCGGACGGCTGGCGGCCTTAGCTCCGACGCCGGTGCAAAGTTGGCTATTAATGTGCGGCTGGTAGATCCTGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCGGGCACTGGTGTGGCGCGGAGAAAATAGCGAtt                                                                                                                                                                                                                                     <  2:517099/238‑1 (MQ=255)
        cTCCGCAGTTATGGCTGATACCGGACAGCTGGCTTCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATACCTGAGAGCGTCCGGCGCAAAATAAATGCCATTTGCTCCCAGCTATCCTGGCAGCGGCGGGCATCCACCGGCTTGTAACCCCGGCACGCGTTCCCCATGGTGAACTGCACATCCCCGGCCCCGGGTGTGGCGGGGGCaaaaaacggatgcg                                                                                                                                                                                                                                   >  2:508492/1‑229 (MQ=255)
         tCCGCAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTCGGGTGGGCCGGGGAAAAAACGGCTTCggg                                                                                                                                                                                                                                 >  1:378055/1‑239 (MQ=255)
           cGCAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGCCAGCGGCGGCATTCAACTGACTGATACGCCAGGCACACCTCCCCGATGGTGAACTGCAAATCCCCGACCCTGGTTTGGGGCGGAGAAAATAAGGGTTCggggt                                                                                                                                                                                                                               >  1:117598/1‑238 (MQ=255)
             cAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATTCCGGAGAGCCGCGGGCGGAAGATACTTGCCCTTGGCTCACAGCCATCCTGACAGCGGCGCGATTTAAATGACGCATACGCCAGGCCACCGTCCCCCCCGGGTAACGCCAAAACCCCGGACCTGGTGGGGGCGGGGGGCAAAAGCGAGTCgggggca                                                                                                                                                                                                                             >  2:687935/1‑236 (MQ=255)
                  aTGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCCATCCTGACAGCGGGGGGCTTCAAATGACAGGTACGCCCGGGACCCCTCCCCCAGGGGGAAACGCCAATCCACGGCCCCGGTTGGGGGCGGGGGAAAAAGGGAATCGGGGGaaac                                                                                                                                                                                                                           >  2:519615/1‑234 (MQ=255)
                        gATTCCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTgcgcgc                                                                                                                                                                                                                   <  1:832265/238‑1 (MQ=255)
                        gATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTgcgcgc                                                                                                                                                                                                                   <  1:835611/238‑1 (MQ=255)
                               gACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGGGTGGCGCGGAGCAAATAGCGGTTCGGGGCAAGCCTGCGCCCAGGCGGaa                                                                                                                                                                                                           >  1:903733/1‑239 (MQ=255)
                                    ctggctgCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGGGAAAATAGCGATTCGGGGCAAGCCTGCCCGCAGACGTAACAGaa                                                                                                                                                                                                      >  1:853008/1‑239 (MQ=255)
                                    ctggctgCACATGCTCCACATCCGCTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGGTAACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCCCGCGTCCCCGATGGGGAACTGCAAAACCCAAGCCCCGGGGTGGCGCGGGGAAAAAAACGGATCTGGGGACGCCCGCGCGCggggggaaacaaa                                                                                                                                                                                                      >  2:886012/1‑229 (MQ=255)
                                      ggttggccaTGCGCCCCATTCGTTGGAAGGTTGGCTTATAAGTGGCGGTGGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCAGTGGCTCACAGCTAGCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGc                                                                                                                                                                                                    <  2:68497/232‑1 (MQ=255)
                                              aTGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACACCTATCCTGACAGCGGGGGCCTTCAACTGACTTATAACCCATGCCACCGTCCCCGCTGGTGAACTGCAAAACCCCAGCCCCGGCTGTGCCCGGGGCAAAATCGCGTTTGGGGCAAACCCGCGCGCAGGGGTAACACAAGccggacacg                                                                                                                                                                                            >  2:669711/1‑231 (MQ=255)
                                                     acaTCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACCCGTACtttt                                                                                                                                                                                     >  1:32815/1‑239 (MQ=255)
                                                      caTCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGGGAACATAGCGCTTCGGGGCAAGCCTGCGCGCAGGCGTAACCGAAACTGCAATCGTTAttttt                                                                                                                                                                                    >  1:434838/1‑239 (MQ=255)
                                                            ttGCAAAGTGTGCTATTGAGTTGCGGCTGGTAAATCCCGGAGAGCTGCGGGCGAAAGATAATTGCCATTGGCTCACAGCCATCCTGACAGCGGCGGCATTCAACTTACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAAtt                                                                                                                                                                              <  2:34418/239‑1 (MQ=255)
                                                                 agtgtGGCTATTAAGTGGCGGGTGGAACATCCCGGGGAGCGTCGGGCGAAAGAAAATTGCAGTTGGCTACCAGGTACCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCCCGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGc                                                                                                                                                                       <  2:32815/235‑1 (MQ=255)
                                                                   tttGGGTTTTTAGTGGCGGGTGGTAGATCCCGGGGAGGGTGGGGCGAAAAAGAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACCGACTGAAACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCGGGCCCTGGTGTGGGGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTc                                                                                                                                                                     <  2:378055/238‑1 (MQ=255)
                                                                       gCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCCCTCGTTATTTTTGAAATTAA‑‑CGCTCGGaa                                                                                                                                                                 >  1:767991/1‑238 (MQ=255)
                                                                             aaGTTGCGGCTGTTAGATCCCGGAGGGCGTGGGGCGAAAGATAATTGCCATGGGCCCACAGCTATCCGGACAGCGGCGGCCTTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAgccgc                                                                                                                                                           <  2:903733/239‑1 (MQ=255)
                                                                                     gCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGGGCACTTCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGCTTCGGGGCAAGCCTGCGCGCCGGCGTGACAGAAACTGCACTCGTTATTTTTTAAATTAA‑‑CGCCCGGATAACCTCccgcgccc                                                                                                                                                   >  2:98723/1‑238 (MQ=255)
                                                                                            agtCCCGGAGACCGTCGGGCGAAAGTTAATTGCTATTGGCTCACAGCTCTCCTGACAGCGGCGGCATTCAACTGACTGATACGGCAGGCACGGGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCa                                                                                                                                            <  2:166785/237‑1 (MQ=255)
                                                                                                        cGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTCCAGAATGTTGttt                                                                                                                                 >  1:288956/1‑238 (MQ=255)
                                                                                                        cGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGGGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGttt                                                                                                                                 >  1:308166/1‑238 (MQ=255)
                                                                                                           cGGGCGAAAGATAATTGCCTTTGGCTCACAGCTATCCGGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGGGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTca                                                                                                                              <  2:279962/238‑1 (MQ=255)
                                                                                                           cGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTca                                                                                                                              <  1:519615/238‑1 (MQ=255)
                                                                                                            ggggGAAAAATAATTGCCATGGGCCAACAGGTATCCGGCAAGCGGGGGCATTAAACTGACGGAGCCGCCGGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGGGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTcaca                                                                                                                            <  2:6872/239‑1 (MQ=255)
                                                                                                            gggCGAAAAATTATTGCCTTTGGCTCCCAGGTTCCCTGGCAGGGGCGGCATTCAACTGACTGATCCTCCCGGCAGGGGGCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGGGGGGAGGAGAAAATAGCGATTCGGGGCAAGACTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTcaca                                                                                                                            <  2:541953/239‑1 (MQ=255)
                                                                                                             ggCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTcaca                                                                                                                            <  1:111253/238‑1 (MQ=255)
                                                                                                             ggCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTcaca                                                                                                                            <  1:846262/238‑1 (MQ=255)
                                                                                                             ggCGAAAGATAAGTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTcaca                                                                                                                            <  1:277876/238‑1 (MQ=255)
                                                                                                                 aaaGATAATTTCCATTGGCCCCCCGCTATCCTGACGGGGGCGGCATTTAACTGACTTATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCAt                                                                                                                       <  2:770288/239‑1 (MQ=255)
                                                                                                                  aaGATAAGTGCCATTGGCTCAAAGCTATCCTGACAGGGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCAt                                                                                                                       <  1:99701/238‑1 (MQ=255)
                                                                                                                            ccATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACTCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCCCGCAGGCGTAACAGAAGCTGCACTCTGTATTTTTGGAATTAA‑‑CGCTCGGATAGCCCCCCCAGCCACGATCCCGAATGGTGTTTTTACAACCAtacaaaaagcgc                                                                                                            >  2:587453/1‑228 (MQ=255)
                                                                                                                                                tgcaggcggcgGCTTTAAACGTACTGAGCCGCCAGGCACGCGTCCCCGATGGTGAACTGCAAAACCCAGGCCCGGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGt                                                                                       <  2:434838/234‑1 (MQ=255)
                                                                                                                                                cgACAGCGGCGGCAATCAACTGACTGATACGCAAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGt                                                                                       <  1:669711/238‑1 (MQ=255)
                                                                                                                                                  aCAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGga                                                                                                                                                                                                                                                   <  1:873745/85‑1 (MQ=255)
                                                                                                                                                  aCAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGga                                                                                                                                                                                                                                                   >  2:873745/1‑85 (MQ=255)
                                                                                                                                                  aCAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATACCGATTCGGGGCAAGCCTGCGCGCCGGCGCAACAGAAGTTGACTTCGTTATTTTTGCACTTAA‑‑CGCCCGGATAGCCGCCCCCGCCGCATTGCAGAATGGTGGTTTCACAAGCATAAATAAAGCCCTCCCACCGGGGGAAAgggggaa                                                                                     >  2:692273/1‑236 (MQ=255)
                                                                                                                                                      cggcggCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGaaa                                                                                 >  1:756085/1‑239 (MQ=255)
                                                                                                                                                      cggcggCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGCCGCAGCGGGTGCAATTGGTCAGaaa                                                                                 >  1:541504/1‑239 (MQ=255)
                                                                                                                                                              ttCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGAtt                                                                         <  1:3023/239‑1 (MQ=255)
                                                                                                                                                               tCAACTGACTGATGCGCAAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATtc                                                                        <  1:98723/239‑1 (MQ=255)
                                                                                                                                                                    tggctgaTCCGCCAGGCCCGCTTCCCCAAGGGGGAACGGCAAATCCCAGGCCCGGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCAc                                                                   <  2:853008/236‑1 (MQ=255)
                                                                                                                                                                    tgacggaTACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCAc                                                                   <  1:22775/239‑1 (MQ=255)
                                                                                                                                                                          gtccGCCAGGCCCGCGTCCCCGAGGGTGAACTGCAAATCCCAGGCCCGGGGGTGGCGCGGGGAAAATAGCGGTTCGGGGCAAGCCTGCGCGCAGGGGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGAcc                                                             <  2:288956/236‑1 (MQ=255)
                                                                                                                                                                           tACGCCAGGCACGGGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGGGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGAcc                                                             <  1:178636/238‑1 (MQ=255)
                                                                                                                                                                           tACCCCGGGCACGGCGCCCCGATGGTGAACTGGAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGAcc                                                             <  2:308166/238‑1 (MQ=255)
                                                                                                                                                                              gCCAGGCACGGGTCCCCGGTGGTGAACGGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAggg                                                         <  1:145569/239‑1 (MQ=255)
                                                                                                                                                                                cAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAACGCGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAggg                                                         >  1:24863/1‑239 (MQ=255)
                                                                                                                                                                                  ggCCCGGGTCCACGGTGGTGAACTGAAAATCCCGGGCCCTGGGTTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGcgga                                                     <  2:756085/239‑1 (MQ=255)
                                                                                                                                                                                  ggCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGcgga                                                     <  1:886012/239‑1 (MQ=255)
                                                                                                                                                                                  ggCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCCCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACACGCATTAATACCGGCCTCGCACGCGGGGCAAAGGGGGAGAAAATGAGATTAAACAACGGACCAggggggg                                                     >  2:847033/1‑238 (MQ=255)
                                                                                                                                                                                     ccgggcccccGATGGGGAACTTCAAATCCCAGGCCCTGGGGTGGCCGGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGcggacg                                                   <  2:916735/232‑1 (MQ=255)
                                                                                                                                                                                     cacgttgcctgATGGTTAACCGAAACGCCCAGGCCCTGGTGTGGGGGGGGGAAAATGGCGATTCGTGTCAAGCCGGCGCGCAGGCGTAACGGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGAATCATCGCCGGACCAGGGcggacg                                                   <  2:117598/228‑1 (MQ=255)
                                                                                                                                                                                        cgTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGAATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGAt                                                >  1:916735/1‑238 (MQ=255)
                                                                                                                                                                                              cGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCCCAGCGGGGGCAATGGGTCAGAAAAATAGATTCATCACCGGCACAGGGCGGGCGGAAACCGtac                                         >  2:8446/1‑237 (MQ=255)
                                                                                                                                                                                                                 cccAGGCCCTGGGGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGa                       <  1:506929/239‑1 (MQ=255)
                                                                                                                                                                                                                      gCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGCTAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGCGCAATGGGTCAGAAAACGAGATTCATCACCGGACAAGGGGGGCAGGATACCGGTACGGGCTTTTGCCAAGCGACCAgg                  >  2:321830/1‑238 (MQ=255)
                                                                                                                                                                                                                       gacagggTGTGGGCGGGGAAAAATTGCGATGCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCCTCACTTGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCa                    <  2:424958/232‑1 (MQ=255)
                                                                                                                                                                                                                           tGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGGGAGt              >  1:672530/1‑239 (MQ=255)
                                                                                                                                                                                                                                gtgGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATg           <  1:692273/238‑1 (MQ=255)
                                                                                                                                                                                                                                  ggCGGGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGcc         <  1:226961/238‑1 (MQ=255)
                                                                                                                                                                                                                                   gcgcGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGcccc       >  1:44138/1‑239 (MQ=255)
                                                                                                                                                                                                                                   gcgcGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGGGCAATGGGTCAGAAAATGAGATTCATAACCGGACCAGGGCGGACGAATACCGTTACTGGCTTTGCGCCACCGACCAGGGAGGAGGcccc       >  2:226961/1‑239 (MQ=255)
                                                                                                                                                                                                                                     gcgGAGAAAATTGCTGTTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCATTCGTTATTTTTGAAATTAC‑‑ACCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGa     <  2:617152/239‑1 (MQ=255)
                                                                                                                                                                                                                                      cgGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGAc    >  1:8231/1‑239 (MQ=255)
                                                                                                                                                                                                                                      aggAGAAAATAGCCGTTGGGGGCAAGGCTGGGCGCAGGCGTAAGAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGGTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGa     <  1:172287/237‑1 (MQ=255)
                                                                                                                                                                                                                                        gagaAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGc                                                                                                                         <  2:416264/121‑1 (MQ=255)
                                                                                                                                                                                                                                        gagaAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGc                                                                                                                         >  1:416264/1‑121 (MQ=255)
                                                                                                                                                                                                                                        gagaAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACCGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGc          >  1:424958/1‑231 (MQ=255)
                                                                                                                                                                                                                                        gagaAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGCCGGATACCGTTACTGGCTTTGCGCCAGCGCCCAGTGAGTATGCCCCGACgg  >  1:84616/1‑239 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 
TGGAGATACTCCGCAGTTATGGCTGATACCGGACAGCTGGCTGCACATGCTCCACATCCGTTGCAAAGTTGGCTATTAAGTTGCGGCTGGTAGATCCCGGAGAGCGTCGGGCGAAAGATAATTGCCATTGGCTCACAGCTATCCTGACAGCGGCGGCATTCAACTGACTGATACGCCAGGCACGCGTCCCCGATGGTGAACTGCAAATCCCAGGCCCTGGTGTGGCGCGGAGAAAATAGCGATTCGGGGCAAGCCTGCGCGCAGGCGTAACAGAAGCTGCACTCGTTATTTTTGAAATTAA‑‑CGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGG  >  NC_000913/2303013‑2303480

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