Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,674,403 T→C 6.2% intergenic (+28/‑383) yhjD → / → yhjE inner membrane putative BrbK family alternate lipid exporter/putative MFS transporter; membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,674,4030TC 6.2% 375.8 / 11.5 129intergenic (+28/‑383)yhjD/yhjEinner membrane putative BrbK family alternate lipid exporter/putative MFS transporter; membrane protein
Reads supporting (aligned to +/- strand):  ref base T (62/59);  new base C (5/3);  total (67/62)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

Only 100 of 134 total aligned reads displayed.

CCGTGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTCATTCC  >  NC_000913/3674254‑3674552
                                                                                                                                                     |                                                                                                                                                     
ccGTGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTGTGACGCTGTTTTGCGCGGCGGGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCtt                                                                                                                                                       <  2:664594/150‑1 (MQ=255)
    gCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                                   >  2:51016/1‑150 (MQ=255)
      tGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCGGTTTTGCGCGGCGTGTATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                                 <  2:1457277/150‑1 (MQ=255)
        gggCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATAAGGCGTTTACGCCGCATCCGGcc                                                                                                                                               <  1:1676193/150‑2 (MQ=255)
        gggCGGTTGGCAGTATCCTACTGCTTCGCCCGTGGGACGCTGTTTTGCGCGGCGGGGATTGCCACCTCCGAATATAAAGACGACCCGGGAATGCCGGGGAAATCGGAGCCTTAAAATAGGCCGGATAAGGCGTTTACGCCGCATCCGGcc                                                                                                                                               <  2:243146/150‑2 (MQ=255)
         ggCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCAt                                                                                                                                              <  1:1233788/150‑1 (MQ=255)
          gCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTTGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATa                                                                                                                                             <  1:1166247/150‑1 (MQ=255)
             gATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGa                                                                                                                                          >  2:1330706/1‑150 (MQ=255)
             gATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGa                                                                                                                                          >  1:280224/1‑150 (MQ=255)
                 gCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGAtttt                                                                                                                                      >  1:933315/1‑150 (MQ=255)
                 gCATTTTTCTACTTCTTCGCCCGTTTGACGCTGCTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGAtttt                                                                                                                                      <  1:1142913/150‑1 (MQ=255)
                   ctctttCTACTTCTTCGCCCGTCTGACGCCGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCTGCGTAAACGCCTTATCCGGCATACGATTTTAt                                                                                                                                    <  1:1144389/147‑1 (MQ=255)
                     ttttCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTc                                                                                                                                  <  1:1330706/150‑1 (MQ=255)
                      tttCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCg                                                                                                                                 <  1:322661/150‑1 (MQ=255)
                           aCTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTg                                                                                                                            >  1:594462/1‑150 (MQ=255)
                           aCTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTg                                                                                                                            >  1:42114/1‑150 (MQ=255)
                            cttcttCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTgg                                                                                                                           >  2:1150736/1‑150 (MQ=255)
                             ttcttcGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGc                                                                                                                          >  1:1424683/1‑150 (MQ=255)
                              tcttcGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAACCGCAGCCTTAAACTCGGCCGGATGCGGCGTAAACGCCGTATCCGGCATACGATTTTATTCGGCCTGGCt                                                                                                                         <  1:291776/150‑1 (MQ=255)
                              tcttcGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCt                                                                                                                         >  1:255744/1‑150 (MQ=255)
                               cttcGCCCGTAAGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGAGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTc                                                                                                                        >  2:937337/1‑150 (MQ=255)
                                ttcGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAACGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTcc                                                                                                                       <  1:1254000/150‑1 (MQ=255)
                                 tcGCCCGTTTGACGCTGGTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTccc                                                                                                                      <  2:306540/150‑1 (MQ=255)
                                 tcGCCCGTTTGACGCGGTGTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTccc                                                                                                                      <  1:174774/150‑1 (MQ=255)
                                  cGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTcccc                                                                                                                     <  2:1421708/150‑1 (MQ=255)
                                   gCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCg                                                                                                                    >  2:1701540/1‑150 (MQ=255)
                                    cccGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGt                                                                                                                   >  2:773236/1‑150 (MQ=255)
                                      cGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAg                                                                                                                 >  1:945364/1‑150 (MQ=255)
                                        tttGCCGCGGTTTTGCGCGGCGTGTATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGTGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGc                                                                                                               <  2:1061015/150‑1 (MQ=255)
                                              gCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATa                                                                                                         >  1:247512/1‑150 (MQ=255)
                                                 gTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGTAAAACGCACCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGa                                                                                                      >  2:1624888/1‑150 (MQ=255)
                                                 gTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGa                                                                                                      <  2:748875/150‑1 (MQ=255)
                                                 gTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGa                                                                                                      <  1:763463/150‑1 (MQ=255)
                                                    ttGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATgc                                                                                                   <  1:1196725/150‑1 (MQ=255)
                                                       cgcgGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGcc                                                                                                <  2:1233859/150‑1 (MQ=255)
                                                            cgTGGATTGCCACCGCCGAATATAAAGACTACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCAGGATAAGATGCGCCAgcat                                                                                           >  1:1257677/1‑150 (MQ=255)
                                                              gggATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATCGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAgcatcg                                                                                         <  2:1630633/149‑1 (MQ=255)
                                                              gggATTGCCACCGCCGAATATAAAGACGACCCGCGAAAGCCGGGGAAAACGCAGCCTAAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAgcatcg                                                                                         <  1:957088/149‑1 (MQ=255)
                                                               ggATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAgcatcgc                                                                                        <  2:1651775/150‑1 (MQ=255)
                                                                 aTTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCAGGCTCCCCGTAGGCCGGATAAGATGCGCCAgcatcgcat                                                                                      >  1:145000/1‑150 (MQ=255)
                                                                   tGCCACCGCCGAATATACAGACGCCCAGCGCATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATcc                                                                                    <  2:1671828/150‑1 (MQ=255)
                                                                    gCCCCCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATACg                                                                                   >  1:512482/1‑150 (MQ=255)
                                                                     ccACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCgg                                                                                  >  2:628479/1‑150 (MQ=255)
                                                                     ccACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCgg                                                                                  >  2:1680497/1‑150 (MQ=255)
                                                                      cACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGc                                                                                 >  1:207250/1‑150 (MQ=255)
                                                                       aCCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATGCGATGCTGGCGCATCTTATCCGGcc                                                                                <  1:75384/150‑2 (MQ=38)
                                                                         cgccgAATAAACAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCtat                                                                              <  2:274765/150‑1 (MQ=255)
                                                                          gccgAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAATTAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCTGCATCGCATCCTGCtata                                                                             <  2:1109704/150‑1 (MQ=255)
                                                                                 ataAAGACGCCCCGCGAATGCCGGGGAAAACCCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGca                                                                      <  1:111029/150‑1 (MQ=255)
                                                                                 ataAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGca                                                                      <  1:666785/150‑1 (MQ=255)
                                                                                  taAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGcac                                                                     >  2:670368/1‑150 (MQ=255)
                                                                                   aaaGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATTCGcaca                                                                    >  2:305094/1‑150 (MQ=255)
                                                                                      gacgacCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATaa                                                                 >  1:363636/1‑150 (MQ=255)
                                                                                        cgacCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAAcc                                                               <  1:1414949/150‑1 (MQ=255)
                                                                                           cccGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACctct                                                            >  2:717798/1‑150 (MQ=255)
                                                                                               cgAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGaa                                                        <  2:66789/150‑1 (MQ=255)
                                                                                                  aTGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACt                                                     <  2:452753/150‑1 (MQ=255)
                                                                                                  aTGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACt                                                     <  1:137346/150‑1 (MQ=255)
                                                                                                  aTGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACt                                                     <  2:1538759/150‑1 (MQ=255)
                                                                                                    gCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCg                                                   <  2:1693991/150‑1 (MQ=255)
                                                                                                    gCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCg                                                   <  1:1518840/150‑1 (MQ=255)
                                                                                                    gCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCg                                                   <  1:1610627/150‑1 (MQ=255)
                                                                                                     acGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCCCATAACCTCTTGAAACTCGt                                                  >  1:1316511/2‑150 (MQ=255)
                                                                                                      cGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTc                                                 >  1:1508793/1‑150 (MQ=255)
                                                                                                       ggggAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTcc                                                >  1:1049720/1‑150 (MQ=255)
                                                                                                       ggggAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTcc                                                >  2:1170620/1‑150 (MQ=255)
                                                                                                        gggAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTccc                                               <  2:1621452/150‑1 (MQ=255)
                                                                                                             aaCGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCagag                                          <  1:1595377/150‑1 (MQ=255)
                                                                                                              aCGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGc                                         >  2:25982/1‑150 (MQ=255)
                                                                                                               cGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGcc                                        <  1:1002963/150‑1 (MQ=255)
                                                                                                                    ccTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTc                                   >  2:1590504/1‑150 (MQ=255)
                                                                                                                    ccTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTc                                   >  2:1152068/1‑150 (MQ=255)
                                                                                                                     cTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCa                                  >  1:641956/1‑150 (MQ=255)
                                                                                                                     cTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCa                                  >  1:1149526/1‑150 (MQ=255)
                                                                                                                     cTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCa                                  >  1:457814/1‑150 (MQ=255)
                                                                                                                     cTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCGTcc                                  >  1:861369/1‑149 (MQ=255)
                                                                                                                         aaaTAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCa                              <  1:51016/150‑1 (MQ=255)
                                                                                                                             aGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATGTa                          >  2:1109854/1‑150 (MQ=255)
                                                                                                                              tgCCGGATAAGGCGTTTACGCCGCATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCatatat                         >  2:1160732/2‑150 (MQ=38)
                                                                                                                              tgCCGGATAAGGCGTTTACGCCGCATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTAt                         >  2:1545918/2‑150 (MQ=39)
                                                                                                                                  ggataaggcgTTTACGCCGCATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATAGCATCCGGCTATAAAGCGCACATAACCTCTTGAACCTCGTCCCCAGAGCCTCGTCAGCCATCTATTTgg                     >  1:555289/7‑150 (MQ=38)
                                                                                                                                   gATGCGGCGTAAACGCCTTCTCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTggg                    <  2:880049/150‑1 (MQ=255)
                                                                                                                                   gATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGCGCCTCTTCAGCCATCTATTTggg                    >  2:523893/1‑150 (MQ=255)
                                                                                                                                    aTGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGa                   <  1:506141/150‑1 (MQ=255)
                                                                                                                                        ggcgTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCaa               >  1:886700/1‑150 (MQ=255)
                                                                                                                                           gTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACa            >  1:24943/1‑150 (MQ=255)
                                                                                                                                             aaaCGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGCAACTCGTCCCCAGCGCCTCTTCAGCCATCTAGTTGGGAGCAAACAAt          >  2:1120424/1‑150 (MQ=255)
                                                                                                                                             aaaCGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAAt          >  2:1404492/1‑150 (MQ=255)
                                                                                                                                             aaaCGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAAt          <  2:533460/150‑1 (MQ=255)
                                                                                                                                             aaaCGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAAt          >  2:536102/1‑150 (MQ=255)
                                                                                                                                              aaCGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCACCATCGCATCAGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCGCTTCAGCCATCTATTTGGGAGCAAACAAtt         >  1:795288/1‑150 (MQ=255)
                                                                                                                                               aCGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCCGCGCCTCTTCAGCCATCTATTTGGGAGCAAACAAttt        >  2:1166566/1‑150 (MQ=255)
                                                                                                                                                cGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAAttta       >  2:1713657/1‑149 (MQ=255)
                                                                                                                                                 gCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGACCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTCa      <  2:28979/150‑1 (MQ=255)
                                                                                                                                                  ccTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTCAt     >  1:572911/1‑150 (MQ=255)
                                                                                                                                                  ccTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTCAt     <  1:149516/150‑1 (MQ=255)
                                                                                                                                                   cTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTCAtt    <  1:1310700/150‑1 (MQ=255)
                                                                                                                                                     tATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTCATTcc  <  1:670368/150‑1 (MQ=255)
                                                                                                                                                     tATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTCATTcc  <  2:1286745/150‑1 (MQ=255)
                                                                                                                                                     tATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTCATTcc  >  1:155227/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGTGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTCATTCC  >  NC_000913/3674254‑3674552

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