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

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: