Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,218,511 A→G 5.2% L16P (CTG→CCG)  yehY ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,218,5110AG 5.2% 334.4 / 6.8 115L16P (CTG→CCG) yehYputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (55/54);  new base G (4/2);  total (59/56)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-02

Only 100 of 115 total aligned reads displayed.

AAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTTGCTGCAATGT  >  NC_000913/2218367‑2218659
                                                                                                                                                |                                                                                                                                                    
aaaCAGGCCGTCAGCCAGGCGCAACCAACGCCCCCCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAAc                                                                                                                                                 >  2:1299375/1‑150 (MQ=255)
aaaCAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAAc                                                                                                                                                 <  2:67755/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)
              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)
                            cGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATAc                                                                                                                     >  2:1460420/1‑150 (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)
                                    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)
                                           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:752375/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                                                                                           <  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)
                                                                      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)
                                                                           cGCCCCTCACCCGAAACTAAACGATTAGGCGAGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCATCATCATCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAgg                                                                      <  1:635145/150‑1 (MQ=255)
                                                                            gCCCCTCACCCGAAACTAACCGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGt                                                                     <  2:68394/150‑1 (MQ=255)
                                                                            gCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGt                                                                     <  1:119008/150‑1 (MQ=255)
                                                                                cTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTc                                                                 <  1:1460420/150‑1 (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)
                                                                                             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)
                                                                                                      ggCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTACTACGGAAATAAGTCCCGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATcc                                           <  1:1299375/150‑1 (MQ=255)
                                                                                                          cgTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTc                                       <  1:1183631/150‑1 (MQ=255)
                                                                                                            tAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCt                                     <  2:427790/150‑1 (MQ=255)
                                                                                                            tAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCt                                     <  1:633900/150‑1 (MQ=255)
                                                                                                            tAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCt                                     >  1:1322719/1‑150 (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)
                                                                                                                                       aTTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgc          <  2:1256561/150‑1 (MQ=255)
                                                                                                                                         tGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgctg        <  1:1008562/150‑1 (MQ=255)
                                                                                                                                         tGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgctg        <  1:24979/150‑1 (MQ=255)
                                                                                                                                         tGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgctg        <  1:1161152/150‑1 (MQ=255)
                                                                                                                                         tGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgctg        <  1:867595/150‑1 (MQ=255)
                                                                                                                                          gCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgctga       <  1:176049/150‑2 (MQ=255)
                                                                                                                                          gCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTAGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTtgctgc       >  1:1651804/1‑150 (MQ=255)
                                                                                                                                           ccGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTTGCTGCa      >  1:1234055/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)
                                                                                                                                               cAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTTGCTGCAATGt  >  1:870727/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTTGCTGCAATGT  >  NC_000913/2218367‑2218659

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: