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,1120CT96.1% 140.3 / NA 54A461V (GCA→GTA) narQsensor histidine kinase NarQ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (20/29);  minor base G (1/1);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTC  >  NC_000913/2586866‑2587355
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
acttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAAAGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGTACCGCGACTGCCGTATGCCAAACCTGGCaagg                                                                                                                                                                                                                                                   >  2:484495/1‑248 (MQ=255)
    acGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGGACTGGATg                                                                                                                                                                                                                                              >  1:259489/1‑251 (MQ=255)
      gTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATgcg                                                                                                                                                                                                                                            >  1:465923/1‑251 (MQ=255)
             gTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGc                                                                                                                                                                                                                                     <  1:480795/251‑1 (MQ=255)
                     tACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTAGCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCa                                                                                                                                                                                                                              <  1:547962/250‑1 (MQ=255)
                                 cGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAatt                                                                                                                                                                                                                 >  1:197566/1‑251 (MQ=255)
                                    tACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGACTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAattatt                                                                                                                                                                                                              >  1:205880/1‑251 (MQ=255)
                                        ggAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTcgcg                                                                                                                                                                                                          >  1:432126/1‑251 (MQ=255)
                                               aaCGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCggg                                                                                                                                                                                                   >  1:161118/1‑250 (MQ=255)
                                                aCGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTg                                                                                                                                                                                                  <  2:197566/251‑1 (MQ=255)
                                                aCGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCCTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTg                                                                                                                                                                                                  <  2:170321/251‑1 (MQ=255)
                                                        gctcataGTATCATGGCCGTTTGTTACAAGGCTTTGATGGATTCTTAGCGGCAGTGCCTCGAGTCGTTCACTACTGTCCGACTCACGCGGCAGCAGGCGTATCTCCCCGCCGCTTTGAGGGAAATGCTGGATACGTTACAAAATCAAACCTGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCTTTGTACTTGATGCGCAAAGTCAGGTTCATTTGTTGCAACTTATTCGCTACGCGCTGCTGAATGc                                                                                                                                                                                         <  2:356140/245‑1 (MQ=255)
                                                                   cATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacg                                                                                                                                                                              <  1:377314/251‑1 (MQ=255)
                                                                        ccGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacgccaac                                                                                                                                                                         >  1:523619/1‑251 (MQ=255)
                                                                        ccGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacgccaac                                                                                                                                                                         <  1:26092/251‑1 (MQ=255)
                                                                          gATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacgccaacgc                                                                                                                                                                       >  1:554813/1‑251 (MQ=255)
                                                                            tttttCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGATCTTTTGTCTACTTTCCGCCTGACGCTGCAACAGGCGGATCTCCCCTCCGCATTGGGGGAAATGCTGGATACGTTACAAAAACAAACCAGCGCCAAACTGACCCTCCACTGCCGTCTGCCAACACTGGTACCGGAGGCGCAAAAGCAGGTGCATTTGGTTCAAATTCTTCGAGAAGCGTGGCTGAATGCGATGAAGCacgcccagaccc                                                                                                                                                                     >  2:38976/1‑244 (MQ=255)
                                                                                    agtcattGAATGATTCTTATCGGGAGTTTCGCGACCTGTTGACTTCTTTCCTCCTGTCGTTGCATCAGGCGGATCCCCCCGCCGCATTGAGGTAAATGTTGGATTCGTGACAAAAGCATGCCAGCGCCAAACTGACCCTCTACTGCCGTCTGCCAATCCTGGTACTTGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATCCGTGAAGCGGTGCTAGATGCGATGAAGCACGCCCACGTCAGCGAAATc                                                                                                                                                             <  2:87320/246‑1 (MQ=255)
                                                                                     ggCGTTGAATGATTCTTATCGGCAGTTACGCGAGCTGGTGACTACTTGCCGCCGGACGCTGCAGCAGGCGGATCACCCCTCCGCATTGAGGGAAATGCTGGACACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAAGTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCg                                                                                                                                                            <  2:205880/251‑1 (MQ=255)
                                                                                           gAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTc                                                                                                                                                       <  1:161804/250‑1 (MQ=255)
                                                                                            aaTGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTc                                                                                                                                                       >  1:162528/1‑249 (MQ=255)
                                                                                                           cagtgtcgcgAGCTGTTGACTACTGTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCTCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGcc                                                                                                                                      <  2:444330/245‑1 (MQ=255)
                                                                                                                aCGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACg                                                                                                                                 <  1:386329/251‑1 (MQ=255)
                                                                                                                         gTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATcac                                                                                                                         <  1:79214/250‑1 (MQ=255)
                                                                                                                           tGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCacac                                                                                                                       <  1:217423/250‑1 (MQ=255)
                                                                                                                           tGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTAGTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCacac                                                                                                                       <  1:58625/250‑1 (MQ=255)
                                                                                                                            gACTACCTGCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTCAGGGAAATGCTGGATACGTTACAAAATCAATCCAGCGCCAACCTGACCCTCGTCTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTGGCCAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACgg                                                                                                                     <  2:137162/251‑1 (MQ=255)
                                                                                                                             aagactTTCCGCCTGACTCGGCAGCAGGCGGAACCCCCCGCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAACCTGACCCTCGACTGCCGTCTGCCAACCCGGGGACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGt                                                                                                                    <  2:162528/248‑1 (MQ=255)
                                                                                                                                             cGCTGCAGCAGGCGGATCTCCCCACCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTTATAAc                                                                                                    <  2:345953/251‑1 (MQ=255)
                                                                                                                                                 gcagcaGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACggg                                                                                                 <  1:372241/250‑1 (MQ=255)
                                                                                                                                                  cagcagGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGCATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGAt                                                                                               >  1:342196/1‑251 (MQ=255)
                                                                                                                                                    gcagGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTg                                                                                             <  1:431084/251‑1 (MQ=255)
                                                                                                                                                    gcagGCGGATCTCCCCACCGCATTGAGGAAAATGCTGGATACGGGAAAAAAACAAACCAGCGCCAAAATGACCACAGACTGCCGTCTGCCAACCCTGGTACTGGGAGCGCCAATGAAGGGGCATTTGTTGCAAATTATTCGCGCAGCGTGTCTGAATACGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCATTTGCGTAACCGCGCCGGACGGCACACACCCGGTTTATATCCTTGTTAACGGGATTg                                                                                             >  2:274721/1‑251 (MQ=255)
                                                                                                                                                      agGCGGATCTCCCCTCCGCATTGAGGAAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGt                                                                                           >  1:150596/1‑251 (MQ=255)
                                                                                                                                                        gCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGCCCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGtct                                                                                         >  1:185580/1‑249 (MQ=255)
                                                                                                                                                                  ccTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCg                                                                               <  1:520189/251‑1 (MQ=255)
                                                                                                                                                                      cGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTTCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAg                                                                           <  2:366132/251‑1 (MQ=255)
                                                                                                                                                                       gCATTGAGGGAAAGTCTGTGTTCGTTACAAAATCAATCCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCTCAAATGCAGTTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCGCCGCGCCGGACGTCAAGCACACGGTGTATATCCGTGATAACGGGATTGTTATCGTTGAACTGAAAGa                                                                          <  2:520454/251‑1 (MQ=255)
                                                                                                                                                                            gAGGGAAATGCTGGATACGTTTCAAAATCAACCCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTCCTTGCTGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGGGAACCGAAAGATccc                                                                      <  2:554813/250‑1 (MQ=255)
                                                                                                                                                                                  aaTGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTc                                                               >  1:426342/1‑251 (MQ=255)
                                                                                                                                                                                      cTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGAATGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGGTAATATCCGTGATAACGGGATTGGTATCGGTGGACCGacaatacccaaaggacctg                                                            >  2:315594/1‑232 (MQ=255)
                                                                                                                                                                                        ggATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTAt                                                          <  1:257564/250‑1 (MQ=255)
                                                                                                                                                                                        ggATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATg                                                         <  2:237309/251‑1 (MQ=255)
                                                                                                                                                                                                   aaaaTCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGa                                                                                                                                                                 >  1:229605/1‑136 (MQ=255)
                                                                                                                                                                                                   aaaaGCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGa                                                                                                                                                                 <  2:229605/136‑1 (MQ=255)
                                                                                                                                                                                                      agcAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTTCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCACTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGTACCCGAAGGTCATTATGGTCTGAATATCATg                                           <  2:212819/249‑1 (MQ=255)
                                                                                                                                                                                                                     aaCTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGa                             >  1:266121/1‑250 (MQ=255)
                                                                                                                                                                                                                       cTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCCGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCAAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACCCGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACg                          >  2:277789/1‑251 (MQ=255)
                                                                                                                                                                                                                          aCCCGCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCAGTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTAGCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGc                        <  1:489800/250‑1 (MQ=255)
                                                                                                                                                                                                                           cacTCGACTGCAGTCGGCCAACACTGGTATGGGATGCGGAAATGACGGTGCATTTGTTGCACATTATTCGCGAAGCGGTTCAGAATGAGCTGCAGCACGACAACGCCAGCGAAATCGCCGTCCGTTGATACTCCGCGCCGGACGGCAATCACCCGGCTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCAGAAGGTCATTATGGTCTGAATTTCATGCGCGAACGCGCGGAACGGCTa                      <  2:494402/249‑1 (MQ=255)
                                                                                                                                                                                                                                gACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGCTTATATCCGTGATAACGGGATTGGTCTCGGTGACCCGAAAGAACCCGACGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTgg                 >  2:437957/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cgggtttGCCAATCCTGGACTTGGATGCGCAAAGTCATGTGCATTTGTTGCAAATGTTTCGCGTAGCGGTGCTGAATGCGATGGAGCTCGACAACTCCAGCGAAATCGCCGTCAGTTGCTTCTCCGCGCCGGACGGGAATCACACGGTTTATAACCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAGCCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCTCGGAACGGCTAGGTGGGaa              <  2:202227/244‑2 (MQ=255)
                                                                                                                                                                                                                                        tCTGCCAACCCTGGAACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGATAGGAGGGACGCTg          >  1:126544/1‑250 (MQ=255)
                                                                                                                                                                                                                                               aCCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGATTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTc  >  1:313570/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
ACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTC  >  NC_000913/2586866‑2587355

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: