Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,218,509 T→A 5.1% T17S (ACG→TCG) ‡ yehY ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,218,5090TA 5.1% 334.3 / 3.9 113T17S (ACG→TCG) ‡yehYputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (54/53);  new base A (3/2);  total (58/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

Only 100 of 116 total aligned reads displayed.

GGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTTGCTGCAATG  >  NC_000913/2218360‑2218658
                                                                                                                                                     |                                                                                                                                                     
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)
                   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)
                          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)
                                    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                                                                                                                                       >  1:701919/1‑123 (MQ=255)
                                           aGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACagcagcagc                                                                                                                                       <  2:701919/123‑1 (MQ=255)
                                           aGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATaa                                                                                                            >  2:891419/1‑150 (MQ=255)
                                            gCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAg                                                                                                           <  2:1533331/150‑1 (MQ=255)
                                             cATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGt                                                                                                          >  1:1387823/1‑150 (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)
                                                         ttGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTa                                                                                              <  2:1713124/150‑1 (MQ=255)
                                                          tGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGCACAGGATAAATACGGAAATAAGTCACGGGAACTTAc                                                                                             >  2:1627214/1‑150 (MQ=255)
                                                          tGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTAc                                                                                             >  2:746114/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                                                                                          >  2:1535743/1‑150 (MQ=255)
                                                             ggCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTc                                                                                          <  1:863549/150‑1 (MQ=255)
                                                                  caGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCa                                                                                     <  1:1128404/150‑1 (MQ=255)
                                                                    gCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCAcc                                                                                   >  2:633900/1‑150 (MQ=255)
                                                                     cTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCAccc                                                                                  >  1:9313/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)
                                                                                  cGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAgg                                                                     >  2:1161152/1‑150 (MQ=255)
                                                                                   gCCCCTCACCCGAAACTAACCGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGt                                                                    <  2:68394/150‑1 (MQ=255)
                                                                                   gCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGt                                                                    <  1:785626/150‑1 (MQ=255)
                                                                                   gCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGt                                                                    <  1:119008/150‑1 (MQ=255)
                                                                                       cTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTc                                                                <  1:1460420/150‑1 (MQ=255)
                                                                                          aCCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGc                                                             >  2:494540/1‑150 (MQ=255)
                                                                                           cccGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTg                                                                            >  1:1092835/1‑134 (MQ=255)
                                                                                           cccGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTg                                                                            <  2:1092835/134‑1 (MQ=255)
                                                                                           cccGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCa                                                            >  1:918515/1‑150 (MQ=255)
                                                                                            ccGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAg                                                           >  2:556232/1‑150 (MQ=255)
                                                                                             cGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGc                                                          >  1:431477/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:941032/150‑1 (MQ=255)
                                                                                                                   tAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCt                                    >  1:1322719/1‑150 (MQ=255)
                                                                                                                   tAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCt                                    <  1:633900/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)
                                                                                                                               ggCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTAGTCCACCCTTTTTGTTTAAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATg                        >  2:409312/1‑150 (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:1715461/1‑150 (MQ=255)
                                                                                                                                    cgccgcTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGc                   >  1:788719/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)
                                                                                                                                          tGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAg             <  1:1263221/150‑1 (MQ=255)
                                                                                                                                             gATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGTAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGAAGCCACTTTATTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtg          >  1:1526206/1‑150 (MQ=255)
                                                                                                                                              aTTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgc         <  2:1387823/150‑1 (MQ=255)
                                                                                                                                                tGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgctg       <  1:1161152/150‑1 (MQ=255)
                                                                                                                                                tGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgctg       <  1:1008562/150‑1 (MQ=255)
                                                                                                                                                tGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgctg       <  1:24979/150‑1 (MQ=255)
                                                                                                                                                 gCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATCCGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgctgc      <  1:252243/150‑1 (MQ=255)
                                                                                                                                                 gCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgctga      <  1:176049/150‑2 (MQ=255)
                                                                                                                                                 gCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTAGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgctgc      >  1:1651804/1‑150 (MQ=255)
                                                                                                                                                   cGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGATGCTGCaa    >  1:275316/1‑150 (MQ=255)
                                                                                                                                                    gTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTTGCTGCAAt   >  2:599683/1‑150 (MQ=255)
                                                                                                                                                     tCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTTGCTGCAATg  <  2:214411/150‑1 (MQ=255)
                                                                                                                                                     tCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTTGCTGCAATg  <  2:1322719/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTTGCTGCAATG  >  NC_000913/2218360‑2218658

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: