Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,218,502 A→G 5.8% I19T (ATC→ACC) ‡ yehY ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,218,5020AG 5.8% 303.7 / 6.8 103I19T (ATC→ACC) ‡yehYputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (51/46);  new base G (4/2);  total (55/48)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 105 total aligned reads displayed.

TTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTTGC  >  NC_000913/2218354‑2218651
                                                                                                                                                    |                                                                                                                                                     
tttACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGAt                                                                                                                                                      <  1:851667/150‑1 (MQ=255)
      gggAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCAACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGCCGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCTCCGCTGCGATTGCCGt                                                                                                                                                >  2:800299/1‑150 (MQ=255)
       ggAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTc                                                                                                                                               <  2:913645/150‑1 (MQ=255)
             aaaCAGGCCGTCAGCCAGGCGCAACCAACGCCCCCCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAAc                                                                                                                                         >  2:1299375/1‑150 (MQ=255)
             aaaCAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAAc                                                                                                                                         <  2:67755/150‑1 (MQ=255)
                cAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAAcagc                                                                                                                                      <  2:1152610/150‑1 (MQ=255)
                cAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCAATTGCCGTCAGCAAcagc                                                                                                                                      <  2:530957/150‑1 (MQ=255)
                cAGGCCGTCAGCAAGGCGCAACCAACGCCCACCAGCATCCAGTTCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTATACGTTTAGGCGCGTAACTGATAAACGGCAGCGACGCTGCGATTGCCGTCAGCAAcagc                                                                                                                                      <  2:1608824/150‑1 (MQ=255)
                  ggCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAAcagcag                                                                                                                                    <  1:10254/150‑1 (MQ=255)
                  ggCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAAcagcac                                                                                                                                    >  2:1008562/1‑149 (MQ=255)
                    ccGTCAGCCAGGCGCAACCAACGCCCACCCGCATCCAGATCGTTTGCGGCCACAGCTGACAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAACCGGCAGCGCCGCTGCGATTGCCGTCAGCACcagcagca                                                                                                                                  <  2:924017/150‑1 (MQ=255)
                       tCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCTGCTGCTGTTGCTGACGGCAATCGCAgcggcg                                                                                                                               <  2:1672783/150‑1 (MQ=38)
                         aGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGGTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGTGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGcc                                                                                                                             <  1:172792/150‑1 (MQ=255)
                          gCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCa                                                                                                                            >  1:1035712/1‑150 (MQ=255)
                          gCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAg                                                                                                                           >  1:1360416/1‑150 (MQ=255)
                           ccAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCTGCTGCTGTTGCTGACGGCAATCGCAGCGgcgctgc                                                                                                                           <  2:140717/150‑4 (MQ=38)
                            cAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGCAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCCGCGCCAGa                                                                                                                          >  1:439101/1‑150 (MQ=255)
                                cgcAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAg                                                                                                                      >  1:1072971/1‑150 (MQ=255)
                                 gcAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCCCCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAgg                                                                                                                     >  2:1263221/1‑150 (MQ=255)
                                   aaccaacGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGAt                                                                                                                   >  1:523315/1‑150 (MQ=255)
                                        acGCCCACCAGCATCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATAc                                                                                                             <  2:1559167/150‑1 (MQ=255)
                                         cGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATAc                                                                                                             >  2:1460420/1‑150 (MQ=255)
                                          gCCCACCAGCATCCAGATCGTTTGCGCCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACg                                                                                                            <  2:220859/150‑1 (MQ=255)
                                             caccaGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGaa                                                                                                         >  1:1313153/1‑150 (MQ=255)
                                               ccaGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAAt                                                                                                       <  2:1017008/150‑1 (MQ=255)
                                                caGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAAAAGCAGCAGCGCCAGAACAGGATTAATACGGAAATa                                                                                                      <  1:1486793/150‑1 (MQ=255)
                                                 aGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACagcagcagc                                                                                                                                <  2:701919/123‑1 (MQ=255)
                                                 aGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACagcagcagc                                                                                                                                >  1:701919/1‑123 (MQ=255)
                                                 aGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATaa                                                                                                     >  2:891419/1‑150 (MQ=255)
                                                  gCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAg                                                                                                    <  2:1533331/150‑1 (MQ=255)
                                                      ccAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCAc                                                                                                >  2:1183631/1‑150 (MQ=255)
                                                       cAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACg                                                                                               >  1:427790/1‑150 (MQ=255)
                                                        aGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACgg                                                                                              <  1:1099295/150‑1 (MQ=255)
                                                        aGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACgg                                                                                              >  1:1256561/1‑150 (MQ=255)
                                                        aGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGAAAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACgg                                                                                              >  2:5975/1‑150 (MQ=255)
                                                             gTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATCCGGAAATAAGTCACGGGAAAt                                                                                         >  1:1330732/1‑150 (MQ=255)
                                                               ttGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTa                                                                                       <  2:1713124/150‑1 (MQ=255)
                                                                tGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGCACAGGATAAATACGGAAATAAGTCACGGGAACTTAc                                                                                      >  2:1627214/1‑150 (MQ=255)
                                                                tGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTAc                                                                                      >  2:746114/1‑150 (MQ=255)
                                                                tGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTAc                                                                                      >  2:752375/1‑150 (MQ=255)
                                                                 gCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACt                                                                                     <  1:1228931/150‑1 (MQ=255)
                                                                  gcgCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACtt                                                                                    <  2:1360416/148‑1 (MQ=255)
                                                                   ggCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTc                                                                                   <  2:867016/150‑1 (MQ=255)
                                                                   ggCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTc                                                                                   <  1:863549/150‑1 (MQ=255)
                                                                   ggCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTc                                                                                   >  2:1535743/1‑150 (MQ=255)
                                                                        caGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCa                                                                              <  1:1128404/150‑1 (MQ=255)
                                                                          gCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCAcc                                                                            >  2:633900/1‑150 (MQ=255)
                                                                           cTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCAccc                                                                           >  1:9313/1‑150 (MQ=255)
                                                                             gCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCtt                                                                         <  2:1072971/150‑1 (MQ=255)
                                                                                agagATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCttttt                                                                      >  1:1528374/1‑150 (MQ=255)
                                                                                   taTGACGCCCCTCACCCGAAACTAAACGATTAGGCTCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGtt                                                                   <  1:1308827/149‑1 (MQ=255)
                                                                                   gaTGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGtt                                                                   <  1:1567753/150‑1 (MQ=255)
                                                                                   gaTGACGCCCCTCACCAGAAACTAAACGATTAGTCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGTACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGtt                                                                   <  1:256797/150‑1 (MQ=255)
                                                                                      gACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCa                                                                >  2:1154780/1‑150 (MQ=255)
                                                                                        ctccccTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCCGCGCCGCTGCGATTGCCGTCAGCTCCCGCAGCAGCGCCCGCACAGGATTAATACGGAAATACGTCACGGGAAATTACTACGTCCCCCCTTTTTGTTACAgg                                                              >  1:46744/3‑150 (MQ=255)
                                                                                        cGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAgg                                                              <  1:800299/150‑1 (MQ=255)
                                                                                        cGCCCCTCACCCGAAACTAAACGATTAGGCGAGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCATCATCATCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAgg                                                              <  1:635145/150‑1 (MQ=255)
                                                                                         gCCCCTCACCCGAAACTAACCGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGt                                                             <  2:68394/150‑1 (MQ=255)
                                                                                         gCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGt                                                             <  1:119008/150‑1 (MQ=255)
                                                                                         gCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGt                                                             <  1:785626/150‑1 (MQ=255)
                                                                                             cTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTc                                                         <  1:1460420/150‑1 (MQ=255)
                                                                                                aCCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGc                                                      >  2:494540/1‑150 (MQ=255)
                                                                                                 cccGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTg                                                                     <  2:1092835/134‑1 (MQ=255)
                                                                                                 cccGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTg                                                                     >  1:1092835/1‑134 (MQ=255)
                                                                                                 cccGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCa                                                     >  1:918515/1‑150 (MQ=255)
                                                                                                  ccGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAg                                                    >  2:556232/1‑150 (MQ=255)
                                                                                                   cGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGc                                                   >  1:431477/1‑150 (MQ=255)
                                                                                                          aaaCGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACtttt                                            >  1:214411/1‑150 (MQ=255)
                                                                                                            aCGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACtttttt                                          >  1:202229/1‑150 (MQ=255)
                                                                                                             cGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTg                                         >  2:941032/1‑150 (MQ=255)
                                                                                                                ttACGCGCGTCACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCTCGGGAAATTCCTGCGTCCACCCTTTTTGTTGCAGGTAGGCGGCAGCCACTTTTTTGGCa                                      <  2:46744/150‑1 (MQ=255)
                                                                                                                  aGGCGCGTAACTGATAAACGGCCGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATc                                    <  1:1440355/150‑1 (MQ=255)
                                                                                                                  aGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATc                                    >  1:470393/1‑150 (MQ=255)
                                                                                                                   ggCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTACTACGGAAATAAGTCCCGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATcc                                   <  1:1299375/150‑1 (MQ=255)
                                                                                                                       cgTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTc                               <  1:1183631/150‑1 (MQ=255)
                                                                                                                         tAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCt                             >  1:1322719/1‑150 (MQ=255)
                                                                                                                         tAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCt                             <  1:941032/150‑1 (MQ=255)
                                                                                                                         tAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCt                             <  2:427790/150‑1 (MQ=255)
                                                                                                                            cTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTcc                          <  2:439101/150‑1 (MQ=255)
                                                                                                                              gATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCcac                        >  1:951521/1‑150 (MQ=255)
                                                                                                                                  aaCGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCa                    >  1:618946/1‑150 (MQ=255)
                                                                                                                                   aCGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCaa                   >  1:1180915/1‑150 (MQ=255)
                                                                                                                                     ggCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATg                 <  1:494540/150‑1 (MQ=255)
                                                                                                                                      ctggcgcTGCTGCTGTTGCTGACGGCAATCGCAGCGGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGc                >  1:1548116/4‑150 (MQ=38)
                                                                                                                                       cAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCt               <  2:1528374/150‑1 (MQ=255)
                                                                                                                                        aGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTg              >  2:291870/1‑150 (MQ=255)
                                                                                                                                         gcgcCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTgg             <  1:1154780/150‑1 (MQ=255)
                                                                                                                                          cgctgcTGCTGTTGCTGACGGCAATCGCAGCGGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGc            >  2:1009556/1‑150 (MQ=38)
                                                                                                                                          cgctgcTGCTGTTGCTGACGGCAATCGCAGCGGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCAACCAGTCCTTCCACAGCAATGCTGGc            >  1:706075/1‑150 (MQ=38)
                                                                                                                                          cgccgcTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGc            >  1:788719/1‑150 (MQ=255)
                                                                                                                                          cgccgcTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGc            >  1:1715461/1‑150 (MQ=255)
                                                                                                                                           gctgctgctGTTGCTGCCGGCAATCGCAGCGGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCa           >  2:1160059/4‑150 (MQ=255)
                                                                                                                                           gccgcTGCGATTGCCGTCAGCAACAGCAGCAGCGCAAGAACAGGATTAATAAGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATACATTCATTCCACAGCAATGCTGGCa           >  2:176049/1‑150 (MQ=255)
                                                                                                                                             cgcTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCAtt         >  2:475036/1‑150 (MQ=255)
                                                                                                                                              gcTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTc        >  2:379039/1‑150 (MQ=255)
                                                                                                                                               cTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCAGGGGAAATTACTTCGTCCTCCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCGTTCCACAGCATTGCTGGCATTCa       >  1:1133550/1‑150 (MQ=255)
                                                                                                                                                tGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAg      <  1:1263221/150‑1 (MQ=255)
                                                                                                                                                   gATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGTAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGAAGCCACTTTATTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtg   >  1:1526206/1‑150 (MQ=255)
                                                                                                                                                    aTTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgc  <  2:1256561/150‑1 (MQ=255)
                                                                                                                                                    aTTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgc  <  2:1387823/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTTGC  >  NC_000913/2218354‑2218651

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