Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 854,796 G→A 7.5% intergenic (+31/+28) ybiR → / ← ldtB putative transporter/L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913854,7960GA 7.5% 371.3 / 12.1 120intergenic (+31/+28)ybiR/ldtBputative transporter/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base G (63/48);  new base A (7/2);  total (70/50)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

Only 100 of 122 total aligned reads displayed.

GCTACGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTT  >  NC_000913/854648‑854945
                                                                                                                                                    |                                                                                                                                                     
gcTAAGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCGCTATTCAATACAGATGCTGTTGTGGGAGGCGTTGGTGGGATATTTTTTGTTAGTTATACTCCAGCCCAACTATTGTCTGTCAAATAAAAAAGGCGGATTATGAGt                                                                                                                                                      <  2:1671185/150‑1 (MQ=255)
   aCGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCg                                                                                                                                                   <  2:549842/150‑1 (MQ=255)
   aCGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCg                                                                                                                                                   >  1:1247407/1‑150 (MQ=255)
     ttATGGCGACCGATCGCCGCAGCTGGTGGCGTTTCCATCTATATTCAAAACCGAAGCTGTTGAGGTCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCTCGGCCTACTAGGGTCTGGCAAAAAAAAAAGGCGGATTATGAGTCCGcc                                                                                                                                                 <  1:272518/149‑1 (MQ=255)
       aTGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCAtt                                                                                                                                               >  2:504838/1‑150 (MQ=255)
       aTGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGATATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCtt                                                                                                                                               >  2:949691/1‑150 (MQ=255)
        tGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCttt                                                                                                                                              <  1:1289018/150‑1 (MQ=255)
           cGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGc                                                                                                                                           <  2:786251/150‑1 (MQ=255)
              aCGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCttt                                                                                                                                        >  1:203603/1‑150 (MQ=255)
               cGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGACTCATAATCCGCCttttttattt                                                                                                                                       <  1:968178/150‑7 (MQ=255)
                 aTCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATAGCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAg                                                                                                                                     <  1:739752/150‑1 (MQ=255)
                   cgccgcATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGtt                                                                                                                                   <  1:1590576/150‑1 (MQ=255)
                   cgccgcATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGtt                                                                                                                                   >  2:271853/1‑150 (MQ=255)
                     ccgcATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTGGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGtttt                                                                                                                                 >  2:1165245/1‑150 (MQ=255)
                     ccgcATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGtttt                                                                                                                                 <  1:1028201/150‑1 (MQ=255)
                       gcATCTGGTGGCGTTTCCATCTCTATTCAATACCGCTGCTGGTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGTCGGATTATGAGTACGCCTTTTGCTTTAAGTTTTGa                                                                                                                               >  2:1682790/1‑150 (MQ=255)
                       gcATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGa                                                                                                                               >  1:39544/1‑150 (MQ=255)
                        gaTCTGGTGTCGTTTCCATATCTATTCAATACCGATTCTGTTGGGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGaa                                                                                                                              <  2:694187/149‑1 (MQ=255)
                         aTCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAg                                                                                                                             <  1:579500/150‑1 (MQ=255)
                                  cGTTTGCATCTGTTATCATTAAAGATGCTGTTTTGGGCGGCGTTGGTGGTTTATGTTTTGTTAGTTATACTCCCGGCCAACTTGTGTATGTCAAATAAAAAAGGCGGAGTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCa                                                                                                                    <  2:919502/150‑1 (MQ=255)
                                  cGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCa                                                                                                                    >  1:1705284/1‑150 (MQ=255)
                                  cGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTACCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTAAAGTTTTGAAGATTAATTCa                                                                                                                    >  2:744974/1‑150 (MQ=255)
                                  cGTTTCCATCTCAATTCAAGACCGAGGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCa                                                                                                                    <  2:231308/150‑1 (MQ=255)
                                      tCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACg                                                                                                                <  1:465856/150‑1 (MQ=255)
                                      tCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGACTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACg                                                                                                                <  1:619826/150‑1 (MQ=255)
                                          tctcTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAAcc                                                                                                            >  2:771579/1‑150 (MQ=255)
                                              tATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCa                                                                                                        <  1:176467/150‑1 (MQ=255)
                                                  cAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATccc                                                                                                    >  1:410424/1‑150 (MQ=255)
                                                      aCCGCTGCTGGTGGTGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAg                                                                                                <  2:311414/150‑1 (MQ=255)
                                                        cGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCg                                                                                              >  1:858512/1‑150 (MQ=255)
                                                        cGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCg                                                                                              <  1:966996/150‑1 (MQ=255)
                                                         gATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCgg                                                                                             >  2:1642845/1‑150 (MQ=255)
                                                         gATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCgg                                                                                             >  1:1096594/1‑150 (MQ=255)
                                                             cTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGtttt                                                                                         >  1:1661863/1‑150 (MQ=255)
                                                              tgttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGttttt                                                                                        <  2:1247407/150‑1 (MQ=255)
                                                               gttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGATGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTa                                                                                       <  1:747766/150‑1 (MQ=255)
                                                                ttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTaa                                                                                      <  1:1367238/150‑1 (MQ=255)
                                                                 tgttGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTTGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATGCAGACGAACCGGCATCCCGGAGCGGTTTTTAAt                                                                                     <  1:740634/150‑1 (MQ=255)
                                                                 tgtgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAAt                                                                                     <  1:1272400/150‑1 (MQ=255)
                                                                  gtgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATc                                                                                    <  1:1618031/150‑1 (MQ=255)
                                                                   tgtgcggtgTTGGGGGGATATGGTTTGGTAGTTATACTACCGTCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTTCGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCTGAGCGGTTTTTACTCg                                                                                   <  2:285709/147‑1 (MQ=255)
                                                                   tgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCg                                                                                   >  1:503709/1‑150 (MQ=255)
                                                                     ggcggcGTTGGTGGGATATGTTTTGTTGGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCt                                                                                 <  1:1165245/150‑1 (MQ=255)
                                                                     ggcggcGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGCTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCt                                                                                 <  1:300735/150‑1 (MQ=255)
                                                                          cgTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTtcatc                                                                            <  1:498453/150‑1 (MQ=255)
                                                                           gTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGACCCGGCATCCCGGAGCGGTTTTTAATCGCTtcatca                                                                           >  2:1546826/1‑150 (MQ=255)
                                                                           gTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTtcatca                                                                           >  2:1540003/1‑150 (MQ=255)
                                                                           gTTGGTGGGATATGTATTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTtcatca                                                                           <  2:302803/150‑1 (MQ=255)
                                                                                 gggATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGaaca                                                                     >  2:926905/1‑150 (MQ=255)
                                                                                  ggATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGaacaa                                                                    >  1:956119/1‑150 (MQ=255)
                                                                                   gATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCGTTTTGGTTTAAGTTTTGATGATTAATTAAGACGAACCGGCATCCCGTAGCGGTTTTTAATCGGTTAATCAAGaacaac                                                                   >  1:392987/1‑150 (MQ=255)
                                                                                     tatGTTTTGTTATTTATAATCCAGTCCAACTAGGGTCTGGAAAATAAAAAATGAGCATTATGAGTCCGCCTTTTGCTTTAAGTTTTTAAGATTAATTCAGACTAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCt                                                                 <  1:924735/150‑1 (MQ=255)
                                                                                     tatGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCt                                                                 >  2:568089/1‑150 (MQ=255)
                                                                                        gTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGt                                                              >  2:1268260/1‑150 (MQ=255)
                                                                                          tttGTTAGTTATACTCCCGGCCAACTATGGTCTGGCAAATAAAAAAGTCGGATTATGAGTCAGCCTTTTGCTTTAAGTTTTGAATATTAATTCAGACGAACCGGCATCCCGGATCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCa                                                            <  1:405550/150‑1 (MQ=255)
                                                                                           ttGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCaa                                                           >  2:371565/1‑150 (MQ=255)
                                                                                             gttagttaTACTCCCGGCGAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACa                                                         >  1:587266/1‑150 (MQ=255)
                                                                                                agttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCt                                                      <  1:803714/150‑1 (MQ=255)
                                                                                                   tataCTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGc                                                   >  1:1180948/1‑150 (MQ=255)
                                                                                                   tataCTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGc                                                   <  2:1556907/150‑1 (MQ=255)
                                                                                                     taCTCCCGTCCATCTAGGGTCTGGCAAAAAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGAGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTg                                                 <  1:949691/150‑1 (MQ=255)
                                                                                                       cTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGAc                                               >  1:554657/1‑150 (MQ=255)
                                                                                                       cTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGAc                                               >  1:561448/1‑150 (MQ=255)
                                                                                                            ggCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTc                                          >  2:1696688/1‑150 (MQ=255)
                                                                                                             gCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCa                                         >  1:1337130/1‑150 (MQ=255)
                                                                                                              ccAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGATTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCAc                                        <  1:1642845/150‑1 (MQ=255)
                                                                                                               cAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCATCATCTGGCTGACCGGTCACt                                       >  2:696927/1‑150 (MQ=255)
                                                                                                               cAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACAGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCAcg                                       >  2:670505/1‑149 (MQ=255)
                                                                                                                   tAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCt                                   <  2:373198/150‑1 (MQ=255)
                                                                                                                    aGGGTCTGTCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTg                                  <  1:271853/150‑1 (MQ=255)
                                                                                                                      ggTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCGTTTGCTTTATGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGATTCATCAAGAACAACCTGGTAAACATCTGGCTGACCGGTCTCTGTCTGCa                                >  1:150510/1‑150 (MQ=255)
                                                                                                                         cTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGc                             >  2:1395167/1‑150 (MQ=255)
                                                                                                                         cTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGc                             <  2:1705284/150‑1 (MQ=255)
                                                                                                                           ggCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGTTCAACATCTGGCTGAAATGTCACTGTCTGCACGctc                           >  2:390914/1‑150 (MQ=255)
                                                                                                                           ggCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGctc                           >  2:1632710/1‑150 (MQ=255)
                                                                                                                           ggCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGctc                           >  1:1709813/1‑150 (MQ=255)
                                                                                                                            gCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGctct                          >  2:303146/1‑150 (MQ=255)
                                                                                                                             cAAATATAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCtt                         <  2:785635/150‑1 (MQ=255)
                                                                                                                               aaagcaaaaGGCGGACTCATAATCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCg                       >  1:478825/6‑150 (MQ=255)
                                                                                                                               aaTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCg                       <  1:1482291/150‑1 (MQ=255)
                                                                                                                                aagcaaaaGGCGGACTCATAATCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGt                      >  2:494008/5‑150 (MQ=255)
                                                                                                                                aTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGt                      >  1:512886/1‑150 (MQ=255)
                                                                                                                                 tAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTc                     >  2:1281877/1‑150 (MQ=255)
                                                                                                                                  gcaaaaGGCGGACTCATAATCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCa                    >  2:730217/3‑150 (MQ=255)
                                                                                                                                   caaaaGGCGGACTCATAATCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAg                   >  1:1576302/2‑150 (MQ=255)
                                                                                                                                   caaaaGGCGGACTCATAATCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAg                   >  1:990869/2‑150 (MQ=255)
                                                                                                                                    aaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAgg                  >  1:890456/1‑150 (MQ=255)
                                                                                                                                    aaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCACCAAGAACAACCTCGTCAACATCTGGCTGACCGCTCACTGTCTGCACGCTCTTCGTCAgg                  >  1:329208/1‑150 (MQ=255)
                                                                                                                                      aaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACATGGTCAACATCTGGCTGAACGGTCACTGTCTGCACGCTCTTCGTCAGGGt                >  1:1059940/1‑150 (MQ=255)
                                                                                                                                        ggCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTaa              >  1:334960/1‑150 (MQ=255)
                                                                                                                                           ggATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTg           >  1:474244/1‑150 (MQ=255)
                                                                                                                                           ggATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTg           >  1:467724/1‑150 (MQ=255)
                                                                                                                                            gATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTgg          >  1:911888/1‑150 (MQ=255)
                                                                                                                                               tATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATAAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTTTCTGCACGCTCTTCGTCAGGGTAATTGGcac       >  1:767808/1‑150 (MQ=255)
                                                                                                                                                   aGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGCCCGGGCCCTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAtt   >  1:1675877/1‑150 (MQ=255)
                                                                                                                                                   aGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAtt   <  2:561448/150‑1 (MQ=255)
                                                                                                                                                   aGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAtt   <  2:554657/150‑1 (MQ=255)
                                                                                                                                                   aGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAtt   >  2:824273/1‑150 (MQ=255)
                                                                                                                                                   aGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAtt   >  1:1669513/1‑150 (MQ=255)
                                                                                                                                                    gTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAttt  <  2:1180948/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTACGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTT  >  NC_000913/854648‑854945

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: