Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,169,382 C→G 6.4% intergenic (+52/+30) bhsA → / ← ldtC biofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,169,3820CG 6.4% 282.4 / 10.3 106intergenic (+52/+30)bhsA/ldtCbiofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base C (49/48);  new base G (5/2);  total (55/51)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

Only 100 of 106 total aligned reads displayed.

GCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGA  >  NC_000913/1169237‑1169525
                                                                                                                                                 |                                                                                                                                               
gcAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTATTTTATAAATAAGCATAAACCGTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCgt                                                                                                                                              >  2:700501/1‑149 (MQ=255)
gcAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTTACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGACgt                                                                                                                                              >  1:872055/1‑149 (MQ=255)
  aaaaaGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTACCCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATTAGGATTAACCCTCATTTAGGCCTGCTACTGCTGATTTTTTCCCCGCGACATGGCgtgt                                                                                                                                            >  1:914518/1‑149 (MQ=255)
  aaaaaGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCAGGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCGGCTTTTTTCCCCGCGACAGGCCgggt                                                                                                                                            >  2:859435/1‑149 (MQ=255)
    aaaGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTcg                                                                                                                                          >  1:1237700/1‑149 (MQ=255)
      aGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTcgcg                                                                                                                                        <  1:1348859/149‑1 (MQ=255)
        cGGATGAGATGGGCGCAAAAGCTTTCCGGATTACGTCTGTAACCGGTCCGACTACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCggg                                                                                                                                      <  2:1142698/149‑1 (MQ=255)
                aTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttttt                                                                                                                              <  1:1226718/149‑1 (MQ=255)
                  gggCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTa                                                                                                                            >  1:401033/1‑149 (MQ=255)
                  gggCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTa                                                                                                                            >  1:1180238/1‑149 (MQ=255)
                       cAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCCTTAATGCCTGCTACTACTGATTTTCTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCgg                                                                                                                       <  2:331027/149‑1 (MQ=255)
                        aaaaTCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCggg                                                                                                                      >  2:47790/1‑149 (MQ=255)
                            tCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTTTAAATAAGCATTAACCCTCAGTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCGGTGTCGCCGGGATTTTTTTATCCggtcgtt                                                                                                                  >  2:1516805/1‑144 (MQ=255)
                            tCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCttt                                                                                                                  <  1:614318/149‑1 (MQ=255)
                               ttCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACa                                                                                                               <  1:935263/149‑1 (MQ=255)
                                  cGTATTACTTCTGTACCCGGTCCGAATACCCTCCATGGAACAGCAGTAAGTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGATGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTGATACGGGCTGTACAGCg                                                                                                            >  2:1049493/1‑149 (MQ=255)
                                   tatttACTTTTGAAACCGTACCTAAAACCCTCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGtt                                                                                                          <  1:1031536/146‑1 (MQ=255)
                                     aTTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGttt                                                                                                         <  1:808889/149‑1 (MQ=255)
                                      ttACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTtg                                                                                                        >  2:586914/1‑149 (MQ=255)
                                      ttACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTtg                                                                                                        >  2:1440771/1‑149 (MQ=255)
                                       tACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTtgt                                                                                                       <  2:1091356/149‑1 (MQ=255)
                                        aCTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTtgtg                                                                                                      <  1:1449797/149‑1 (MQ=255)
                                            cTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCt                                                                                                  >  1:314546/1‑149 (MQ=255)
                                               tAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTACCCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATTCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCAc                                                                                               >  2:1089652/1‑149 (MQ=255)
                                               tAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCAc                                                                                               >  1:11277/1‑149 (MQ=255)
                                               tAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCAc                                                                                               <  2:1220083/149‑1 (MQ=255)
                                                  ccGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAGTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtg                                                                                            >  2:494924/1‑149 (MQ=255)
                                                   cGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgt                                                                                           >  2:202051/1‑148 (MQ=255)
                                                   cGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                                                           <  1:26133/149‑1 (MQ=255)
                                                   cGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATAAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGGCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                                                           <  2:1126659/149‑1 (MQ=255)
                                                     gTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCGTTACAGCGTTTGTGGGCGCACTtgatg                                                                                         >  1:698548/1‑149 (MQ=255)
                                                      tCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatga                                                                                        <  2:1237700/149‑1 (MQ=255)
                                                         gAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAcgg                                                                                     <  1:47790/149‑1 (MQ=255)
                                                           aTACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAcggcg                                                                                   <  1:575063/149‑1 (MQ=255)
                                                              cccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGGGTCGCGGGGATTTTTTTCTCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGAc                                                                                <  2:464800/149‑1 (MQ=255)
                                                               ccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGGGTTTGTGGGCTCACTTGATGCCGGCGGACa                                                                               >  2:670443/1‑149 (MQ=255)
                                                                   cATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTACCCCGCGACATGCGGGGTCGCGGGGTTTTTTTTTGCCGGGCTTAACGGCGTTTGTGGGCGCCATTGATGCCGGCGGCCAGccc                                                                           >  1:470096/1‑148 (MQ=255)
                                                                     tGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCAcc                                                                         >  1:1272800/1‑149 (MQ=255)
                                                                      ggAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCg                                                                        <  1:1440771/149‑1 (MQ=255)
                                                                      ggAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGGGTCGCGGGGATTTTTTTATCCGGGATTTACCGCGTTTGTGGGCTCGCTTGCTGACGGCGGACATCCACCg                                                                        >  1:845874/1‑149 (MQ=255)
                                                                       gAACAGCAGTAATTTATAACTAAGCATTCACCCTCATTAATGCCTGCTACTGCTGATTTATTCCCCGCGCCATGCCGTGTCGCGGGGATTTTTTTATCCGGGCGTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCgg                                                                       <  2:905476/149‑1 (MQ=255)
                                                                         ccagcagTAATTTATAAATAAGCATTCACCCTCATTACTGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGGCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCa                                                                     <  1:479382/148‑1 (MQ=255)
                                                                           agcagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACCGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTATACAGCGTTTGGGGGCTCCCTTGATGACGGCGGACATCCCCCGACagc                                                                   >  2:1357787/1‑147 (MQ=255)
                                                                             cagTAATTTCTAAATACGCATGTAGCCGCTTTAAACCCGGCTTCAGTTGATTGTTTCCGGGCGTAATTCATGTTCGAGGTGATTTTTTATTCAGGTCTTTACAGCGTTTGTGGGCTCTCTTGATGACGGCGGACTTCCACCGGCATccc                                                                 <  2:275509/149‑1 (MQ=255)
                                                                              agTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGTATTTTTTTATCCGGGCTTTAAAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCg                                                                >  1:377999/1‑149 (MQ=255)
                                                                              agTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCg                                                                >  1:503253/1‑149 (MQ=255)
                                                                                   tttATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACg                                                           <  1:1078192/149‑1 (MQ=255)
                                                                                     tataAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa                                                         <  2:335667/149‑1 (MQ=255)
                                                                                        aaataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGActt                                                      >  1:1225769/1‑147 (MQ=255)
                                                                                        aaataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcat                                                      <  1:859435/149‑1 (MQ=255)
                                                                                        aaataaGCATTAACCCTCAGTAATGCCTGCTACTGCTTATTTTTTCCCCGCGCCCTGCCGTGTCGCGGGGAGTGTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCCCATGAGGACGGCGGACAGCCACCGGCCTCCCGGAGCGGAcat                                                      >  2:1566333/1‑149 (MQ=255)
                                                                                         aataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatc                                                     <  2:5326/149‑1 (MQ=255)
                                                                                          ataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatcc                                                    >  2:497458/1‑149 (MQ=255)
                                                                                            aaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccat                                                  <  1:885302/149‑1 (MQ=255)
                                                                                             aGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACATCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccatc                                                 <  2:137285/149‑1 (MQ=255)
                                                                                             aGCATTAACCCTCATTAATGCCTGCTACTGCTGATGTTTTCCCCGCGACATGCCGGGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccatc                                                 <  2:1259463/149‑1 (MQ=255)
                                                                                              gCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATca                                                <  2:572913/149‑1 (MQ=255)
                                                                                               cATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATcac                                               >  1:695933/1‑149 (MQ=255)
                                                                                                 ttAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACAt                                             <  2:1515288/149‑1 (MQ=255)
                                                                                                    aCCCTCATTAATGCCTGCTACTGCTGATTTTATCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGtt                                          <  2:914518/149‑1 (MQ=255)
                                                                                                    aCCCTCATTAATGCCTGATACTGCTGATTTTTTCCCCTCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGtt                                          <  1:126774/149‑1 (MQ=255)
                                                                                                     cccTCAATACTGCCTGCTACTGCTGATTTTTTCCCCGCGACCTGCCGTGGCGCGGTGATTTTTTTATCCGGGCTTTACCGCGTTTGTGGGCTCACTTGATGACGGCGGACAGCCACCGGCATCCCGGAACGGACATCCATCACATGTTg                                         <  2:1329383/149‑1 (MQ=255)
                                                                                                      ccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATTCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTga                                        >  2:682673/1‑148 (MQ=255)
                                                                                                      ccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCTGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGc                                        >  2:515392/1‑149 (MQ=255)
                                                                                                      ccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGCCATCCCCCGGCATCCCGGAACGGAGATACATGACATGTTGc                                        >  1:1062963/1‑149 (MQ=255)
                                                                                                       cTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGATTTACAGCGTTTGTGGGCTCCCTTGATGACGGCGGACATCCACCGGCATCCCGGAACGCACAGCCGTCCCATGTTGCa                                       >  2:280657/1‑149 (MQ=255)
                                                                                                       cTCAATTATGCCTGCTACTTCTGATTTTTTAACAGGGAGATTCCGTGTAGAGGGGATTTTTTTATCCGGGATTTAATTCGTTTGTTGGCTCAATAGATGACTGCGGACATCCATCGTCATCCCGGAACGTACATCCATCACATGTTGCa                                       <  2:411225/149‑1 (MQ=255)
                                                                                                          aTTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCa                                    >  1:459670/1‑149 (MQ=255)
                                                                                                          aTTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCa                                    >  2:501784/1‑149 (MQ=255)
                                                                                                             aaTGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACtt                                 >  2:913184/1‑149 (MQ=255)
                                                                                                             aaTGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACtt                                 <  2:55442/149‑1 (MQ=255)
                                                                                                              aTGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGCACGGACATCCATCACATGTTGCATCACTTc                                >  1:1338683/1‑149 (MQ=255)
                                                                                                                     gtACTCCTGATTTTTTGCCCGCGACAGGCCGGGTCGCGGGGATTTTTTTAAGCGGGCTTTACAGCGTTTGTGGGCTGACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCTTCACATGTTGCATCACTTCAGCGTCa                         <  2:549313/148‑1 (MQ=255)
                                                                                                                     cTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCa                         <  2:1272800/149‑1 (MQ=255)
                                                                                                                     cTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCa                         <  1:134166/149‑1 (MQ=255)
                                                                                                                      tACTGCTGATTTTTTCCCCGCGACATTCCGTGGCGATGGGATTTTTTTATCCTGGGTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGCCATACACCGGCATCCCGGACCGGACATCCCTCACATGGTGCCGCCCTTCAGCGTCAg                        >  1:1049330/1‑149 (MQ=255)
                                                                                                                        ctgctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGtt                      >  1:1175110/1‑149 (MQ=255)
                                                                                                                         tgctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGttt                     <  2:1180238/149‑1 (MQ=255)
                                                                                                                          gctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtg                    >  2:331792/1‑149 (MQ=255)
                                                                                                                           ctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt                   <  2:902272/149‑1 (MQ=255)
                                                                                                                             gATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACAGCCATCACATGTTGCATGACTTCAGCGTCTGTTtgtgg                 >  1:506321/1‑148 (MQ=255)
                                                                                                                             gATTTTGTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCTGTTTTtgc                 <  2:178966/149‑1 (MQ=255)
                                                                                                                              aTTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACGTCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgct                >  1:847042/1‑149 (MQ=255)
                                                                                                                               ttttttCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgctg               <  1:700501/149‑1 (MQ=255)
                                                                                                                               aaaaatCCCCGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgctg               >  2:590438/6‑149 (MQ=255)
                                                                                                                               aaaaatCCCCGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGCCTTTACAGCGTTTGAGTGCTGACTTGATGAAGGCGGACATCCACCGGCATCCCGGAACGGACATCCAGCACATGTTGCATCCCTTCAGAGTCAGTTTGagctg               >  2:476904/6‑149 (MQ=255)
                                                                                                                                tttttCCCCGCTACATGACTTGTCGCGTGTATTTTTTTATTCGGGCTTTAGTGCTTTTGTGGGATCACTTGATGACGGCTTACATCCACCGGCATCCCGGAACGGACATCGATCACATTTTGCATCACTTCAGCGTCATTTTGtgctgc              <  1:82784/149‑1 (MQ=255)
                                                                                                                                tttttCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTGGTGGGCTCACTTGATGACGGCGGACATCCGCCGGCATCCCGGAACGGGCGTCCATCAGAGGTTGCATGCCGTCCGCGTGAGTTGGtgcggg              >  2:1534983/1‑146 (MQ=255)
                                                                                                                                  tttCCCCGCGACATGCCGTGGCGCGGGGATTTTTTTAACCGGGCTTTAGAGCGTTTGTTGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCAt            <  1:754972/149‑1 (MQ=255)
                                                                                                                                     ccccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCCCTTCAGCGTCAGTTTGTGCTGCATCtt         >  1:877906/1‑149 (MQ=255)
                                                                                                                                      cccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCCTCACTTCAGCGTCAGTTTGTGCTGCATCttt        >  1:729871/1‑149 (MQ=255)
                                                                                                                                       ccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa       >  2:473901/1‑149 (MQ=255)
                                                                                                                                       ccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTATACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTCGATCACTTCAGCGTCAGTTTGTGCTGCTTCTTTa       >  2:1092737/1‑149 (MQ=255)
                                                                                                                                       ccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGCCGGCGGACATCCACCGGCATCCCGGCACGGACATCCATCACATGTTGCCTCACTTCAGCGTGAGTTTGTGCTGCATCTTTa       >  2:123104/1‑149 (MQ=255)
                                                                                                                                         gcgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaaa     >  1:801721/1‑149 (MQ=255)
                                                                                                                                         gcgACATGCCGGGTCGGGGGGATTTTTTTATCCGGGCTTTACAGCGTGTGTGGGCTCCCTTGATGACGGCGGACATCCACCGGCATCCCGTAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaaa     <  2:845874/149‑1 (MQ=255)
                                                                                                                                          cgACATGCCGTGTCGCGGGGATTTTTTTATCCGTGCTTTACAGCGTTTGTGGGCTGACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATAACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAAt    <  1:911454/149‑1 (MQ=255)
                                                                                                                                          cgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTTTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAAt    <  2:872055/149‑1 (MQ=255)
                                                                                                                                          cgACATGCCGTGCCGGGGGGATTTTTTTATCCGGGCTTTACAGCGGTTGGGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAAt    <  2:428197/149‑1 (MQ=255)
                                                                                                                                            aCATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTTAATGa  >  2:451638/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGA  >  NC_000913/1169237‑1169525

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: