Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 385,724 T→A 5.9% W165R (TGG→AGG)  tauA → taurine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913385,7240TA 5.9% 758.8 / 26.0 216W165R (TGG→AGG) tauAtaurine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (93/110);  new base A (6/7);  total (99/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

Only 100 of 218 total aligned reads displayed.

CAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGC  >  NC_000913/385576‑385872
                                                                                                                                                    |                                                                                                                                                    
        aaGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTCCAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCAGGCAGCga                                                                                                                                              <  1:936617/149‑2 (MQ=255)
         aGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCAGGCAGCgat                                                                                                                                             <  2:1338893/149‑3 (MQ=255)
          gATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCAGGCAGCgata                                                                                                                                            <  1:144731/149‑4 (MQ=255)
           aTCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGga                                                                                                                                           <  4:1018579/149‑1 (MQ=255)
           aTCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCAGGCAGCgataa                                                                                                                                           <  1:972300/149‑5 (MQ=255)
            tCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGgag                                                                                                                                          >  2:571243/1‑149 (MQ=255)
              tGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGat                                                                                                                                        >  3:1136529/1‑149 (MQ=255)
              tGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGat                                                                                                                                        >  3:483696/1‑149 (MQ=255)
               gATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGata                                                                                                                                       >  4:455437/1‑149 (MQ=255)
                  tGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGata                                                                                                                                       <  2:718890/146‑1 (MQ=255)
                   ggCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGa                                                                                                                                   <  1:1243640/149‑1 (MQ=255)
                   ggCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGa                                                                                                                                   <  1:495741/149‑1 (MQ=255)
                    gCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGAt                                                                                                                                  >  1:1288324/1‑149 (MQ=255)
                     cAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATg                                                                                                                                 <  3:303051/149‑1 (MQ=255)
                     cAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATg                                                                                                                                 <  1:1433995/149‑1 (MQ=255)
                      aaaCGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATgg                                                                                                                                >  2:45354/1‑149 (MQ=255)
                      aaaCGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATgg                                                                                                                                >  4:1354807/1‑149 (MQ=255)
                       aaCGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGt                                                                                                                               <  1:626291/149‑1 (MQ=255)
                        aCGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTg                                                                                                                              <  3:65152/149‑1 (MQ=255)
                         cGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGc                                                                                                                             <  2:305210/149‑1 (MQ=255)
                               gCCGTACCGTTTATCTCCACCACCCACTAAAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGt                                                                                                                       >  3:1185435/1‑149 (MQ=255)
                                  gTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTg                                                                                                                    <  4:373254/149‑1 (MQ=255)
                                   tACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTgg                                                                                                                   >  4:808992/1‑149 (MQ=255)
                                      cGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGAt                                                                                                                                  <  4:756430/131‑1 (MQ=255)
                                      cGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGAt                                                                                                                                  >  3:756430/1‑131 (MQ=255)
                                       gTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCAc                                                                                                               >  3:222995/1‑149 (MQ=255)
                                       gTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCCGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCAc                                                                                                               >  2:398467/1‑149 (MQ=255)
                                        tttATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCAcc                                                                                                              >  3:536187/1‑149 (MQ=255)
                                         ttATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACcg                                                                                                             >  1:101380/1‑149 (MQ=255)
                                           aTCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACcggc                                                                                                           >  4:1142281/1‑149 (MQ=255)
                                            tctcCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCa                                                                                                                                                  <  2:872941/108‑1 (MQ=255)
                                              tcCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGt                                                                                                        >  4:1411843/1‑149 (MQ=255)
                                              tcCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGt                                                                                                        <  3:1436523/149‑1 (MQ=255)
                                              tcCACCACCCACCACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCATCCGCCCGCGGTTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATCTCTAGGCACCGGCGGt                                                                                                        >  4:73781/1‑149 (MQ=255)
                                                 accaccCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTaa                                                                                                     <  3:1411843/149‑1 (MQ=255)
                                                  ccacccacTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAAc                                                                                                    <  3:558201/149‑1 (MQ=255)
                                                    acccacTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGc                                                                                                  >  1:205421/1‑149 (MQ=255)
                                                    acccacTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGc                                                                                                  >  1:382313/1‑149 (MQ=255)
                                                     cccacTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGcc                                                                                                 <  1:651568/149‑1 (MQ=255)
                                                       cacTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCt                                                                                               <  2:1150712/149‑1 (MQ=255)
                                                          tACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGa                                                                                            >  3:891129/1‑149 (MQ=255)
                                                             aGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGaaaa                                                                                         <  3:360639/149‑1 (MQ=255)
                                                              gCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGaaaaa                                                                                        <  1:45354/149‑1 (MQ=255)
                                                                 tgctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGAc                                                                                     <  3:455437/149‑1 (MQ=255)
                                                                  gctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACg                                                                                    <  3:14281/149‑1 (MQ=255)
                                                                  gctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACg                                                                                    <  2:664523/149‑1 (MQ=255)
                                                                   ctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACgg                                                                                   <  4:604065/149‑1 (MQ=255)
                                                                   ctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACgg                                                                                   <  3:376318/149‑1 (MQ=255)
                                                                       cggcACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAg                                                                               <  1:322887/149‑1 (MQ=255)
                                                                        ggcACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCTCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAgg                                                                              >  2:1297105/1‑149 (MQ=255)
                                                                         gcACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGgt                                                                             <  4:1196993/149‑1 (MQ=255)
                                                                          cACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGgtg                                                                            <  2:205421/149‑1 (MQ=255)
                                                                           aCTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTg                                                                                                                    >  2:510675/1‑108 (MQ=255)
                                                                           aCTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGgtgt                                                                           >  3:1416205/1‑149 (MQ=255)
                                                                           aCTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGgtgt                                                                           >  4:43767/1‑149 (MQ=255)
                                                                           aCTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGgtgt                                                                           >  1:1149132/1‑149 (MQ=255)
                                                                            cTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGtt                                                                          <  4:222995/149‑1 (MQ=255)
                                                                               aaaCACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                       <  1:686197/149‑1 (MQ=255)
                                                                                    cTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGAtt                                                                  <  2:441595/149‑1 (MQ=255)
                                                                                    cTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGAtt                                                                  <  1:745486/149‑1 (MQ=255)
                                                                                     tGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTc                                                                 >  2:786914/1‑149 (MQ=255)
                                                                                      ggggCATTAAACCCGGCCAAGTGGAGATTGTGAACCTGCAGCCGCCAGCGATTATAGCTGCCTGACAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCt                                                                <  4:870322/149‑1 (MQ=255)
                                                                                       gggCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGgt                                                                             >  1:507358/1‑135 (MQ=255)
                                                                                         gCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGaa                                                             <  3:164974/149‑1 (MQ=255)
                                                                                            ttAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                       <  2:159628/136‑1 (MQ=255)
                                                                                            ttAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                       >  1:159628/1‑136 (MQ=255)
                                                                                                  ccGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTAt                                                                                                                         >  4:169550/1‑80 (MQ=255)
                                                                                                  ccGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCggg                                                    >  3:73657/1‑149 (MQ=255)
                                                                                                      gCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGt                                                <  3:238682/149‑1 (MQ=255)
                                                                                                          gTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTgggg                                            >  2:783574/1‑149 (MQ=255)
                                                                                                           tGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                       >  2:169699/1‑121 (MQ=255)
                                                                                                               gaTTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGgcgcgc                                       <  2:978791/149‑1 (MQ=255)
                                                                                                               gaTTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGgcgcgc                                       >  3:95513/1‑149 (MQ=255)
                                                                                                               gaTTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGACCGCGGGGGGATATTGATGGTGCTTATGTATGGGCACCGGCGGATAACGCCCTGGAAAAAGACGGCCAGGTGTTGGCCGAATCTGAACAGGTCGGGCAGTGGGgcgcgc                                       >  2:72591/1‑149 (MQ=255)
                                                                                                                     gAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGc                                 >  2:802382/1‑149 (MQ=255)
                                                                                                                       aCCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTg                               <  1:293822/149‑1 (MQ=255)
                                                                                                                          tGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGAc                            <  3:808992/149‑1 (MQ=255)
                                                                                                                          tGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGAc                            <  3:1204878/149‑1 (MQ=255)
                                                                                                                                cgcccgcGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTgg                      <  3:1354807/149‑1 (MQ=255)
                                                                                                                                 gcccgcGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTggg                     <  4:536187/149‑1 (MQ=255)
                                                                                                                                  cccgcGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGggt                    >  2:1041542/1‑149 (MQ=255)
                                                                                                                                   ccgcGATCATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGggtg                   >  1:1083996/1‑149 (MQ=255)
                                                                                                                                    cgcgATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACgg                                                                                   >  1:1062893/1‑84 (MQ=255)
                                                                                                                                    cgcgATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGggtgg                  >  2:965995/1‑149 (MQ=255)
                                                                                                                                    cgcgATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGggtgg                  >  2:1145957/1‑149 (MQ=255)
                                                                                                                                       gATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTgc               <  3:56950/149‑1 (MQ=255)
                                                                                                                                        aTTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCaa                                                                                >  1:857549/1‑83 (MQ=255)
                                                                                                                                         ttATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTgcgc             >  4:1279753/1‑149 (MQ=255)
                                                                                                                                         tccccGCTGCCAGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTgcgc             >  4:1081443/5‑149 (MQ=255)
                                                                                                                                          ccccGCTGCCAGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCa            >  1:574552/4‑149 (MQ=255)
                                                                                                                                           cccGCTGCCAGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGgcgcgc                                       >  4:680507/3‑121 (MQ=255)
                                                                                                                                           cccGCTGCCAGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGgcgcgc                                       <  3:680507/119‑1 (MQ=255)
                                                                                                                                           cccGCTGCCAGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGggtggt                 >  4:1372757/3‑143 (MQ=255)
                                                                                                                                           cccGCTGCCAGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCaa           >  2:1406804/3‑149 (MQ=255)
                                                                                                                                            tCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGAtt                                                                  >  2:1282857/1‑93 (MQ=255)
                                                                                                                                            tCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCaaa          <  1:802382/149‑1 (MQ=255)
                                                                                                                                            tCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCaaa          <  1:254546/149‑1 (MQ=255)
                                                                                                                                               cTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGAt       <  1:786914/149‑1 (MQ=255)
                                                                                                                                                   cTGGCATCGGGGAGATATTGATGTTGCTTATGTCTGGGCACCGGCGGTTAACGCCTTGGAAAAAGACGGCAAGGTGGTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTTGTGCGCAAAGATTTTg   <  2:1112875/149‑1 (MQ=255)
                                                                                                                                                   cTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTg   <  4:1420199/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGC  >  NC_000913/385576‑385872

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: