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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,587,1120CT96.6% 176.8 / NA 59A461V (GCA→GTA) narQsensor histidine kinase NarQ
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (22/35);  total (24/35)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑CACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAG  >  NC_000913/2586866‑2587336
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                  
acttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACt                                                                                                                                                                                                                                 <  1:159459/250‑1 (MQ=255)
acttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTg                                                                                                                                                                                                                                <  1:81750/251‑1 (MQ=255)
    acGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGAt                                                                                                                                                                                                                             <  1:262526/250‑1 (MQ=255)
      gTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATCCGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATgcg                                                                                                                                                                                                                          >  1:111662/1‑251 (MQ=255)
      gTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTGTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATgcg                                                                                                                                                                                                                          <  1:287301/251‑1 (MQ=255)
       tATCCAGATGACGTTACTGCAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCCTTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATgcgc                                                                                                                                                                                                                         >  2:170118/1‑251 (MQ=255)
        aTCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGCCTCCTTTCCGCCTGACGCTGCAGCAGGCGGGTCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCACACTGACCCTCGACTGCCGGCCGCCAACCCTGG‑TACTGGACgcgc                                                                                                                                                                                                                         >  2:107272/1‑250 (MQ=255)
          ccAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGATCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACGGGATGCGCaaa                                                                                                                                                                                                                      <  2:54191/251‑1 (MQ=255)
              ttGACGTCACTGCTGCGTTCGCTACAGGAAGATAACGCCCCCGCACACCGTCTCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGAAAGTTACCCGAGCTGTTGACTACTTTACCCCTGACGCTGCAGCAGGCGGATCTCCCCTACGCCTTGAGGGAAATGCTGGAAACGTTACCAACTCAAACCCACGCCAAACTGACCCTAGACTGCCGTCTGCCAACCCTGG‑TCCTGGGTGCTCAAATGc                                                                                                                                                                                                                   >  2:198888/1‑250 (MQ=255)
                     tACTGAAGCGTTCGCTACCGGAAGATAAAGCCACCGCACAAAGTATCATGGCCGATTTTTACCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTAAAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCACACCCTGGTCA‑TGGATGCGCAAATGCAGGTGCA‑t                                                                                                                                                                                                          >  2:376135/1‑251 (MQ=255)
                           aGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑Tttgttg                                                                                                                                                                                                    <  2:172139/251‑1 (MQ=255)
                                   aTACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAattat                                                                                                                                                                                            >  1:184609/1‑251 (MQ=255)
                                          aaGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGaa                                                                                                                                                                                     <  1:305458/251‑1 (MQ=255)
                                                    cACCGCACAAAGTATCATGTCCGATGTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGCAAATGCTGTATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCTACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGa                                                                                                                                                                          <  1:361425/251‑1 (MQ=255)
                                                           caAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGAt                                                                                                                                                                   <  1:250759/251‑1 (MQ=255)
                                                                         ggATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTAAGCGCGCTGTTGCCTCCTTTCCGCCTGCCGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGACATGCTGGTTACGTTACCAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGGCTGCACACCCTGG‑CACTGGATGCGCAAATGCAGG‑GAACTTTTTTGCAAATTATTCGCGAAGCGCGT‑CTGAAGGCGATGAAGCacgccaacg                                                                                                                                                     >  2:216948/2‑251 (MQ=255)
                                                                          gATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCacgccaacgc                                                                                                                                                    >  1:334124/1‑251 (MQ=255)
                                                                             ttttCCTAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATGATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCacgccaacgcca                                                                                                                                                  <  1:354047/250‑1 (MQ=255)
                                                                              tttCCCAGGCGTGGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATGATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACTCCAACGCCAg                                                                                                                                                 <  2:375921/250‑1 (MQ=255)
                                                                                    aGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATc                                                                                                                                          >  1:172839/1‑251 (MQ=255)
                                                                                     ggCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCg                                                                                                                                         >  1:216072/1‑251 (MQ=255)
                                                                                           gAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACGTCTTTCCGCCTGACGCTGCAGCAGTCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCTGTCGGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCa                                                                                                                                   <  1:262972/251‑1 (MQ=255)
                                                                                              tgatgCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGGGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGtt                                                                                                                                >  1:328738/1‑251 (MQ=255)
                                                                                                  gCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTAGTCGCGAAGCG‑GTGCTGAATGCTATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCg                                                                                                                             <  2:178184/250‑1 (MQ=255)
                                                                                                  gCTTATCGGCAGTTACGCGAGCTGTGGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCCCCGCATTGAGGGAAAGGCTGGATACGTTACAAAAGCAAACCAGCGCCAAACTGATCCTCGACTGCCGTCTGCCACCCCTGG‑TACTGGAGGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGTACGCCAACGCCAGCGAAATCGCCGTCAGTTGCg                                                                                                                             <  1:174372/250‑1 (MQ=255)
                                                                                                       tCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACcg                                                                                                                       >  1:194468/1‑251 (MQ=255)
                                                                                                       tCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCAAAACTGACCCTCGACTGCCGTCTGCCAACACTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCCGAATGCGATGAAGCACGCCAACGCCAGCGACATCGCCGTAAGTTGCGTCACcg                                                                                                                       >  1:183744/1‑251 (MQ=255)
                                                                                                         ggCAGTTACGCGAGCTGTGGACTACTTTCCGCCTGAAGCTGCAGCAGGCGTATCGCCCCTCCGCATGGAGGGAAATGCTGGATACGTTACAAAAGCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCAAACCCTGG‑TAATGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑TTGCTGAATGAGATGAAGCACGCAAACGCCAGTGAAATCGCCTGCAGTCGAGTCACcgcg                                                                                                                     <  2:236313/251‑1 (MQ=255)
                                                                                                           cAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGcc                                                                                                                   <  1:170118/251‑1 (MQ=255)
                                                                                                              ttACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGCAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCgg                                                                                                                 >  1:198353/1‑250 (MQ=255)
                                                                                                               tACGCGAGCTGTTGACAACTTTCCGCCTGACGCTGAAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCggcc                                                                                                               >  1:100497/1‑249 (MQ=255)
                                                                                                                 cgcgAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCCGTTGCGTCACCGCGCCGGACgg                                                                                                             >  2:338289/1‑251 (MQ=255)
                                                                                                                 cgcgAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACgg                                                                                                             >  1:200802/1‑251 (MQ=255)
                                                                                                                    gAGCTGTTGACTTCTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCaa                                                                                                          <  2:145808/251‑1 (MQ=255)
                                                                                                                    gAGCTGTTGACTACTTACCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATGCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCaa                                                                                                          <  1:198888/251‑1 (MQ=255)
                                                                                                                      gCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATc                                                                                                        <  1:338289/251‑1 (MQ=255)
                                                                                                                         gTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATccca                                                                                                     <  1:16542/251‑4 (MQ=255)
                                                                                                                            gACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑Tttg                                                                                                                                                                                                       <  2:115747/151‑1 (MQ=255)
                                                                                                                            gACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑Tttg                                                                                                                                                                                                       >  1:115747/1‑151 (MQ=255)
                                                                                                                                actTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTta                                                                                              <  1:376135/251‑1 (MQ=255)
                                                                                                                                actTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCAATGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGttt                                                                                               <  2:194468/250‑1 (MQ=255)
                                                                                                                                    tCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTatat                                                                                           <  1:158200/250‑1 (MQ=255)
                                                                                                                                       gCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCCGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGt                                                                                       >  1:203056/1‑251 (MQ=255)
                                                                                                                                           gACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGAt                                                                                    <  1:229501/250‑1 (MQ=255)
                                                                                                                                             cGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAACTCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATacc                                                                                 >  2:358517/1‑249 (MQ=255)
                                                                                                                                                     cagGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATAGCCGTGATAACGGGATTgg                                                                         <  2:146903/251‑1 (MQ=255)
                                                                                                                                                           gATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATACCGGGATTGGTATCgg                                                                   >  2:84409/1‑251 (MQ=255)
                                                                                                                                                                               ggAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAg                                               <  2:302584/251‑1 (MQ=255)
                                                                                                                                                                                            aCGTTACAAAAGCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTGGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGGCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTc                                   <  2:372956/250‑1 (MQ=255)
                                                                                                                                                                                                    aaaTCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATc                           >  2:90387/1‑250 (MQ=255)
                                                                                                                                                                                                        cAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATgcg                      <  1:90387/251‑1 (MQ=255)
                                                                                                                                                                                                        cAAACCAGCGCCAAACTGACCCTCGACTGCCGGCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATAGCATgc                       <  2:184609/250‑1 (MQ=255)
                                                                                                                                                                                                         aaaCCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATgcgc                     >  2:126596/1‑251 (MQ=255)
                                                                                                                                                                                                                cgcCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCTCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAAcgcg              <  1:401719/251‑1 (MQ=255)
                                                                                                                                                                                                                   cAAACTGCCCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATGATGGTCTGAATATCATGCGCGAACGCGCgg           <  2:73822/251‑1 (MQ=255)
                                                                                                                                                                                                                        tGACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACgg      >  1:386403/1‑251 (MQ=255)
                                                                                                                                                                                                                         gACCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACg       <  2:97495/249‑1 (MQ=255)
                                                                                                                                                                                                                          aCCCTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAAGTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCt    <  1:251530/251‑1 (MQ=255)
                                                                                                                                                                                                                            ccTCGACTGCCGTCTGCCAACCCTGG‑TACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAg  <  1:56740/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                  
ACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGG‑CACTGGATGCGCAAATGCAGGTGCA‑TTTGTTGCAAATTATTCGCGAAGCG‑GTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAG  >  NC_000913/2586866‑2587336

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