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

TACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGAC  >  NC_000913/2586869‑2587349
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             
tacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGTGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGa                                                                                                                                                                                                                                         <  2:497083/250‑1 (MQ=255)
 acGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTACGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATg                                                                                                                                                                                                                                       <  1:100245/251‑1 (MQ=255)
    tATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCAGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATgcgc                                                                                                                                                                                                                                    >  1:348672/1‑251 (MQ=255)
        caggtgatgtgaCTGAATCGTTCGATACCTGAATATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCTGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAAt                                                                                                                                                                                                                                <  2:430788/240‑1 (MQ=255)
         aGTTGACGTTACTGAAGCGTTCGATCCCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGCCCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATg                                                                                                                                                                                                                               <  1:503497/251‑1 (MQ=255)
               cGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTg                                                                                                                                                                                                                         <  1:402500/251‑1 (MQ=255)
                   aCTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCACTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCAtt                                                                                                                                                                                                                     <  1:447371/251‑1 (MQ=255)
                         gCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGCCTGCCAACCCTGGCCCTGGATGCGCAAATGCAGGTGCATTTTTTGc                                                                                                                                                                                                               >  2:85151/1‑251 (MQ=255)
                             tCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGGTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGGACTGGATGCGCCAATGCAGGTGCATTTGTTGCAAat                                                                                                                                                                                                           >  1:487002/1‑251 (MQ=255)
                                   ccGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTATTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCGCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTcg                                                                                                                                                                                                     <  2:175874/251‑1 (MQ=255)
                                       aaGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCATGGTACTGGATGCGCAAATGCCGGCGCATTTGTTGCCAattatttgccg                                                                                                                                                                                                  >  1:292412/1‑245 (MQ=255)
                                            aaCGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGATCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGGGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCggg                                                                                                                                                                                            >  1:233862/1‑250 (MQ=255)
                                                ccACCGCACAAAGTATCATGGCCGATTTTTCCCCGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTTACTACTTTCCGCCTGACGCTGCAGCCGGCGGATCTCCCCTCCGCCTTGAGGGGAATGCTGGATACGTTACAAAATCAAACCAGCGCCCAACTGACCCTCGACTGCCGTCTGCCCACCCTGGTACTGGATGCGCCAAAGCAGGTGCCTTTGTTGCAAATTATTCTCGAAGCGGGGCTg                                                                                                                                                                                        >  2:28852/1‑251 (MQ=255)
                                                       acaAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCgg                                                                                                                                                                                 >  1:486886/1‑250 (MQ=255)
                                                          aaGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCCGGCGGATCTCCCCTCCGCATTGAGGGCAATGCTGGATACGTTACAAAATCAAACCCGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATg                                                                                                                                                                               >  1:266068/1‑250 (MQ=255)
                                                           aGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGa                                                                                                                                                                              >  1:410248/1‑250 (MQ=255)
                                                                 aTGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGGCTACTTTCCGCCTGACGCTGCAGCCGGCGGATATCCCCTCCGAATTGAGGGAAATGCTGGGTACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGGAATGGATGCGCAAATGCAGGGGCCTTTGTTGCAAATTATTCGCCGAGCGGTGCTGAAAGCGGTGAAGccagc                                                                                                                                                                       >  2:93127/1‑247 (MQ=255)
                                                                       gATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGGACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGCTGAAGCacgccaacgc                                                                                                                                                                 >  1:209963/1‑251 (MQ=255)
                                                                            ttCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCTAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGc                                                                                                                                                             <  2:263203/250‑1 (MQ=255)
                                                                               ccAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAACTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGATGGATACGGTACAAAATCAAACCAGCGCCAAACTGACCCTCAACTGCCGACTGCCCACCCTAGTCCTGGCTGCGCAAATGCAGGTCAATTTGTTTGAAATTATTCGGGAAGCGGGGCTGGATTCGGTGAAGAaagccaaagaacacgaa                                                                                                                                                          >  2:301990/1‑240 (MQ=255)
                                                                                         aaTGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCCCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAAACAGCGCCAAACTGCCCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCCAATGACGGTACATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATACGATGAAGCacggcaacgccaccggaaacgcccgc                                                                                                                                                 >  2:154719/1‑235 (MQ=255)
                                                                                                   aTCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCAcc                                                                                                                                     <  2:31206/251‑1 (MQ=255)
                                                                                                     cGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGAGCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACcg                                                                                                                                    <  2:457537/250‑1 (MQ=255)
                                                                                                      ggcggTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCCACAGGCGGATCTCCCCTCCGCATTGAGGGAACTGCTGGATACGTTACAAAATCAACCCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCACCCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGCAGCAGTGATGAATGCGATGAAGCCCACCAACGCCAGCGAAATCGCCGGCAGTTGCGTCAccccg                                                                                                                                  >  2:393636/5‑248 (MQ=255)
                                                                                                        aaGTTACGCGGGCTGTTGACTACTTCCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGcc                                                                                                                                <  2:280150/250‑1 (MQ=255)
                                                                                                             aCGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACg                                                                                                                           >  1:87136/1‑251 (MQ=255)
                                                                                                             aCGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGAGGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGAc                                                                                                                            <  1:93127/250‑1 (MQ=255)
                                                                                                                         gtcttctTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCCTTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGCCGGCAATCACACg                                                                                                                <  2:299829/245‑1 (MQ=255)
                                                                                                                             actTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGttt                                                                                                            >  1:214402/1‑250 (MQ=255)
                                                                                                                               tttCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTtat                                                                                                          >  2:64596/1‑250 (MQ=255)
                                                                                                                                             tGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACggg                                                                                           >  1:354366/1‑251 (MQ=255)
                                                                                                                                                agcagGCTGACCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGAtt                                                                                        <  1:338915/251‑1 (MQ=255)
                                                                                                                                                agcagGCGAATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGAACTGGATGCGCAAATGCAGGAGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATACGTGATAACGGGAtt                                                                                        >  2:306476/1‑251 (MQ=255)
                                                                                                                                                     gCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTAt                                                                                   <  1:529759/251‑1 (MQ=255)
                                                                                                                                                             ccccTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAAc                                                                           <  2:87136/251‑1 (MQ=255)
                                                                                                                                                               ccTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCCCACGGTTTATATCCGTGATAACGGGATTGGAATCGGTGAACCg                                                                         >  2:350993/1‑251 (MQ=255)
                                                                                                                                                                                   cTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCGACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCTGTGAACCGAAAGAACCCGAAGGTCAtt                                                      >  2:162620/1‑250 (MQ=255)
                                                                                                                                                                                             gaCAAAAACAATCCAGGGCTAAACTGACCCTCTGCTGCCGTCTGCCAACCCTGTTTCGGGTTGCGCAAATGCAGGTGCATTTGTTGCAAAGTATTCGCTAATCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGTGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGaa                                            <  2:396577/249‑1 (MQ=255)
                                                                                                                                                                                              aCAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGACTGCGATGAAGCACGCCAGCGCGAGCGACATCGCCGTCAGCTGCGTCCCCGAGCCGGAAGGCAATCACACGGTTTATCTCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGACGGACCTTATGGTCTGAatt                                          >  2:205921/1‑250 (MQ=255)
                                                                                                                                                                                                    tCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCTCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATgc                                    >  1:434104/1‑251 (MQ=255)
                                                                                                                                                                                                      aaaCCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATgcgc                                  >  1:523472/1‑251 (MQ=255)
                                                                                                                                                                                                                cAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATAGTCTGAATATCATGCGCGAACGCGCgg                        >  2:368251/1‑251 (MQ=255)
                                                                                                                                                                                                                  caCTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGCAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGACGGACGGCAATCACACGGTTTATCTCCGTGATACCGGGATTGGTATCGGTGAACCGAAAGAACCCTAAGGTCATTATGGTCTGAATATCATGCGCGAACTCGCGGa                       >  2:95725/2‑250 (MQ=255)
                                                                                                                                                                                                                  aaCTGACCCTCCACTGTCGTCTGCGAACCCAGGTACTGGATGCGCAAATTGATGTGCAATTGTTGCAAATTATTCGCGAAGCGGGGCTGAATGCGATGAAGCGCGCCAACGCCAGCGAAATCGCCGGCAGTTGCGTCACCGCGCCGGACGGCAATCACCCGGTTTATATCCGTTATAACGGGATTGGTATCGGTGAACCGAAAGAACACGAAGGTCAATAGGGTCTGAAGATCATGCGCGACCGCGCGGaa                      <  1:324439/251‑1 (MQ=255)
                                                                                                                                                                                                                      gccccTCGGCTGCCGTCTGGCAACCCTGGTACTGGAGGCGCAAATGCAGGTGCATTTGTAGCAAATTATTCGCGAAGCGGGGCTGAAGGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTGGCGGCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGTTCTGAATATCATGCGCGAACGCGCGGAACgg                   <  2:487002/248‑1 (MQ=255)
                                                                                                                                                                                                                      gACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACgg                   <  1:202216/250‑1 (MQ=255)
                                                                                                                                                                                                                      gACCCCCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGc                  <  2:171319/251‑1 (MQ=255)
                                                                                                                                                                                                                            cgccTGCCGTCTGCCAACCCTGGTACTGGAGGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCTCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGt             <  2:434104/247‑1 (MQ=255)
                                                                                                                                                                                                                              ccTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGCGCATTTGTTGCAAATTATTAGCGCAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGGTTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGGTAGGTggg          >  2:57612/2‑251 (MQ=255)
                                                                                                                                                                                                                                ggCCGTCTGCCATCCCAGGTACTGGATGCGCAAATGCAGGTTCATTTGTTGCAAAGTATTCGCGAAGCGGTGCTGAATGCGAGGAAGCACGCCAACGCCAGCGAAATCGGCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGAc        <  2:249001/250‑1 (MQ=255)
                                                                                                                                                                                                                                      cagCCAACCCCGGTACCGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGAc  <  2:486886/249‑1 (MQ=255)
                                                                                                                                                                                                                                       tGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATAAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTg            <  2:139801/240‑1 (MQ=255)
                                                                                                                                                                                                                                       tGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATAAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTg            >  1:139801/1‑240 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             
TACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGAC  >  NC_000913/2586869‑2587349

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