Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 854,794 G→T 7.5% intergenic (+29/+30) 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,7940GT 7.5% 342.2 / 12.3 114intergenic (+29/+30)ybiR/ldtBputative transporter/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base G (60/45);  new base T (7/2);  total (67/47)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 115 total aligned reads displayed.

GCTACGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAT  >  NC_000913/854648‑854943
                                                                                                                                                  |                                                                                                                                                     
gcTAAGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCGCTATTCAATACAGATGCTGTTGTGGGAGGCGTTGGTGGGATATTTTTTGTTAGTTATACTCCAGCCCAACTATTGTCTGTCAAATAAAAAAGGCGGATTATGAGt                                                                                                                                                    <  2:1671185/150‑1 (MQ=255)
   aCGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCg                                                                                                                                                 >  1:1247407/1‑150 (MQ=255)
   aCGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCg                                                                                                                                                 <  2:549842/150‑1 (MQ=255)
    cGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGc                                                                                                                                                <  1:963649/150‑1 (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)
           cGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATAACGATGCTGTTGTGGGCTGCTTTGGTGGGATATGTATTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGACTCATAATCCGCCtttttt                                                                                                                                         <  1:219906/150‑3 (MQ=255)
             aaCGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGACTTTTGCtt                                                                                                                                       >  1:442272/1‑150 (MQ=255)
              aCGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCttt                                                                                                                                      >  2:803714/1‑150 (MQ=255)
              aCGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCttt                                                                                                                                      >  1:203603/1‑150 (MQ=255)
               cGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGACTCATAATCCGCCttttttattt                                                                                                                                     <  1:968178/150‑7 (MQ=255)
                 aTCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATAGCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAg                                                                                                                                   <  1:739752/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)
                      cgcATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGttttt                                                                                                                              >  1:373198/1‑149 (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)
                                     ttCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGAc                                                                                                               <  1:1617916/150‑1 (MQ=255)
                                      tCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACg                                                                                                              <  1:465856/150‑1 (MQ=255)
                                      tCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGACTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACg                                                                                                              <  1:619826/150‑1 (MQ=255)
                                              tATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCa                                                                                                      <  1:176467/150‑1 (MQ=255)
                                                  cAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATccc                                                                                                  >  1:410424/1‑150 (MQ=255)
                                                        cGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCg                                                                                            <  1:966996/150‑1 (MQ=255)
                                                        cGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCg                                                                                            >  1:858512/1‑150 (MQ=255)
                                                         gATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCgg                                                                                           >  2:1642845/1‑150 (MQ=255)
                                                         gATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCgg                                                                                           >  1:1096594/1‑150 (MQ=255)
                                                          aTGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGt                                                                                          >  2:48757/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)
                                                                  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)
                                                                           gTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTtcatca                                                                         >  2:1482291/1‑150 (MQ=255)
                                                                           gTTGGTGGGATATGTATTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTtcatca                                                                         <  2:302803/150‑1 (MQ=255)
                                                                                tggGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGaac                                                                    >  1:501168/1‑150 (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)
                                                                                               tagttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCagccag                                                     >  1:1368272/1‑145 (MQ=255)
                                                                                                agttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCt                                                    <  1:803714/150‑1 (MQ=255)
                                                                                                   tataCTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGc                                                 >  1:1180948/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                         >  2:307729/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)
                                                                                                                                 agcaaaaGGCGGACTCATAATCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGAGGATTTTAATCGATTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTc                   >  2:711924/4‑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)
                                                                                                                                     aaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAggg               >  1:1222902/1‑150 (MQ=255)
                                                                                                                                     aaaGGCGGACTCATAATCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAggg               >  2:1647958/1‑150 (MQ=255)
                                                                                                                                      aaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACATGGTCAACATCTGGCTGAACGGTCACTGTCTGCACGCTCTTCGTCAGGGt              >  1:1059940/1‑150 (MQ=255)
                                                                                                                                           ggATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTg         >  1:474244/1‑150 (MQ=255)
                                                                                                                                            gATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTgg        >  1:911888/1‑150 (MQ=255)
                                                                                                                                               tATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGcac     <  1:1696688/150‑1 (MQ=255)
                                                                                                                                               tATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATAAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTTTCTGCACGCTCTTCGTCAGGGTAATTGGcac     >  1:767808/1‑150 (MQ=255)
                                                                                                                                                  gAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAt  <  1:568089/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCTACGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAT  >  NC_000913/854648‑854943

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: