Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,744,088 A→T 12.6% Y211F (TAC→TTC)  mdtK → multidrug efflux system transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,744,0880AT12.6% 175.6 / 2.1 64Y211F (TAC→TTC) mdtKmultidrug efflux system transporter
Reads supporting (aligned to +/- strand):  ref base A (21/32);  new base T (7/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 6.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCGTTGTTGTGGGGCGCGCCGGGATATCTGTTCTTCCAGGTTGCCCGTAACCAGTGTGAAGGTCTGGCAAAAACCAAGCCGGGTATGGTAATGGGCTTTATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTCGCAGGACACCAGATTGCCCTGAACTTTAGTTCACTAATGTTCGTGCTTCCAATGT  >  NC_000913/1743851‑1744324
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                
gtgcgtTGTTGTGGGGCGCGCCGGGATATCTGTTCTTCCAGGTTGCCCGTAACCAGTGTGAAGGTCTGGCAAAAACCAAGCCGGGTATGGTAATGGGCTTTATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGGTATGTTCCTTGCCCCGGTTTTTTAc                                                                                                                                                                                                                                                 >  2:1236789/1‑239 (MQ=255)
 tgcgtTGTTGTGGGGCGCGCCGGGATATCTGTTCTTCCAGGTTGCCCGTAACCAGTGTGAAGGTCTGGCAAAAACCAAGCCGGGTATGGTAATGGGCTTTATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGGGGTTTTGGTCATTTTCCTTTCCCTTGGTTCTTTCa                                                                                                                                                                                                                                                >  2:985173/1‑239 (MQ=255)
       gttgtGGGGCGCGCCGGGATATCTGTTCTTCCAGGTTGCCCGTAACCAGTGTGAAGGTCTGGCAAAAACCAAGCCGGGTATGGTAATGGGCTTTATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTaaa                                                                                                                                                                                                                                           <  2:84973/238‑1 (MQ=255)
       gttgtGGGGCGCGCCGGGATATCTGTTCTTCCAGGTTGCCCGTAACCAGTGTGAAGGTCTGGCAAAAACCAAGCCGGGTATGGTAATGGGCTTTATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGGGGTTTTGGGCATGTTTCCTTCCCTTGGTTTTTTCCTTTaa                                                                                                                                                                                                                                            >  2:345029/1‑238 (MQ=255)
           tggtgtgctccGGGTTTTTTGTTCTTCCTGGTTGCCTGTAACCAGTGTGAAGGTGTGGCAAAAATTAAGCCGGGTATGGTAATGGGGTTTATCGGCCTGGTGGTGAACATCCCGGTTGACTATATCTTTATTTATGGTCATTCCGGTATGCCTTAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACgcgc                                                                                                                                                                                                                                      <  1:128377/230‑1 (MQ=255)
           gggggCGGCCGGGATATCGGGTTGTCGAGGTTGCCCGGAACCTGTGGGAAGGTCTGGCAAAAGCCAAGCCGGGTATGGTAATGGGCTTGATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGcc                                                                                                                                                                                                                                     <  1:799642/238‑1 (MQ=255)
                 gcgcCGGGATATCTGTTCTTCCAGGTTGCCCGTAACCAGTGTGAAGGTCTGGCAAAAACCAAGCCGGGTATGGTAATGGGCTTTATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCt                                                                                                                                                                                                                                 <  2:584183/238‑1 (MQ=255)
                  cgcCGGGATATCTGTTCTTCCAGGTTGCCCGTAACCAGTGTGAAGGTCTGGCAAAAACCAAGCCGGGTATGGTAATGGGCTTTATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTTCGGCGGGGGTTTGGGGTGTGTTCCTTGCCCTGGGTTTTTTCATTTAAAGCGCCCGccc                                                                                                                                                                                                                                >  2:789528/1‑236 (MQ=255)
                  cgcCGGGATATCTGTTCTTCCAGGTTGCCCGTAACCAGTGTGAAGGTCTGGCAAAAACCAAGCCGGGTATGGTAATGGGCTTTATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGGGTATTGGGGCATGTTTCTTTCCCTGGGTTTTTTCATTTAACCGGccccgtc                                                                                                                                                                                                                                >  2:920167/1‑234 (MQ=255)
                  cgcCGGGATATCTGTTCTTCCAGGTTGCCCGTAACCAGTGTGAAGGTCTGGCAAAAACCAAGCCGGGTATGGTAATGGGCTTTATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGGGGTTTGGGGTAAGTTCCTTGCCATGGGTTTTTCAATTAAACGCCCCCGCTc                                                                                                                                                                                                                                >  2:723947/1‑238 (MQ=255)
                             cTGTTCTTCCAGGTTGCCCGTAACCAGTGTGAAGGTCTGGCAAAAACCAAGCCGGGTATGGTAATGGGCTTTATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGatat                                                                                                                                                                                                                    <  2:103883/239‑1 (MQ=255)
                                ttcttcCAGGTTGCCCGTAACCAGTGTGAAGGTCTGGCAAAAACCAAGCCGGGTATGGTAATGGGCTTTATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGcc                                                                                                                                                                                                                                     <  2:1130408/219‑1 (MQ=255)
                                            tcccTTAACCATTGTGAAGGTTTGGCAAAAACCAAGCCGGGTATGGTAATGGGTTTTATCGGCCTGCGGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAgg                                                                                                                                                                                                     <  1:399674/238‑1 (MQ=255)
                                                        tgtgAAGGTCTGGCAAAAACCAAGCCGGGTATGGTAATGGGCTTTATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTTCGGCGGGGTATTTGGGCCTGTTTCCTGCCATGGGTTCTTAAATTTAACGCGCCCCCCCCCTGCGCCAAATTCGTAACGAAAAAGGGACCCCCAAAc                                                                                                                                                                                          >  2:904375/1‑238 (MQ=255)
                                                            aaGGTCTGGCAAAAACCAAGCCGGGTATGGTAATGGGCTTTATCGGCCTGCTGGTGAACATCCCGGAGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTTGGTGCTGTTCCTTGCCCTGGGTTCTTTCCTTTAACGCGCCCCCTCCATGCGCGGTATTCGTAACGAAAAAAGCACCGCCAAACCCGa                                                                                                                                                                                      >  2:166449/1‑238 (MQ=255)
                                                               gTCTGGCAAAAACCAAGCCGGGTATGGTAATGGGCTTTATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGGGGATTGGGTCATGTTCCTTGCCATTGTTTCTTACATTTAACGCGCCCGCTCCATTCGCGATATTCGGAACGAAAAAGGCAACGCAAAACCCCATcc                                                                                                                                                                                   >  2:735057/1‑238 (MQ=255)
                                                                 ttagaaaaaaCCACGCGGGGTAGGTTAAGGGGTTTTTTGGGGCTGCGGGGTAACATCCTGGTGAATTATATCTTTATTTATGGGCATTTCGGTATGCCTGAGCTCGGTGGGGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTg                                                                                                                                                                                 <  1:1236789/233‑1 (MQ=255)
                                                                                 cGGGTTTTGTAATGGGCTTTATCGGCCTGTGGGGGAACCTCCCGGTGAACTATATCTTTATTTATGGGCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGa                                                                                                                                                                   <  1:904375/236‑1 (MQ=255)
                                                                                                    tATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGGGTATTGTGGCATGTTCCTTTGCCTGGGTTCTTTCAATTAACGCGCCCGCCCCATTCGCCATATTCGGTAGCAAAAAGTCACCCCCAAACCCCCTCCCGCGGTCTTGTAAAAACCGATTCAAccccgtgttac                                                                                                                                               >  1:922402/1‑227 (MQ=255)
                                                                                                    tATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCAGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGGGGATTGGGTCCAGTTTCTTGCCCAGGTTTCTTACAATAAAAGAGCCCGCTCCATTCGCGCTATTCCTAACGGAAAAAGTAACCCCAAACCCGAACCCGCGGTTAGTAAACGACAGATTaaaaaccatatgccgg                                                                                                                                             >  1:491532/1‑226 (MQ=255)
                                                                                                     gtCGGTCTGCTTGTGAACTTCCCGGTGAACTTTATCTGTATTTATGGTCCTTGCGGTAGGTCTGAGTTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGAt                                                                                                                                            <  1:789528/238‑1 (MQ=255)
                                                                                                              ctgGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGGTTCTTTCATTAAACGCGCCCGCTCCATGCGCGGGATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGGTTTGAAACGACTGATTTACATCGGGTTTCCCGTTTCGCTggg                                                                                                                                  >  2:1030586/1‑238 (MQ=255)
                                                                                                                  tGAACATCCCGGTGAACTATATCTTTATTTATGGGCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGGGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACt                                                                                                                               <  2:922402/238‑1 (MQ=255)
                                                                                                                       aTCCCGGTTAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCtt                                                                                                                        <  2:1014461/239‑1 (MQ=255)
                                                                                                                                tcagatatCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGa                                                                                                               <  2:1113648/234‑1 (MQ=255)
                                                                                                                                              ttatGGTCATTTCGGTATGCCCGAGCGCGGTGGCGTTGGTTGTGGCGTGGCTTCTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCcgtcg                                                                                                <  2:311810/239‑1 (MQ=255)
                                                                                                                                               tatGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTTGGGCATGGTCCTTGCCATGGGTTCTTTAATTAAATGCGCCCGCTCCATGCGCCATTTTCGTAACGGAAAAGGCACCGCAAAAACCGCTCCTTCGGCTATGAAATGAACTATTCACACCGGTTTTCCCAGTTCGCGGGGCCTGTGTCTTTGAAGTAAACCTGTTTGcccg                                                                                                  >  1:1113648/1‑235 (MQ=255)
                                                                                                                                               tatGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCcgtcgt                                                                                               <  2:693458/239‑1 (MQ=255)
                                                                                                                                                                  cTGGGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCg                                                                            <  2:208414/239‑1 (MQ=255)
                                                                                                                                                                   tGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGc                                                                           <  2:856931/239‑1 (MQ=255)
                                                                                                                                                                    ggtcccGGTGGTTTTGGTTGTGGGGTGGCTTCTGGGGCGGTGTTTTGGGGACTGTTCCTTGCCATGGTTTCTTACATTAAACGCACCCGCGCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGTTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCt                                                                          <  1:985173/234‑1 (MQ=255)
                                                                                                                                                                    gAGCTTGGTGGGTTTGGTTGGGGGGTGGCTTCTGGGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCt                                                                          <  1:166449/239‑1 (MQ=255)
                                                                                                                                                                    gAGCTCGGTGGCGTTGGTTTGGGCTTGGCTACTGCGGCTGTGTATTGGGTCATGTTCCTTTCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCt                                                                          <  2:326739/239‑1 (MQ=255)
                                                                                                                                                                       cTCGGGGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCcgtcgt                                                                                               <  2:602735/215‑1 (MQ=255)
                                                                                                                                                                                 ttggttggGGCGTGGCTACTGGGGCGGTGTTTTGGGTCATGTTCCTTGCCATGGGTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGGTTGCCGTCGTGGCTCTGTTAGGGTCTCCGCTCGGTATTGTTGa                                                              <  2:621739/238‑1 (MQ=255)
                                                                                                                                                                                 ttggttggGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGa                                                              <  2:557060/238‑1 (MQ=255)
                                                                                                                                                                                 ttggttgTGGCGTGGCTACTGGGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGa                                                              <  2:526364/238‑1 (MQ=255)
                                                                                                                                                                                             tggcTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTCGCAGGaca                                                  <  2:491532/238‑1 (MQ=255)
                                                                                                                                                                                                 tACTGTGGCGGTTTTTTGGGTCATTTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGGTCCATGCGCGATATTCGTAACGAAAAAGGCCCCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTCGCAGGACACCAGa                                            <  1:735057/239‑1 (MQ=255)
                                                                                                                                                                                                  gatggtgataagggacaggTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTCGCAGGACACCAGa                                            <  2:576904/221‑1 (MQ=255)
                                                                                                                                                                                                    tGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTAttgtt                                                                <  1:921825/217‑1 (MQ=255)
                                                                                                                                                                                                    tGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGGAttgtt                                                                >  2:921825/1‑217 (MQ=255)
                                                                                                                                                                                                         cggTGTCTTGGGGCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAAGCAACGCAAAACCCGATCCTGCCGGTATGAAACGACCGAATTAAATCGGTTTGCCGATTTCGCCGGGCACTGTACTTGAAGATGAACTCGTTTCCGTCGTGGCGCTGGTAGGGTCTCCGCCCGGGGCTGGTGATGTTCAGGAGACCCAGATTGccccg                                    >  1:398763/1‑237 (MQ=255)
                                                                                                                                                                                                                    gTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTCGCAGGACACCAGAt                                           <  2:595642/221‑1 (MQ=255)
                                                                                                                                                                                                                     gtctgtttCCTTGCCCTGGTTTCTTTATTTAAACGCGCCCGTTCCATGCGCGATATTCGTAACGAAAAAGGCACCGTAAAACCCGTGCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTCGCAGGACACCAGATTGCCCTGAACTTTAGtt                          <  1:924304/232‑1 (MQ=255)
                                                                                                                                                                                                                      cATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCt                                                                          <  2:1318548/189‑1 (MQ=255)
                                                                                                                                                                                                                      cATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCt                                                                          >  1:1318548/1‑189 (MQ=255)
                                                                                                                                                                                                                      cATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTCGCAGGACCCCAGATTGCCCTGAACTTTAGTTc                         >  2:82226/1‑237 (MQ=255)
                                                                                                                                                                                                                       aTGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTcc                                                         >  1:867721/1‑204 (MQ=255)
                                                                                                                                                                                                                       aTGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTCg                                                         <  2:867721/205‑1 (MQ=255)
                                                                                                                                                                                                                        tGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTATTGTATCCGCTCGGTATTGTTGATGTCGCAGGCCCCCAGATTGCCCTGAACTTTAGTTCCCTa                     >  2:99792/1‑239 (MQ=255)
                                                                                                                                                                                                                        tGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTCGCAGGACACCAGATTGCCCTGAACTTTAGTTCACTa                     >  2:1123043/1‑239 (MQ=255)
                                                                                                                                                                                                                        tGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTCGCAGGACACCAGATTGCCCTGAACTTTAGTTCACTa                     >  2:601605/1‑239 (MQ=255)
                                                                                                                                                                                                                        tGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACGGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTCGCAGGACACCAGATTGCCCTGAACTTTAGTTCACTa                     >  2:309595/1‑239 (MQ=255)
                                                                                                                                                                                                                              cTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGGGTGCCCGCTCGGTATTGTTGATTTCGCAGGGCACCAGATTTCCCTGAACTTTCGTTCATAAATGtt                >  1:92129/1‑239 (MQ=255)
                                                                                                                                                                                                                               ttGCCATGGTTTCTTACATTAAACGCGCCCGCCCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCCCTTTTAAGGTCTCCGCTCGGGATTGTTGATTTCCCCGGACACCAGATTGCCCTGaaacctggtccttaaaacac                >  1:595642/1‑220 (MQ=255)
                                                                                                                                                                                                                               ctgcgtgtgTTTTTTTCATTAAGCGGAGCCTCTCCGTGCGCGCTAGTCGTAACGAAAAAGGGCCGCCAAAACCCGTTCCTGCGGGTATGAAGCGGCTGATTCATGTCTGTTTGCAGATTGCGCTGTCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTTTCTCCTCTCGGTATTGTTGATGTCGCAGGACACCAGATTGCCATGAACTCTAGTTCACTatcgtta               <  1:798951/231‑7 (MQ=255)
                                                                                                                                                                                                                                tGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTCGCAGGACACCAGATTGCCCTGAACTTTAGTTCCCTAATGt                 >  2:1231058/1‑236 (MQ=255)
                                                                                                                                                                                                                                tGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGGGCAACTGTTTGCCGCTGTGGCGCTGTTAGTGTCTCCGCCCGGTATTGTTGATGTCGCAGGCAAACAATTTGCCCCGAAcctgtcctttatacaac                 >  1:576904/1‑220 (MQ=255)
                                                                                                                                                                                                                                   cATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGc                                                                                                                                  <  1:32166/123‑1 (MQ=255)
                                                                                                                                                                                                                                   cATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGc                                                                                                                                  >  2:32166/1‑123 (MQ=255)
                                                                                                                                                                                                                                    aTGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTCGCAGGACACCAGATTGCCCTGAACTTTAGTTCACTAATGTTCGTGc           >  2:124499/1‑238 (MQ=255)
                                                                                                                                                                                                                                     tGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAACACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCACTGTTAGTGTCCCCTCTCGGTATTGTTGGTGCCGCAGGACCCCACATTGCCCCGAAATTTTGTTCACTAATATTCGTGCt          >  1:381827/1‑238 (MQ=255)
                                                                                                                                                                                                                                             tACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTGGCAGGACACCAGATTGCCCTGAACTTTAGTTCACTAATGTTCGTGCTTCCAATGt  >  2:732013/1‑238 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                
GTGCGTTGTTGTGGGGCGCGCCGGGATATCTGTTCTTCCAGGTTGCCCGTAACCAGTGTGAAGGTCTGGCAAAAACCAAGCCGGGTATGGTAATGGGCTTTATCGGCCTGCTGGTGAACATCCCGGTGAACTATATCTTTATTTATGGTCATTTCGGTATGCCTGAGCTCGGTGGCGTTGGTTGTGGCGTGGCTACTGCGGCGGTGTATTGGGTCATGTTCCTTGCCATGGTTTCTTACATTAAACGCGCCCGCTCCATGCGCGATATTCGTAACGAAAAAGGCACCGCAAAACCCGATCCTGCGGTTATGAAACGACTGATTCAACTCGGTTTGCCGATTGCGCTGGCACTGTTCTTTGAAGTGACACTGTTTGCCGTCGTGGCTCTGTTAGTGTCTCCGCTCGGTATTGTTGATGTCGCAGGACACCAGATTGCCCTGAACTTTAGTTCACTAATGTTCGTGCTTCCAATGT  >  NC_000913/1743851‑1744324

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