Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,587,112 C→T A461V (GCA→GTA)  narQ → sensor histidine kinase NarQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,587,1120CT89.5% 101.9 / NA 39A461V (GCA→GTA) narQsensor histidine kinase NarQ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (21/13);  minor base G (4/0);  total (26/13)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCAC  >  NC_000913/2586869‑2587357
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
tacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGa                                                                                                                                                                                                                                                 <  1:53940/250‑1 (MQ=255)
       ccAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCaa                                                                                                                                                                                                                                          >  1:161941/1‑250 (MQ=255)
           ttGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGc                                                                                                                                                                                                                                      >  1:115007/1‑250 (MQ=255)
                       aaGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGAGTTGAATGATGCTTATCGGCAGTTAAGCGAGCTGTTGACTACTTTCCGCCTGACACTGCCGCAGGCGGATCTCCCCTCCGCCTTGAGGGAAATGCTGGGTACGGTACACAATCAAACCAGCGCCAAACAGACCCCCGCATGCCGGCTGCCCACCCTGGGAATGGCTGCGCCAAAGCAGGTGCCTttgtt                                                                                                                                                                                                                         >  2:342226/1‑251 (MQ=255)
                          cGTTCGATACCGGAAGATAACGCCACCGCCCAACGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCCGTTACGCGAGCTGTTGACTACTTTCCCCCTGACGCTGCCGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGGTACGTTGCATAATCAAACCCGCACCAAACTGAACCTCGACTGCCCTCTGCCCACCCCGGGACTGGAGCCGAAAATGCCGGGGCGTTTGTTGcc                                                                                                                                                                                                                      >  2:107705/1‑250 (MQ=255)
                          agTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATATCCCCTCCGCATTGGGGGAAATGCTGGGTACGTTACAAAATCAAACCAACGCCCAAATGACCCTCGACTGACGTCTGCCAACCCTGGGACTGGATGCGACATTGCAGGGGCATTTGTTGc                                                                                                                                                                                                                       >  2:359247/2‑250 (MQ=255)
                                  aCCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTc                                                                                                                                                                                                              <  2:349945/251‑1 (MQ=255)
                                           tAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCgg                                                                                                                                                                                                     >  1:205773/1‑251 (MQ=255)
                                                 cACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGACATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACt                                                                                                                                                                                                                                                    >  2:58915/1‑198 (MQ=255)
                                                 cACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACt                                                                                                                                                                                                                                                    <  1:58915/198‑1 (MQ=255)
                                                               tCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGATACGCGAGCTGTTGACTACTTTCCGGCTGACGCGGCAGCAGGCGGATCTACCCTCCGCATTGAGGGAAATGCTGGATAAGTTACAAAATCAAACCAGTGCCAAAATTAACCTCGACTGCCGTCTGTCAACCCTGGTACTGGAGACGCCAATGCAGGTGCATTTGTTGCAAAATATTCGCGAAGCGGTGGTGCATGCGGTGAAGCa                                                                                                                                                                                  >  1:113312/1‑250 (MQ=255)
                                                                cATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacg                                                                                                                                                                                <  1:54072/251‑1 (MQ=255)
                                                                                 aGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTTCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAAAGCGATGAAGCACGCCAACGCCAGCGAAATc                                                                                                                                                               >  2:388858/1‑251 (MQ=255)
                                                                                       tGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTc                                                                                                                                                         >  1:185922/1‑251 (MQ=255)
                                                                                             atgCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGt                                                                                                                                                                                                                                  >  1:127256/1‑172 (MQ=255)
                                                                                             atgCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGt                                                                                                                                                                                                                                  <  2:127256/172‑1 (MQ=255)
                                                                                                   aTCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAattatt                                                                                                                                                                                                               <  1:70538/185‑1 (MQ=255)
                                                                                                   aTCGGCAGGTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGGGGATCTCCCCTCCGCATTGTGGGAAATGCTGGGTACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGGTATGGATGCGCAAATGCAGGGGCATTTGTTGCAAATTATTCCGGAAGCGGTGCTGAAAGAGAGGAAGCACGCCAACGCCCGCGAAAACGCCGGCAGTTGCGTCAcc                                                                                                                                             >  2:163401/1‑251 (MQ=255)
                                                                                                     cGGCAGTTACGCGTGCTGTTGACTACTTTCCGCCTGCCGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGTAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTCATGGATGCGCAACTGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGATGCTGAATGCGATGAAGCACGCCAACGCAAGAGAAATCGCCGTAAGTTGCGTCAcggc                                                                                                                                           >  2:233996/1‑248 (MQ=255)
                                                                                                     cGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACcgc                                                                                                                                           <  2:205773/251‑1 (MQ=255)
                                                                                                                      gTTGACTCCTTTCCGCCGGATGCGTGAGCAGGGGGATCTCCCCTCCGCATTGATGGAAATGCTGGTGCAGTTAAAAAATCAAACCAGCGCCAAATTGACCCGCGCCCGCCGTCTGCCAACCCTGGTACGGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGTGGTGCTGAATGCGATGAAGCACGTCAATCCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATca                                                                                                                            <  2:47196/249‑1 (MQ=255)
                                                                                                                            tactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCAGTCTGCCAACCCTGGTACTGGATGCGCCAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGttt                                                                                                                    >  1:82803/1‑251 (MQ=255)
                                                                                                                             actTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGttt                                                                                                                    <  1:448050/250‑1 (MQ=255)
                                                                                                                                    gCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGt                                                                                                            >  2:187417/1‑251 (MQ=255)
                                                                                                                                      cTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGAGAAATGCTGGATACGTTACAAAATCAAACCAGCGCAAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCTCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGACGCGGTGCTGAATGCGATGCAGCCCGCCATAGCCAGCGCAATCGCCGGACGTTGCGTCCCCGCGCCGGACGGCAATCACACGGTTTATATCCGTg                                                                                                           >  2:298281/1‑250 (MQ=255)
                                                                                                                                         aCGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATaa                                                                                                       <  1:302711/251‑1 (MQ=255)
                                                                                                                                                                      cATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGCCGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATAGGTGCACCGAAAGa                                                                            >  2:466550/1‑250 (MQ=255)
                                                                                                                                                                      cATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTAGACAGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACAGCGCCGGACGGCAATCACACGGTTTATATCCGTGCTAACGGGATTGGTATCGGTGAACCGAAAGa                                                                            >  1:475522/1‑250 (MQ=255)
                                                                                                                                                                                      ggATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATg                                                           <  1:113640/251‑1 (MQ=255)
                                                                                                                                                                                          aCGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCt                                                       <  1:388858/251‑1 (MQ=255)
                                                                                                                                                                                             gTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAGGGTCATTATGGTCTg                                                      >  1:127593/1‑250 (MQ=255)
                                                                                                                                                                                                aCAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAat                                                   >  1:328612/1‑250 (MQ=255)
                                                                                                                                                                                                                     aCTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGACATTATGGTCTGCATTTCATGCGCGCACGCGCGGaa                              >  2:431724/1‑250 (MQ=255)
                                                                                                                                                                                                                         aCCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCt                         >  2:171792/1‑251 (MQ=255)
                                                                                                                                                                                                                           ccTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCt                         >  2:128074/1‑249 (MQ=255)
                                                                                                                                                                                                                              cGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAgg                      <  1:171792/249‑1 (MQ=255)
                                                                                                                                                                                                                                   gCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTgggccg               >  2:132506/1‑248 (MQ=255)
                                                                                                                                                                                                                                         tGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACt         >  2:397754/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ccTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTcac  >  2:400618/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
TACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCAC  >  NC_000913/2586869‑2587357

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