Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,169,384 G→A 5.9% intergenic (+54/+28) 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,3840GA 5.9% 276.9 / 6.6 103intergenic (+54/+28)bhsA/ldtCbiofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base G (47/48);  new base A (4/2);  total (51/52)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

Only 100 of 107 total aligned reads displayed.

GCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCAT  >  NC_000913/1169237‑1169531
                                                                                                                                                   |                                                                                                                                                   
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)
      aGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTcgcg                                                                                                                                              <  2:1101757/149‑1 (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:1180238/1‑149 (MQ=255)
                  gggCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTa                                                                                                                                  >  1:401033/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:1440771/1‑149 (MQ=255)
                                      ttACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTtg                                                                                                              >  2:586914/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)
                                                       ccGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAc                                                                                             <  1:1557428/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)
                                                               ccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCGCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATGTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACa                                                                                     >  1:77893/1‑149 (MQ=255)
                                                                   cATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTACCCCGCGACATGCGGGGTCGCGGGGTTTTTTTTTGCCGGGCTTAACGGCGTTTGTGGGCGCCATTGATGCCGGCGGCCAGccc                                                                                 >  1:470096/1‑148 (MQ=255)
                                                                     tGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCAcc                                                                               >  1:1272800/1‑149 (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)
                                                                               gTAATTTATAAATAAGCATTAACCCTCATTAATGCATGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGCGATTTTTTAATCCGGGCTTTACAGCGGTTGTGGGCTCACATTATGCCTGAGGAAATCCACCGGCATcaggg                                                                     >  2:189203/1‑145 (MQ=255)
                                                                                   tttATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACg                                                                 <  1:1078192/149‑1 (MQ=255)
                                                                                     tataAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa                                                               >  2:528701/1‑149 (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)
                                                                                                 ttAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACAt                                                   <  2:1515288/149‑1 (MQ=255)
                                                                                                    aCCCTCATTAATGCCTGCTACTGCTGATTTTATCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGtt                                                <  2:914518/149‑1 (MQ=255)
                                                                                                    aCCCTCATTAATGCCTGATACTGCTGATTTTTTCCCCTCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGtt                                                <  1:126774/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)
                                                                                                              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)
                                                                                                                                          cgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTTTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAAt          <  2:872055/149‑1 (MQ=255)
                                                                                                                                          cgACATGCCGTGCCGGGGGGATTTTTTTATCCGGGCTTTACAGCGGTTGGGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAAt          <  2:428197/149‑1 (MQ=255)
                                                                                                                                            aCATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTTAATGa        >  2:451638/1‑149 (MQ=255)
                                                                                                                                                  cGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattgcat  <  2:801721/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCAT  >  NC_000913/1169237‑1169531

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: