Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,169,382 C→G 18.7% 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,3820CG18.7% 247.9 / 52.1 101intergenic (+52/+30)bhsA/ldtCbiofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base C (36/46);  new base G (10/9);  total (46/55)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

Only 100 of 104 total aligned reads displayed.

GTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACA  >  NC_000913/1169283‑1169481
                                                                                                   |                                                                                                   
gTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGc                                                                                                     >  1:2005289/1‑100 (MQ=255)
gTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGc                                                                                                     >  2:225905/1‑100 (MQ=255)
   aCCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCgt                                                                                                  >  1:1788315/1‑100 (MQ=255)
      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                                                                                         <  2:323728/100‑1 (MQ=39)
            aaTACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgggg                                                                                         <  1:2411859/100‑1 (MQ=39)
            aaTACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgggg                                                                                         <  1:1453546/100‑1 (MQ=39)
            aaTACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgggg                                                                                         <  1:1303727/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)
                             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:532783/100‑1 (MQ=255)
                                             taaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTggg                                                        <  2:2005275/100‑1 (MQ=255)
                                             taaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTggg                                                        <  2:1251387/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:69630/1‑100 (MQ=255)
                                                        ccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatg                                             >  1:1393380/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                       <  2:1947598/100‑1 (MQ=255)
                                                                              tgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATc                       >  1:556839/1‑100 (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              >  2:275406/1‑100 (MQ=39)
                                                                                       ccccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACgg              >  1:1504716/1‑100 (MQ=39)
                                                                                       ccccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACgg              >  1:163762/1‑100 (MQ=39)
                                                                                        cccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa             >  1:1472033/1‑100 (MQ=255)
                                                                                        cccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa             >  1:1471724/1‑100 (MQ=255)
                                                                                        cccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa             >  1:2154523/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)
                                                                                                   |                                                                                                   
GTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACA  >  NC_000913/1169283‑1169481

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: