Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,290,728 T→A 6.4% D77E (GAT→GAA narP → response regulator in two‑component regulatory system with NarQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,290,7280TA 6.4% 209.8 / 5.1 63D77E (GAT→GAAnarPresponse regulator in two‑component regulatory system with NarQ
Reads supporting (aligned to +/- strand):  ref base T (34/24);  new base A (2/2);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

TACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCT  >  NC_000913/2290445‑2291003
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                   
tACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGt                                                                                                                                                                                                                                                                                >  1:65158/1‑289 (MQ=255)
         tacgctCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCa                                                                                                                                                                                                                                                                       <  2:271675/289‑1 (MQ=255)
         tacgctCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCa                                                                                                                                                                                                                                                                       <  1:327688/289‑1 (MQ=255)
               cAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGGCTGAAGCAACTCCTTTTCAGGTGATGATTTTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGTATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAattat                                                                                                                                                                                                                                                                 <  2:491349/289‑1 (MQ=255)
                        tGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGa                                                                                                                                                                                                                                                         >  1:70223/1‑288 (MQ=255)
                            tCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTa                                                                                                                                                                                                                                                    >  2:175076/1‑289 (MQ=255)
                            tCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTa                                                                                                                                                                                                                                                    >  2:26436/1‑289 (MQ=255)
                            tCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTa                                                                                                                                                                                                                                                    >  1:249613/1‑289 (MQ=255)
                            tCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTa                                                                                                                                                                                                                                                    >  1:196169/1‑289 (MQ=255)
                                  aaaaCCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGAt                                                                                                                                                                                                                                              <  2:34389/289‑1 (MQ=255)
                                      ccTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCt                                                                                                                                                                                                                                          <  2:243907/289‑1 (MQ=255)
                                          aGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCa                                                                                                                                                                                                                                       >  2:387225/1‑288 (MQ=255)
                                                      tGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGAAGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTgc                                                                                                                                                                                                                          >  2:406663/1‑289 (MQ=255)
                                                       gCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTgcg                                                                                                                                                                                                                         <  2:143021/289‑1 (MQ=255)
                                                           gAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGa                                                                                                                                                                                                                     >  2:118000/1‑289 (MQ=255)
                                                           gAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCGGGCGACTAGACTGGATATCGACGTGATCTTGCTGCACCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTAAATGCCTTGCGCAGGGTTGGCGTTACCGCGCAAATTATTATCCTGACCGTAACCACTGCCTCCAGCGATTTCCTTgacgtga                                                                                                                                                                                                                     >  1:485428/1‑283 (MQ=255)
                                                                  cacCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCg                                                                                                                                                                                                                                                <  1:391697/255‑1 (MQ=255)
                                                                  cacCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCg                                                                                                                                                                                                                                                >  2:391697/1‑255 (MQ=255)
                                                                  cacCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGc                                                                                                                                                                                                              >  2:93905/1‑289 (MQ=255)
                                                                  cacCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGc                                                                                                                                                                                                              >  2:86168/1‑289 (MQ=255)
                                                                        ttCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATcc                                                                                                                                                                                                                                                 <  2:389040/248‑1 (MQ=255)
                                                                        ttCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATcc                                                                                                                                                                                                                                                 >  1:389040/1‑248 (MQ=255)
                                                                                atgatTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCTCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGCACCAGACGGtt                                                                                                                                                                                                >  1:485026/1‑289 (MQ=255)
                                                                                       tgGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGAGGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCTCTGTTAGACGCAGGCGCCGACCGTTATATGtt                                                                                                                                                                                         >  1:314015/1‑289 (MQ=255)
                                                                                              tcatcCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGAc                                                                                                                                                                                  >  2:426784/1‑289 (MQ=255)
                                                                                                atcCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAg                                                                                                                                                                                <  2:70223/289‑1 (MQ=255)
                                                                                                      ttATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAcc                                                                                                                                                                           >  1:195986/1‑288 (MQ=255)
                                                                                                        aTGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCgg                                                                                                                                                                        <  1:387225/289‑1 (MQ=255)
                                                                                                                                         gATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGa                                                                                                                                       >  1:379204/1‑289 (MQ=255)
                                                                                                                                                               ccGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGGCGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCg                                                                                                                  <  1:175076/288‑1 (MQ=255)
                                                                                                                                                                          gACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAAt                                                                                                       <  2:195986/288‑1 (MQ=255)
                                                                                                                                                                            cGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCa                                                                                                     <  1:86168/288‑1 (MQ=255)
                                                                                                                                                                                         aTCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTacgtgaac                                                                                       <  1:93905/289‑1 (MQ=255)
                                                                                                                                                                                                            tGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGga                                                                    >  2:12858/1‑289 (MQ=255)
                                                                                                                                                                                                             ggATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGCGAACTGTTTGGCGCggca                                                                   >  1:319663/1‑287 (MQ=255)
                                                                                                                                                                                                               atatCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGAAGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGgaaga                                                                 >  1:394757/1‑289 (MQ=255)
                                                                                                                                                                                                                   cGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAagaaga                                                              >  2:237640/1‑288 (MQ=255)
                                                                                                                                                                                                                      cGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATc                                                            >  2:458830/1‑287 (MQ=255)
                                                                                                                                                                                                                           tCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAg                                                     <  1:26436/289‑1 (MQ=255)
                                                                                                                                                                                                                           tCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGAAGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAg                                                     <  1:406663/289‑1 (MQ=255)
                                                                                                                                                                                                                               gCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGt                                                  >  2:202488/1‑288 (MQ=255)
                                                                                                                                                                                                                                      ctcAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACgg                                          >  1:8513/1‑289 (MQ=255)
                                                                                                                                                                                                                                         aaTATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAgc                                       >  2:487034/1‑289 (MQ=255)
                                                                                                                                                                                                                                             tGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGa                                   >  1:346822/1‑289 (MQ=255)
                                                                                                                                                                                                                                               aaaGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGa                                   >  1:85413/1‑287 (MQ=255)
                                                                                                                                                                                                                                               aaaGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTCTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGc                                 >  2:373194/1‑289 (MQ=255)
                                                                                                                                                                                                                                                aaGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGc                                 <  2:8513/288‑1 (MQ=255)
                                                                                                                                                                                                                                                aaGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCt                                <  1:176632/289‑1 (MQ=255)
                                                                                                                                                                                                                                                   gTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGa                             >  2:162997/1‑289 (MQ=255)
                                                                                                                                                                                                                                                    tATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGAAGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTgg                              <  2:394757/287‑1 (MQ=255)
                                                                                                                                                                                                                                                     aTGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATg                           >  2:413287/1‑289 (MQ=255)
                                                                                                                                                                                                                                                     aTGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATg                           >  2:499777/1‑289 (MQ=255)
                                                                                                                                                                                                                                                      tGAGTGGCCTGGATACTCTCAATACCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGt                          >  1:87925/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                  aTACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGc               <  2:249613/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aTACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCt              <  1:499777/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCt              <  1:237640/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                      tctcAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGca          <  1:458830/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                      tctcAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGc           >  2:39660/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                       ctcAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGcac         <  1:207305/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                           aTGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTTAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAggg     >  1:510732/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                           aTGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAggg     >  2:476536/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                               cTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCt  <  2:65158/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCt  <  2:196169/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                   
TACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCTCTGAAGTGGTCGCCGAAGCGGGCGACGGCGCGAGCGCTATCGATCTGGCGAATAGACTGGATATCGACGTGATCTTGCTGGATCTCAATATGAAAGGTATGAGTGGCCTGGATACTCTCAATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGGAAGCGATTCGTGCCGGAGCGAAAGGCAGCAAAGTCTTTAGCGAACGCGTCAATCAGTACTTACGTGAACGTGAAATGTTTGGCGCGGAAGAAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCT  >  NC_000913/2290445‑2291003

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: