Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,169,384 G→A 18.8% 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,3840GA18.8% 218.1 / 39.3 101intergenic (+54/+28)bhsA/ldtCbiofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base G (35/47);  new base A (10/9);  total (45/56)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

Only 100 of 104 total aligned reads displayed.

ACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATG  >  NC_000913/1169286‑1169483
                                                                                                  |                                                                                                   
   ggTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTGCCCCGCGACATGCCGTGTc                                                                                                 <  2:2156036/100‑1 (MQ=255)
     tCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTcgc                                                                                               >  2:2275515/1‑100 (MQ=255)
       cGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCgg                                                                                             >  2:141920/1‑100 (MQ=255)
       cGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCgg                                                                                             <  1:1350722/100‑1 (MQ=255)
       cGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                             >  1:1051058/1‑100 (MQ=255)
       cGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                             <  2:1945656/100‑1 (MQ=255)
        gAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCggg                                                                                            <  1:2266957/100‑1 (MQ=255)
        gAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCggg                                                                                            <  1:1250150/100‑1 (MQ=255)
        gAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCggg                                                                                            <  2:656511/100‑1 (MQ=255)
        gAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCggg                                                                                            <  2:1279464/100‑1 (MQ=255)
        gAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCggg                                                                                            <  1:1774167/100‑1 (MQ=255)
         aaTACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgggg                                                                                           <  1:1303727/100‑1 (MQ=39)
         aaTACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgggg                                                                                           <  1:2411859/100‑1 (MQ=39)
         aaTACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgggg                                                                                           <  2:323728/100‑1 (MQ=39)
         aaTACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgggg                                                                                           <  1:1453546/100‑1 (MQ=39)
           tACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAt                                                                                         <  1:1377501/100‑1 (MQ=255)
           tACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAt                                                                                         >  1:858259/1‑100 (MQ=255)
             cccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttt                                                                                       <  1:141922/100‑1 (MQ=255)
             cccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCGGGGaaaa                                                                                       <  2:1051050/100‑4 (MQ=38)
              ccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttt                                                                                      >  2:1033729/1‑100 (MQ=255)
              ccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttt                                                                                      >  1:858439/1‑100 (MQ=255)
              ccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttt                                                                                      >  2:183846/1‑100 (MQ=255)
              ccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCAACATGCCGTGTCGCGGGGAttt                                                                                       >  2:1404326/1‑99 (MQ=255)
                tCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttttt                                                                                    <  2:1788302/100‑1 (MQ=255)
                tCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttttt                                                                                    <  1:2412873/100‑1 (MQ=255)
                 ccATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGTCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttttt                                                                                   >  1:2500154/1‑100 (MQ=255)
                  cATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttttt                                                                                   >  2:1651976/1‑99 (MQ=255)
                  cATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTa                                                                                  >  1:1178676/1‑100 (MQ=255)
                     ggAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATcc                                                                               <  1:1926682/100‑1 (MQ=255)
                        aCAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCggg                                                                            <  1:719898/100‑1 (MQ=255)
                         cagcagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGc                                                                           <  1:1807354/100‑1 (MQ=255)
                          agcagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCt                                                                          <  1:1479110/100‑1 (MQ=255)
                             agTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTa                                                                       <  2:810306/100‑1 (MQ=255)
                              gTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTAc                                                                      >  1:2063439/1‑100 (MQ=255)
                                aaTTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAg                                                                    >  2:830672/1‑100 (MQ=255)
                                  tttATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCg                                                                  >  1:1251395/1‑100 (MQ=255)
                                  tttATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCCCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCg                                                                  >  2:2533042/1‑100 (MQ=255)
                                    tataAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGtt                                                                >  1:1450693/1‑100 (MQ=255)
                                    tataAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGtt                                                                >  1:351261/1‑100 (MQ=255)
                                        aataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTtgtg                                                            <  2:1276407/100‑1 (MQ=255)
                                        aataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTtgtg                                                            <  2:911955/100‑1 (MQ=255)
                                         ataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTgg                                                           <  1:1850007/100‑1 (MQ=255)
                                          taaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTggg                                                          <  2:2005275/100‑1 (MQ=255)
                                          taaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTggg                                                          <  2:1251387/100‑1 (MQ=255)
                                          taaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTggg                                                          <  2:532783/100‑1 (MQ=255)
                                           aaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGc                                                         <  1:225907/100‑1 (MQ=255)
                                               aTTAACCCTCGTTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCAc                                                     <  2:1678020/100‑1 (MQ=255)
                                               aTTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCAc                                                     <  2:436062/100‑1 (MQ=255)
                                                ttAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACt                                                    >  1:428323/1‑100 (MQ=255)
                                                ttAACCCTCATTAATGCCTGCTACTGCTGATTTTTCCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACCGCGATTGTGGGCTCAca                                                    <  2:2300744/100‑2 (MQ=255)
                                                   aCCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                 <  1:661425/100‑1 (MQ=255)
                                                    cccTCATTCATGCCTGCTACTGCTGCTTTTTTCCCCGCGACATGCCGTGTCGCGGGTGTTTTTTTATACGTGCTGTACCGCCTTTGTGGGATCACTtgct                                                >  2:524696/1‑98 (MQ=255)
                                                    cccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgat                                                <  1:114377/100‑1 (MQ=255)
                                                     ccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatg                                               >  1:1393380/1‑100 (MQ=255)
                                                     ccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatg                                               >  1:69630/1‑100 (MQ=255)
                                                      cTCATTAATGCCTGCTACTGCTGATTTTTTTCTTTCTTCATGCCGTGTCGCGGGGATTTTTTTATTCTTTTTTTTTAGCGTTTGTGGGCTCAttt                                                   >  2:2469282/1‑95 (MQ=255)
                                                      cTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatga                                              <  1:667509/100‑1 (MQ=255)
                                                       tCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAc                                             >  1:72806/1‑100 (MQ=255)
                                                        cATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAc                                             >  2:1643895/1‑99 (MQ=255)
                                                         aTTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAcgg                                           <  1:1182767/100‑1 (MQ=255)
                                                          ttAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAcggc                                          >  2:1994537/1‑100 (MQ=255)
                                                           tAATGCCTGCTACTGCTGATTTTTTTTTTTTTTTTTGCCGTGTCGCGGGGATTTTTTTATCCGTGCTTTACAGCGTTTGTGGGCTCACtttttcacgg                                           >  2:759262/1‑90 (MQ=255)
                                                            aaTGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAcggcgg                                        <  2:475338/100‑1 (MQ=255)
                                                              tGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGAc                                      >  1:2483673/1‑100 (MQ=255)
                                                               gCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACa                                     <  2:1873700/100‑1 (MQ=255)
                                                               gCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGATTTACAGCGTTTGTGGGCTCACTTCATGACGGCGGACa                                     <  2:1581221/100‑1 (MQ=255)
                                                                 cTGGTACTGCTGAGTTTTTCCCGGGGACATGCCGTGTCGCGGGGATGTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCAATTGATGACGGCGGAc                                      <  2:1164104/97‑1 (MQ=255)
                                                                     tACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATcc                                  <  2:329514/97‑1 (MQ=255)
                                                                     tACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCAcc                               <  1:1432152/100‑1 (MQ=255)
                                                                       atgctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCgg                             <  1:1529238/99‑1 (MQ=255)
                                                                         gctgATTTTTTCCCCGCGACATGCCGTGTCGCAGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCa                           <  1:386920/100‑1 (MQ=255)
                                                                           tgATTTTTTCCCCGCGACATGCCGTGTCGCGGGTATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATc                         >  2:144559/1‑100 (MQ=255)
                                                                           tgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATc                         >  1:556839/1‑100 (MQ=255)
                                                                           tgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATc                         <  2:1947598/100‑1 (MQ=255)
                                                                            gATTTTTTCCCCGCGACATGCCGTGTCGTGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGCCGGTGGATGTCCGCCGTCATca                        <  2:1522718/100‑2 (MQ=37)
                                                                            gATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGCCGGTGGATGTCCGCCGTCATca                        <  2:50647/100‑2 (MQ=38)
                                                                            gATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATcc                        <  2:1716990/100‑1 (MQ=255)
                                                                               tttttCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCgg                     <  1:875765/100‑1 (MQ=255)
                                                                                ttttCCCCGCGACATGCCGTGTCGCGGGTATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGa                    <  1:144561/100‑1 (MQ=255)
                                                                                ttttCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGCCGGTGGATGTCCGCCGTCATcaagt                     <  1:128977/99‑5 (MQ=37)
                                                                                 tttCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGCCGGTGGATGTCCGCCGTCATcaagtga                   <  2:900391/100‑7 (MQ=25)
                                                                                  ttCCCCGCGACATGCCGGGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGCCGGTGGATGTCCGCCGTCATcaagtgag                  <  2:587012/100‑8 (MQ=21)
                                                                                    ccccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGTGCTTTACAGCGTTTGTGGGCTCACTTTATGACGGCGGGCATCCACCGGCATCCCGGAACg                 >  2:463088/1‑99 (MQ=255)
                                                                                    ccccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCATTTGATGACGGAGGACATGCACCCGCATccc                       >  2:845848/1‑93 (MQ=37)
                                                                                    ccccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACGGGCATCCCGGAACgg                >  1:2518299/1‑100 (MQ=38)
                                                                                    ccccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACgg                >  1:1504716/1‑100 (MQ=39)
                                                                                    ccccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACgg                >  2:275406/1‑100 (MQ=39)
                                                                                    ccccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACgg                >  1:163762/1‑100 (MQ=39)
                                                                                     cccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa               >  1:1472033/1‑100 (MQ=255)
                                                                                     cccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa               >  1:2154523/1‑100 (MQ=255)
                                                                                     cccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa               >  1:1471724/1‑100 (MQ=255)
                                                                                     cccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa               >  1:1647656/1‑100 (MQ=255)
                                                                                     cccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa               >  1:2095634/1‑100 (MQ=255)
                                                                                         cgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatc           <  1:2100498/100‑1 (MQ=255)
                                                                                         cgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatc           <  1:569412/100‑1 (MQ=255)
                                                                                              tGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATca      >  1:718374/1‑100 (MQ=255)
                                                                                               gCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATcac     >  2:206527/1‑100 (MQ=255)
                                                                                                ccGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATcaca    <  2:806865/100‑1 (MQ=255)
                                                                                                  gtgtCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATg  <  2:47454/100‑1 (MQ=255)
                                                                                                  gtgtCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATg  >  1:1511333/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
ACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATG  >  NC_000913/1169286‑1169483

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: