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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,587,1120CT89.6% 113.0 / 5.1 52A461V (GCA→GTA) narQsensor histidine kinase NarQ
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base T (15/28);  total (20/32)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

TACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATG‑CGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCAC  >  NC_000913/2586865‑2587357
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
tacttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCTCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACt                                                                                                                                                                                                                                                     >  1:245554/1‑251 (MQ=255)
              gTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGGTACGCGGGCTGTTGAATACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGGGGGAAATGGTGGATAAGTT‑CAACAAACCAACACAG‑GCCAAACAGGCCCCCGGCTGCCGTCTGCCAAACCTGGGACTGGAATCGCAAATgg                                                                                                                                                                                                                                       >  2:488877/1‑250 (MQ=255)
                  aCGTTACTGAAGCGGTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGAGTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTGCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTTGTACTGGATGCGCAAATGCAGGt                                                                                                                                                                                                                                   <  2:133344/251‑1 (MQ=255)
                         ggAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCAGTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttg                                                                                                                                                                                                                            <  2:234329/250‑1 (MQ=255)
                          gAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCAttt                                                                                                                                                                                                                             <  2:379936/249‑1 (MQ=255)
                                  cGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAatt                                                                                                                                                                                                                   <  2:210176/251‑1 (MQ=255)
                                      aCCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTc                                                                                                                                                                                                               <  1:488877/251‑1 (MQ=255)
                                        cGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTGCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTcgc                                                                                                                                                                                                             <  2:82018/251‑1 (MQ=255)
                                         ggAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTcgcg                                                                                                                                                                                                            <  1:306296/251‑1 (MQ=255)
                                          gAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCCTTGGGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGCTGCCCAAACGCAGGTGCCTTTGTTGCAAATTATTcgcgg                                                                                                                                                                                                           >  2:132723/1‑250 (MQ=255)
                                                                gTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAg                                                                                                                                                                                     <  1:382543/251‑1 (MQ=255)
                                                                  aTCATGGCCGATTTTTCCCAGGCGTTGAATGATGCGTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGc                                                                                                                                                                                    <  2:374585/250‑1 (MQ=255)
                                                                    cATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCac                                                                                                                                                                                  <  1:453543/250‑1 (MQ=255)
                                                                               tttCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGATCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAg                                                                                                                                                                       <  1:479948/250‑1 (MQ=255)
                                                                                   ccAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTTTTTACCACTTTCCGCCTGGATCTGCAGCAGGCGGATCTCCCCTCCGAATTTGGGGAAATGCTGGATACGGTACAA‑AATCACAC‑CAGCGCCAAACTGCCCCTCGCATGCCGTCTTCCCACCCTGGAAATGGATGCACAAATGCAGGCGCATTTGTTGCACATTATTCGAGAAGCGGTGCTGAATGCGATGAAGCACGCCCAAGCCAGccaa                                                                                                                                                                   >  2:164371/1‑247 (MQ=255)
                                                                                        cGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGcc                                                                                                                                                             >  1:415264/1‑251 (MQ=255)
                                                                                        cGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGGGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGcc                                                                                                                                                             >  1:27327/1‑251 (MQ=255)
                                                                                                           gCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGGACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAATTCGCCGCCAGTTGCGTCACcgcac                                                                                                                                          >  2:212395/1‑249 (MQ=255)
                                                                                                                          gTTGACTACTTCCCGCCTGACGCTGCAGCAGGCGGATCGCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATcac                                                                                                                            <  2:72986/250‑1 (MQ=255)
                                                                                                                              actactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGt                                                                                                                       >  1:476938/1‑251 (MQ=255)
                                                                                                                                    ttccgccgggagCTGCAGCAGGCGGAATTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCAGTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTtat                                                                                                                   <  2:367918/238‑1 (MQ=255)
                                                                                                                                       cgacggACGCTGCAGAAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCTAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCg                                                                                                              <  2:125871/246‑1 (MQ=255)
                                                                                                                                           tGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacgccaac                                                                                                                                                                            >  1:257008/1‑185 (MQ=255)
                                                                                                                                           tGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacgccaac                                                                                                                                                                            <  2:257008/185‑1 (MQ=255)
                                                                                                                                              cGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCAAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAAGTATTCGCGAAGCGGTGCTGAATGCGAAGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTTATAAc                                                                                                       <  1:347448/251‑1 (MQ=255)
                                                                                                                                                  gcagcaGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACgg                                                                                                     <  1:204712/249‑1 (MQ=255)
                                                                                                                                                  gcagcaGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACgg                                                                                                     <  2:177686/249‑1 (MQ=255)
                                                                                                                                                          cGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATc                                                                                           >  1:175346/1‑251 (MQ=255)
                                                                                                                                                          cGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATc                                                                                           <  2:175346/251‑1 (MQ=255)
                                                                                                                                                                       cGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGaacg                                                                              >  1:262203/1‑249 (MQ=255)
                                                                                                                                                                         cATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGa                                                                             >  1:314723/1‑250 (MQ=255)
                                                                                                                                                                          aTTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGCTTGGTATCGGTGAACCGAAAGAAc                                                                           >  2:299308/1‑251 (MQ=255)
                                                                                                                                                                            ttaGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAAcc                                                                          <  2:29668/248‑1 (MQ=255)
                                                                                                                                                                            tGAGGGAAATGCTGGATTCTTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCCCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGGATCGGTGAACCGAAAGAAcc                                                                          <  2:86750/250‑1 (MQ=255)
                                                                                                                                                                             gAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCg                                                                        <  1:107278/251‑1 (MQ=255)
                                                                                                                                                                                 gAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGTTTGGTATCGGTGACCCGAAAGCACCCGAAgg                                                                    >  1:229407/1‑251 (MQ=255)
                                                                                                                                                                                    aTGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCATCCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGTTATCGGTGAACCGAAAGAACCCGAAGGTCa                                                                 <  1:449224/251‑1 (MQ=255)
                                                                                                                                                                                        tGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGACTGCCAACCCTGGAACTGGATGCGCCAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGACAACGACAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGACAGAACCCGAAGGACATTAt                                                             >  1:499663/1‑251 (MQ=255)
                                                                                                                                                                                         ggATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGCAAGAACACGAAGGACATTATg                                                            >  1:455737/1‑251 (MQ=255)
                                                                                                                                                                                            cacgctACAA‑AATCAAAC‑CAGCGCCAAAATGTCCCTCGACTGTCGTCTGCCATCCCTGGTACAGGATGCGCAAATGCAGGTGCATTTGGTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCTGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGATCCCGAAGGTCATTATGGTc                                                         <  2:187088/246‑1 (MQ=255)
                                                                                                                                                                                               gTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGaa                                                                                     >  1:160132/1‑220 (MQ=255)
                                                                                                                                                                                               gTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGaa                                                                                     <  2:160132/220‑1 (MQ=255)
                                                                                                                                                                                                 tACAA‑ACTCAAAC‑CAGCGCCAAACTGACCCGCGACTGCCGTCTGCAAACCCTTGTACTGGATGCGCAAATGCAGGTGCAGTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGACAACGCCAGCGAAATCGCCGGCAGTTGCGTCACCGCGCCGGACGGCTATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGaa                                                     <  2:146335/250‑1 (MQ=255)
                                                                                                                                                                                                          cAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGTTTCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTGGCGTCACCGCGCCGGACGGCAATCACGCGGTTTATATCCGTGATAACTGGATTGGTATCGGTGTACCCAAAGAACCCGAAGGTCATGATGGTCTGAATATCATg‑cg                                           <  1:55947/251‑1 (MQ=255)
                                                                                                                                                                                                          cAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGGCTGAATATCATg‑c                                            <  1:394103/250‑1 (MQ=255)
                                                                                                                                                                                                              c‑cAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGAtt                                                                                                 >  2:199019/1‑194 (MQ=255)
                                                                                                                                                                                                              c‑cAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGAtt                                                                                                 <  1:199019/194‑1 (MQ=255)
                                                                                                                                                                                                                   cgcCAAACTGATCCTCGACTGCCTTCTGGCAACCCTGGTACTGGATGTGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCGCCGCGCCGGACTGCAATCACACGGGTTATATCCGTGATAACGGGATTGGTATCGGTGATCCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATG‑CGCGAAcgcg                                   <  2:241487/251‑1 (MQ=255)
                                                                                                                                                                                                                     ccAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGAGCATTTGTTGCAAATTATTCGCGAAGCAGTGCTGAATGCGCTGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCACTCACACGGTTTATATCCGTGCTACCGGGATTGGTATCGGTGAACCGAAGGCACCCTAAGGACATTCTGGGCTGAATATC‑TGTCGCGAAcgcgcg                                 >  2:138373/1‑251 (MQ=255)
                                                                                                                                                                                                                             acacTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATG‑CGCGAACGCGCGGAACGGCt                         <  1:337119/248‑1 (MQ=255)
                                                                                                                                                                                                                                                 aaCCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATG‑CGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACttttt     <  1:382773/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    ccTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATG‑CGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTcac  >  1:354634/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
TACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAA‑AATCAAAC‑CAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATG‑CGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCAC  >  NC_000913/2586865‑2587357

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