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,1120CT98.4% 174.5 / NA 62A461V (GCA→GTA) narQsensor histidine kinase NarQ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (26/34);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

TTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTG  >  NC_000913/2586868‑2587347
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           
ttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTgg                                                                                                                                                                                                                                        >  1:463625/1‑250 (MQ=255)
 tacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGAt                                                                                                                                                                                                                                      >  1:176877/1‑251 (MQ=255)
  acGTATCCAGTTGACGTTACTGAAGCGTTAGATACCGGAAGATAACGACACCGCACAAAGTATCATGGACGCTTTTTCAAAGGCGTTGAATGATGCTTATCGGCAGTTACGAGAGCTGTTGACTACTTTCCGCCTGAAGATGCCGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAAACAGCGCCAAACTGACCATCGACTGCCGTCTGCCAACCCTGGTACTggctg                                                                                                                                                                                                                                     >  1:356321/1‑248 (MQ=255)
         cAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGGCGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACTGTACAAAGTCAAACCAGCGCCAACCTGACCCGCAACTGCCGTCTTCCAACCCTGGTACTGTATGCGCAAAt                                                                                                                                                                                                                              <  2:383855/251‑1 (MQ=255)
          aGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATg                                                                                                                                                                                                                             <  2:463625/251‑1 (MQ=255)
             ggACGTTACTGAAGCGTTCTATTCCGGAAGATAACGCCTCCGCACAAAGTATCATGGCCGAGTTTTCCCAGGCGTTGAATGATGCTTAGCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCGGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACGGACCCTCGACGGCCGTCTGCCAACCCGGGTTCTGGATTCGCAAATGCAg                                                                                                                                                                                                                          <  2:388167/250‑1 (MQ=255)
              gACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCCTGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGCAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGAACCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAgg                                                                                                                                                                                                                         >  1:305651/1‑251 (MQ=255)
              gACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGAGGCGCAAATGCAg                                                                                                                                                                                                                          <  2:117700/250‑1 (MQ=255)
                cGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCTCAAATGCAGGTg                                                                                                                                                                                                                       <  2:427862/251‑1 (MQ=255)
                 gTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGc                                                                                                                                                                                                                      <  2:253199/251‑1 (MQ=255)
                    aCTGAAGCGTGCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACGTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAAGCAAACCAGCGCCAAGCTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCAtt                                                                                                                                                                                                                   <  2:4999/251‑1 (MQ=255)
                       taaGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTACCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCAttt                                                                                                                                                                                                                  <  2:473792/248‑1 (MQ=255)
                            gTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCa                                                                                                                                                                                                            >  1:508341/1‑250 (MQ=255)
                                gATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGGGCATTTGTTGCAAattt                                                                                                                                                                                                       >  2:207223/1‑250 (MQ=255)
                                 aTACCGGAAGATAACGCCAACGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGCGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGGTGCGCAAACGCAGGGGCCTTTGGTGCAAattat                                                                                                                                                                                                      >  2:499171/1‑251 (MQ=255)
                                        aaGATAACGCCATCGCACAAATGATCATGGCCGATTTATCCCAGGCGTTGAATGATGCTTGTCGGCAGTTACGCGAGCTGTTGACGACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCACATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTTACCCTCGACTGCCGTCTGCCACCCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATGATTCGCGaa                                                                                                                                                                                               <  1:113804/251‑1 (MQ=255)
                                           aTAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCg                                                                                                                                                                                            <  1:425848/251‑1 (MQ=255)
                                                  cACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGa                                                                                                                                                                                     <  2:387899/251‑1 (MQ=255)
                                                      gtacaAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCTGTTACGCGAGCTGTTGACTTCTTTCCGCCTGACGCTGCAGCAGGCGTATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCTCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGc                                                                                                                                                                                 <  1:22516/249‑1 (MQ=255)
                                                          aaaGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATg                                                                                                                                                                             <  1:290212/251‑1 (MQ=255)
                                                            aGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGCTACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGAGGTGCTGAATGCGATGa                                                                                                                                                                            >  1:198100/1‑250 (MQ=255)
                                                                    gtccGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacgcc                                                                                                                                                                    <  1:231779/248‑1 (MQ=255)
                                                                           ttttCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacgccaacgcc                                                                                                                                                              >  1:469497/1‑249 (MQ=255)
                                                                             ttCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGc                                                                                                                                                           <  1:211969/250‑1 (MQ=255)
                                                                               cccAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCTGATCTCCCCTCCGCATTGAGGGAAATGCTGGATTCGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGaa                                                                                                                                                        <  1:377873/251‑1 (MQ=255)
                                                                               cccAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGCTGAAGCACGCCAACGCCAGCGa                                                                                                                                                         >  2:337742/1‑250 (MQ=255)
                                                                                   ggCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCg                                                                                                                                                    <  1:583019/251‑1 (MQ=255)
                                                                                      gTTGACTGATGCTTATCGGCGGTAACGCGCGCTGTTGACCCATTTCCGCCTGACGATGCAGCACGCGGATCTCCCCTCCCCATTGAGGGAAATGCTGGATCCGTTACAACATCAACCCAGCGCCAAACTGAACCTCGCCTGCCGCCTGCCAACCCTGGTACTGGCTGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGACTGCGATGCAGCACGCCAACGCCAGCGGAATCGCCg                                                                                                                                                 >  2:27700/1‑251 (MQ=255)
                                                                                             gatgCTTATCGGCAGTTACGCGAGCTGTTGCCTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGCTGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTg                                                                                                                                          >  1:60353/1‑251 (MQ=255)
                                                                                                   tATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCCGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGCTGCAGCACGCCAACGCCAGCGCAATCGCCGTCAGGTGCGTCAc                                                                                                                                    >  1:3189/1‑251 (MQ=255)
                                                                                                      cGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACcg                                                                                                                                  <  2:241282/250‑1 (MQ=255)
                                                                                                                  gAGCTGTTGACTAATTTACACCTGAGGCTGCAGCCGGCGGATCTACCCTACGCATTGCGGGAAATGATGGCTAAGCTACAAACTCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACACTGGTACAGGATGAGCCAAAGCAGGGGCATTTGTTGAAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGAACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGAAGGCa                                                                                                                      >  2:518328/1‑250 (MQ=255)
                                                                                                                       gTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATcac                                                                                                                 >  2:88716/1‑250 (MQ=255)
                                                                                                                         tGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCAGTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCacac                                                                                                               <  2:459212/250‑1 (MQ=255)
                                                                                                                          gACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAAGCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATcacaagg                                                                                                             <  2:126427/251‑4 (MQ=255)
                                                                                                                                 ttCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACCGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCCCACGGTTtat                                                                                                        <  1:74528/249‑1 (MQ=255)
                                                                                                                                       cTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTg                                                                                                 >  1:524558/1‑250 (MQ=255)
                                                                                                                                       cTGACGCAGCCGAAGGCGGATCGCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATAAAAGCAGCGCCAAACTGAGACTCGACTGCCGTCTGCCAACCCTGGTACTGGATGAGCACATGCAGGCACATTTGTTGCCAATTATTCGCGAAGCGGTGCTGAATGCGGTGAAGCACGCCAACGCCAGCGAAATCGCCGTAAGTTGCGTCCCCGCGCCGGACGTCAATCACCCGGCTTAAATCCGTgg                                                                                                >  2:376700/1‑250 (MQ=255)
                                                                                                                                        tGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGAt                                                                                               >  1:97616/1‑251 (MQ=255)
                                                                                                                                           cGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAAc                                                                                            <  1:339102/251‑1 (MQ=255)
                                                                                                                                            gCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACg                                                                                           >  1:189700/1‑251 (MQ=255)
                                                                                                                                               gcagcaGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGa                                                                                        <  1:206882/251‑1 (MQ=255)
                                                                                                                                                cagcagGCGGATCTCCCCTCCGCAGTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATGCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGAt                                                                                       <  2:152558/251‑1 (MQ=255)
                                                                                                                                                                  cccGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGaa                                                                     <  1:499171/250‑1 (MQ=255)
                                                                                                                                                                       aTTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCACACTGAGCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGGGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAAACgaacgaac                                                                >  1:497582/1‑251 (MQ=255)
                                                                                                                                                                       aTTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTta                                                                                                         <  2:503394/210‑1 (MQ=255)
                                                                                                                                                                       aTTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTta                                                                                                         >  1:503394/1‑210 (MQ=255)
                                                                                                                                                                                    cTGGCTACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTATGCCAACCCTGGTACTGGATGCGCAAATGCAGGCGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCAGGGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCAtt                                                    >  1:347997/1‑250 (MQ=255)
                                                                                                                                                                                       gATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTTGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATgg                                                <  2:524558/251‑1 (MQ=255)
                                                                                                                                                                                            gTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAAGCGAAAGAACCCGAAGGTCATTATGGTCTGa                                           <  1:180417/251‑1 (MQ=255)
                                                                                                                                                                                             ttACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCTCGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTtata                                                                                                       <  1:287521/190‑1 (MQ=255)
                                                                                                                                                                                             ttACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTtata                                                                                                       >  2:287521/1‑190 (MQ=255)
                                                                                                                                                                                                       aaaCCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATgcgc                                >  1:330483/1‑251 (MQ=255)
                                                                                                                                                                                                           cAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCATTTGCGCCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGAGAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAAc                            <  1:518225/251‑1 (MQ=255)
                                                                                                                                                                                                            aGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAAcg                           <  2:189700/251‑1 (MQ=255)
                                                                                                                                                                                                              cgcCAAACTGAACATCGACTGCCGCCTGCCAACCCTGGACCTGGACGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCTCGAAGCGGGGCTGAATGCGATGAAGTACGCCAACGACGGCGAAATCGCCGTTAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCGGTGATAACGGGGTTGGTATCGGTACACCGAACGAACACGTAGGTCATTATGGTCCGAATATCATGCGCGAAcgcg                         >  2:489292/1‑251 (MQ=255)
                                                                                                                                                                                                                       gACCCTCGACTGCCGGCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGGTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACg                  <  1:376436/249‑1 (MQ=255)
                                                                                                                                                                                                                                 ggCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGTTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCTCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCGAGGTGGGAc      <  2:284222/250‑1 (MQ=255)
                                                                                                                                                                                                                                    cGTCTGCCAACCATGGCCCCGGATGCGCCAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGCGCTGACTGCGATGAAGCACGCCAACGCCACCGAAATCGCCGTCAGCTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGAccct   >  2:573543/1‑248 (MQ=255)
                                                                                                                                                                                                                                     gTCTGCCACCCCTGGTAATGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGAGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGGCATTATGGTCTGAATATCATGAGCGAACGCGCGGAACGGCTAGTTGGGACGCTg  >  2:124055/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGACGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTg  >  2:468480/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTg  <  1:88716/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           
TTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTG  >  NC_000913/2586868‑2587347

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: