Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 854,793 T→A 8.3% intergenic (+28/+31) 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,7930TA 8.3% 319.4 / 14.4 114intergenic (+28/+31)ybiR/ldtBputative transporter/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base T (60/44);  new base A (7/2);  total (67/47)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 114 total aligned reads displayed.

GCTACGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCAC  >  NC_000913/854648‑854940
                                                                                                                                                 |                                                                                                                                                   
gcTAAGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCGCTATTCAATACAGATGCTGTTGTGGGAGGCGTTGGTGGGATATTTTTTGTTAGTTATACTCCAGCCCAACTATTGTCTGTCAAATAAAAAAGGCGGATTATGAGt                                                                                                                                                 <  2:1671185/150‑1 (MQ=255)
   aCGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCg                                                                                                                                              >  1:1247407/1‑150 (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)
           cGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGc                                                                                                                                      <  2:786251/150‑1 (MQ=255)
           cGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATAACGATGCTGTTGTGGGCTGCTTTGGTGGGATATGTATTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGACTCATAATCCGCCtttttt                                                                                                                                      <  1:219906/150‑3 (MQ=255)
             aaCGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGACTTTTGCtt                                                                                                                                    >  1:442272/1‑150 (MQ=255)
              aCGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCttt                                                                                                                                   <  2:567571/150‑1 (MQ=255)
              aCGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCttt                                                                                                                                   >  2:803714/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)
                   cgccgcATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGtt                                                                                                                              <  1:1590576/150‑1 (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)
                                  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:966996/150‑1 (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)
                                                                 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: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)
                                                                                     tatGTTTTGTTATTTATAATCCAGTCCAACTAGGGTCTGGAAAATAAAAAATGAGCATTATGAGTCCGCCTTTTGCTTTAAGTTTTTAAGATTAATTCAGACTAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCt                                                            <  1:924735/150‑1 (MQ=255)
                                                                                     tatGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCt                                                            >  2:568089/1‑150 (MQ=255)
                                                                                       tGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTgg                                                          >  2:1192247/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)
                                                                                               tagttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCagccag                                                  >  1:1368272/1‑145 (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                                     <  1:745047/150‑1 (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)
                                                                                                                    aGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTg                             <  2:203603/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                      >  1:1709813/1‑150 (MQ=255)
                                                                                                                           ggCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGctc                      >  2:1632710/1‑150 (MQ=255)
                                                                                                                            gCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGctct                     >  2:303146/1‑150 (MQ=255)
                                                                                                                             cAAATATAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCtt                    <  2:785635/150‑1 (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)
                                                                                                                                    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)
                                                                                                                                        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)
                                                                                                                                               tATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGcac  <  1:1696688/150‑1 (MQ=255)
                                                                                                                                               tATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATAAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTTTCTGCACGCTCTTCGTCAGGGTAATTGGcac  >  1:767808/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCTACGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCAC  >  NC_000913/854648‑854940

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: