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,1120CT97.0% 204.8 / NA 68A461V (GCA→GTA) narQsensor histidine kinase NarQ
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (30/34);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCA  >  NC_000913/2586870‑2587356
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
acGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGAt                                                                                                                                                                                                                                                <  2:256783/250‑1 (MQ=255)
    aTCCAGTTGACGTTTCTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGAGCTCCCCGCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATgcgc                                                                                                                                                                                                                                            <  2:572461/250‑1 (MQ=255)
    aTCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCa                                                                                                                                                                                                                                           <  1:153971/251‑1 (MQ=255)
       catttGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAAt                                                                                                                                                                                                                                        <  1:770767/248‑1 (MQ=255)
            gACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAgg                                                                                                                                                                                                                                  <  1:74558/251‑1 (MQ=255)
            gACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAgg                                                                                                                                                                                                                                  >  1:489632/1‑251 (MQ=255)
             aCGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGt                                                                                                                                                                                                                                 <  1:477915/251‑1 (MQ=255)
             aCGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTAACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGt                                                                                                                                                                                                                                 <  1:766513/251‑1 (MQ=255)
                   cggAAGCGTTCGATACCGGAAGATAACGCCAACGCACAAAGAATCAAGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTTACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCGCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACCGGATGCGCAAATG‑CAGGTGCAttt                                                                                                                                                                                                                           <  2:446357/249‑1 (MQ=255)
                   cTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGGGCAttt                                                                                                                                                                                                                           >  1:415392/1‑251 (MQ=255)
                       aGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGGGCTGTTGACTACTTTCCGCCTGACGCTGCAGCCGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGCTACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGGCTGCCGTCCGCCCACCCTGGTGCTGGATGAGCAAATG‑CACGTGCATatgttg                                                                                                                                                                                                                       >  2:677952/1‑251 (MQ=255)
                          gTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGCTTTTTACCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGGTGTTGCCTACTTTCCGCCTGACGCTGCAGCAGGCGGGACTCACCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGCCCCTAGGCTGCCGTCTGCCAACCCTGGCACTGGCTGCGCCAATG‑CAGGTTCATTTGTTGCa                                                                                                                                                                                                                     >  1:463820/1‑250 (MQ=255)
                            tCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCaaa                                                                                                                                                                                                                   <  2:112797/250‑1 (MQ=255)
                            tCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGATAGCTGTTGACCACTTTCCGCCTGACGCTGCAGCAGGCGGGTCTCCCCCCCGCATTGAGGGAAATGCTGGAAACGTTACAAAATCAAACCCGCGCCAAACTGACCCTCGACTGCCGTCTGCCAAACCTGGGCCTGGATGGGCAA‑TGACCGGGTCATTTGTTGCAAat                                                                                                                                                                                                                  >  2:579419/1‑251 (MQ=255)
                               aTACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTGCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAattat                                                                                                                                                                                                               <  1:461902/251‑1 (MQ=255)
                                 aCCGGAAGATAACGCCTCCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTc                                                                                                                                                                                                             <  1:751778/251‑1 (MQ=255)
                                          tAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTTCTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTAACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGCGCCTTTGTTGCAAATTATTTGCGAAGCgg                                                                                                                                                                                                    >  1:37474/1‑251 (MQ=255)
                                                       aaaaGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGAt                                                                                                                                                                                       <  2:443715/250‑1 (MQ=255)
                                                        aaaGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATg                                                                                                                                                                                      >  1:547553/1‑251 (MQ=255)
                                                         aaGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCCGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCCCGCCTGCCGTCTACCAAACCTGGTACTGGCTGCGCAAATG‑CAGGGGCTTTTTTTGCAAATTATTACCGAAGTGGGGCTGatagcgtggg                                                                                                                                                                                     >  2:532785/1‑242 (MQ=255)
                                                               cATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCac                                                                                                                                                                                <  1:58234/250‑1 (MQ=255)
                                                                              ccAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGGTGCGCAAATG‑CAGGTGCATTTGTTGCCAATTATTCGCGCAGCGGTGCTTAATGCGATGAAGCACGCCCACGCCAGCGaa                                                                                                                                                                 >  1:622430/1‑250 (MQ=255)
                                                                                    gTTGAATGATGCTTATCGGAAGTTAAGCGTGATGTTGACTCCTTTCCGCCTGACGCTGCAGCAGGCGGGTATCCCCCCCGCATTTAGGGAAATGCTGGGTACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCCCACATG‑CAGGTGCATTTGTTGCAAATTATTCCCGAAGCGGTGATGGATGCGACGAAGCACGGCAAAGCCAGCGAAAACGccc                                                                                                                                                          >  2:491516/1‑250 (MQ=255)
                                                                                      tGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTc                                                                                                                                                        >  1:347596/1‑251 (MQ=255)
                                                                                            atgCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGCCTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCAAAGCGATGCTGCATTCGATGAAGCACGCCAACGCCAGCGAAATCGACGTCAGTTGc                                                                                                                                                  >  1:508590/1‑251 (MQ=255)
                                                                                              gCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGt                                                                                                                                                >  1:628849/1‑251 (MQ=255)
                                                                                                   tcgacAGTTACGCGAGCTGTTGACTACATTCCGCCTAACGCTGCAGCAGGCGGATCTCCCCTACGCATTGAGGGAACAGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGGCTGCCAACCCTGGAACTGGATGCGCAAATG‑CAGGAGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCAAGCCAACGCCAGCGAAATCGCCGGCGGTTGCGTCACcg                                                                                                                                           >  2:458403/5‑251 (MQ=255)
                                                                                                       cAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCCCCGCGcc                                                                                                                                       >  2:227078/1‑251 (MQ=255)
                                                                                                        aGTTACGCGAGCTGTTGACTCCTTTCCGCCTGACGATGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACACTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGt                                                                                                                                                                                                                                 >  2:361063/1‑160 (MQ=255)
                                                                                                        aGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAAGCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGt                                                                                                                                                                                                                                 <  1:361063/160‑1 (MQ=255)
                                                                                                            ccgcgAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCGCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACg                                                                                                                                  <  2:646733/250‑1 (MQ=255)
                                                                                                            aCGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGAc                                                                                                                                   <  1:522445/250‑1 (MQ=255)
                                                                                                            aCGCGAGCTGTTGACTACTTGCCGCCTGACGCGGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGGTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTAGTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGAc                                                                                                                                   <  2:15307/250‑1 (MQ=255)
                                                                                                               cgAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCa                                                                                                                               <  1:751755/251‑1 (MQ=255)
                                                                                                               cgAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCa                                                                                                                               <  1:699555/251‑1 (MQ=255)
                                                                                                                        gACTACTTTCCGTCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACgg                                                                                                                      <  1:53598/251‑1 (MQ=255)
                                                                                                                        gACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACgg                                                                                                                      <  2:614475/251‑1 (MQ=255)
                                                                                                                         ccttccttCCGCCTGACGTCGCCGCAGACGGATCCCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCTAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGt                                                                                                                     <  2:556892/245‑1 (MQ=255)
                                                                                                                            actTTCCGCCTGACGCTGCAGCAGGCGGATCCCCCCTCCGCATCGAGGGAAATGCGGGATACGTTACAAAATCAAACCAGCGCCAAAATGACACTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAACTATTCGCGAAGTGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGGCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGttt                                                                                                                   <  2:225001/250‑1 (MQ=255)
                                                                                                                               ttCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCACCGCCAGCGAAATCGACGTCAGTTGCGTCCCCGCCCCGGGCGGAAATCACACGGTTTttat                                                                                                               >  2:113474/1‑251 (MQ=255)
                                                                                                                                  cGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTcgc                                                                                                                                                                                                           >  1:783500/1‑156 (MQ=255)
                                                                                                                                  cGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTcgc                                                                                                                                                                                                           <  2:783500/156‑1 (MQ=255)
                                                                                                                                    ccTGACGCTGCATCAGTCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAAACAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTg                                                                                                          <  2:27880/251‑1 (MQ=255)
                                                                                                                                       gACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCACCGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGAt                                                                                                        >  2:766736/1‑250 (MQ=255)
                                                                                                                                       gACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGAt                                                                                                        >  2:281792/1‑250 (MQ=255)
                                                                                                                                             gcagcaGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCCCGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGCAGGCAATCACACGGTTTATATCCGTGATAACgggg                                                                                                 >  2:163623/1‑250 (MQ=255)
                                                                                                                                                    gCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTAt                                                                                          >  1:262652/1‑251 (MQ=255)
                                                                                                                                                    gCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTAt                                                                                          >  1:655577/1‑251 (MQ=255)
                                                                                                                                                           tcCCCTCCGCATTGCGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTACCAACACTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCACATTATTCGCGACGCAGTGCTGCATGCGATGAAGCACGCCAACGCCAACCAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTTGTAACgggcacggtctcaggga                                                                                    >  2:683479/1‑236 (MQ=255)
                                                                                                                                                           cccccACCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGGGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGaa                                                                                   <  2:415392/250‑1 (MQ=255)
                                                                                                                                                             cccTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGGCACCGCGCCGGACGGCAATCACACGGTTTAAATCCGTGATAACGGGATTGGTATCGGTGAAcc                                                                                 >  2:1407/1‑251 (MQ=255)
                                                                                                                                                                      ttGAGGGAAATGCTGGATACGTTACAAACTCAAACCAGCGCCAACCTGACCCTCGACTGCCGTCTGCAAACCCTGGTACTGGATGCGCAAATG‑CAGGAGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCCCCGCGCAAGACGGCAATCACACGGTTTATATCCGTGATAACGTGCTTGGTATCGGTGAACCGAAAGAAcc                                                                        >  2:549117/1‑251 (MQ=255)
                                                                                                                                                                          gggAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGa                                                                     >  1:338708/1‑250 (MQ=255)
                                                                                                                                                                                    ggATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATg                                                          <  1:491516/251‑1 (MQ=255)
                                                                                                                                                                                      aTACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGACATTATGGt                                                        >  1:733927/1‑251 (MQ=255)
                                                                                                                                                                                         cGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAAGGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTAGATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTg                                                     <  1:56478/251‑1 (MQ=255)
                                                                                                                                                                                          gTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGa                                                    <  1:721176/251‑1 (MQ=255)
                                                                                                                                                                                            tACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAat                                                  >  1:379694/1‑251 (MQ=255)
                                                                                                                                                                                                  aTCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATg                                            >  1:128738/1‑251 (MQ=255)
                                                                                                                                                                                                    cAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATgcg                                          >  1:131276/1‑251 (MQ=255)
                                                                                                                                                                                                      aaCCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACGGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGcgcg                                        <  2:495521/251‑1 (MQ=255)
                                                                                                                                                                                                         cAGCGCCAAACTGACCCTCGACTGCCGTCTACCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAAc                                     <  1:532785/251‑1 (MQ=255)
                                                                                                                                                                                                             gcCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAAcgcgc                                 >  1:137104/1‑251 (MQ=255)
                                                                                                                                                                                                                        ccTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAg                      <  2:432927/251‑1 (MQ=255)
                                                                                                                                                                                                                           cgccTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGGCAACGCCAGCTAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTg                   <  1:602532/248‑1 (MQ=255)
                                                                                                                                                                                                                                       gCCAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGAGAACGCGCGGAACGGCTCGGTGGGACGCTGACtt       >  2:493990/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cAACCCTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGGCGCTGACtttt     >  2:92341/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cccTGGTACTGGATGCGCAAATG‑CAGGTGCATTTGTGGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTca  <  2:780230/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
ACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATG‑CAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCA  >  NC_000913/2586870‑2587356

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: ‑