Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,169,382 C→G 5.5% intergenic (+52/+30) bhsA → / ← ldtC DUF1471 domain‑containing multiple stress resistance outer membrane protein BhsA/L,D‑transpeptidase YcfS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,169,3820CG 5.5% 245.6 / 4.3 73intergenic (+52/+30)bhsA/ldtCDUF1471 domain‑containing multiple stress resistance outer membrane protein BhsA/L,D‑transpeptidase YcfS
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base G (2/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

CGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCA  >  NC_000913/1169236‑1169530
                                                                                                                                                  |                                                                                                                                                    
cgcAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCgt                                                                                                                                                   >  2:136042/1‑150 (MQ=255)
        gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCggg                                                                                                                                           <  1:441361/150‑1 (MQ=255)
        gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCggg                                                                                                                                           <  1:933071/150‑1 (MQ=255)
        gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCggg                                                                                                                                           <  2:28907/150‑1 (MQ=255)
         cGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCgggg                                                                                                                                          <  1:893529/150‑1 (MQ=255)
             tGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttt                                                                                                                                      >  1:473462/1‑150 (MQ=255)
              gagaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttt                                                                                                                                     >  1:1096286/1‑150 (MQ=255)
               agaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttt                                                                                                                                    <  2:400011/150‑1 (MQ=255)
               agaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttt                                                                                                                                    <  1:182919/150‑1 (MQ=255)
                gaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttttt                                                                                                                                   >  2:541251/1‑150 (MQ=255)
                gaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttttt                                                                                                                                   <  1:535930/150‑1 (MQ=255)
                            aTCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCttt                                                                                                                       >  2:709076/1‑150 (MQ=255)
                              cTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTAc                                                                                                                     >  1:1391721/1‑150 (MQ=255)
                                    gTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGtt                                                                                                               <  2:473462/150‑1 (MQ=255)
                                    gTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGtt                                                                                                               >  1:399301/1‑150 (MQ=255)
                                        tACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTtgtg                                                                                                           <  2:846536/150‑1 (MQ=255)
                                         aCTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTgg                                                                                                          >  1:270895/1‑150 (MQ=255)
                                          cTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACGGCGTTTGTggg                                                                                                         >  2:195612/1‑150 (MQ=255)
                                          cTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTggg                                                                                                         >  2:985539/1‑150 (MQ=255)
                                             cTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTc                                                                                                      >  2:542322/1‑150 (MQ=255)
                                                 aaCCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACtt                                                                                                  >  2:290033/1‑150 (MQ=255)
                                                 aaCCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACtt                                                                                                  >  2:673659/1‑150 (MQ=255)
                                                  aCCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCGGGGa                                                                                                                                         <  2:74054/110‑1 (MQ=255)
                                                  aCCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCGGGGa                                                                                                                                         >  1:74054/1‑110 (MQ=255)
                                                   ccGGTGCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                                                                <  1:185783/150‑1 (MQ=255)
                                                    cGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgat                                                                                               <  2:1267162/150‑1 (MQ=255)
                                                    cGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgat                                                                                               <  1:154012/150‑1 (MQ=255)
                                                            aTACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGAcggcgg                                                                                       <  2:399301/150‑1 (MQ=255)
                                                              aCCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGAc                                                                                     >  1:1237348/1‑150 (MQ=255)
                                                              aCCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGAc                                                                                     <  2:828906/150‑1 (MQ=255)
                                                               cccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCgg                                                                                                                                            >  1:73400/1‑94 (MQ=255)
                                                               cccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCgg                                                                                                                                            <  2:73400/94‑1 (MQ=255)
                                                                 cTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATc                                                                                  <  2:270895/150‑1 (MQ=255)
                                                                  tCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCggg                                                                                                                                           >  2:518580/1‑92 (MQ=255)
                                                                  tCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCggg                                                                                                                                           <  1:518580/92‑1 (MQ=255)
                                                                  tCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATcc                                                                                 <  2:616255/150‑1 (MQ=255)
                                                                   ccATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCa                                                                                >  1:1043506/1‑150 (MQ=255)
                                                                      tGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCg                                                                             >  1:706600/1‑150 (MQ=255)
                                                                           cagcagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATc                                                                        >  2:401552/1‑150 (MQ=255)
                                                                            agcagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGCCGGTGGATGTCCGCCGTCATca                                                                       <  1:136042/150‑2 (MQ=255)
                                                                              cagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCTTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCg                                                                     <  2:1366514/150‑1 (MQ=255)
                                                                              cagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGCCGGTGGATGTCCGCCGTCATcaag                                                                     <  2:603095/150‑4 (MQ=255)
                                                                                    tttATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACg                                                                >  2:1126971/1‑149 (MQ=255)
                                                                                    tttATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACg                                                                <  1:1126971/149‑1 (MQ=255)
                                                                                     ttATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa                                                              <  1:484139/150‑1 (MQ=255)
                                                                                     ttATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGa                                                              <  2:654600/150‑1 (MQ=255)
                                                                                          aataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatcc                                                         <  2:1096286/150‑1 (MQ=255)
                                                                                              aGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATca                                                     <  1:541251/150‑1 (MQ=255)
                                                                                                  ttAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATg                                                 <  1:1323559/150‑1 (MQ=255)
                                                                                                   tAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGt                                                <  2:516753/150‑1 (MQ=255)
                                                                                                    aaCCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGtt                                               >  2:35220/1‑150 (MQ=255)
                                                                                                      cccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGc                                             >  2:302109/1‑150 (MQ=255)
                                                                                                          cATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatcca                                                        >  1:70942/1‑135 (MQ=255)
                                                                                                          cATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatcca                                                        <  2:70942/135‑1 (MQ=255)
                                                                                                              aaTGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTc                                     >  1:300767/1‑150 (MQ=255)
                                                                                                                  ccTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATc                                          <  2:943921/141‑1 (MQ=255)
                                                                                                                  ccTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATc                                          >  1:943921/1‑141 (MQ=255)
                                                                                                                       tACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGt                            >  1:1222043/1‑150 (MQ=255)
                                                                                                                           gctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt                        >  2:339265/1‑150 (MQ=255)
                                                                                                                            ctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTggg                                                                                                         >  2:605851/1‑68 (MQ=255)
                                                                                                                            ctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTggg                                                                                                         <  1:605851/68‑1 (MQ=255)
                                                                                                                            ctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgtg                       >  1:645426/1‑150 (MQ=255)
                                                                                                                            ctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgtg                       <  2:1391721/150‑1 (MQ=255)
                                                                                                                              gATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgct                     <  2:706600/150‑1 (MQ=255)
                                                                                                                                ttttttCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgctgc                   >  2:752002/1‑150 (MQ=255)
                                                                                                                                 tttttCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCa                  <  2:1043506/150‑1 (MQ=255)
                                                                                                                                    ttCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCt               <  1:290033/150‑1 (MQ=255)
                                                                                                                                       cccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTa            >  1:182378/1‑150 (MQ=255)
                                                                                                                                          gcgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAAt         <  2:1185415/150‑1 (MQ=255)
                                                                                                                                          gcgACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAAt         >  1:1105264/1‑150 (MQ=255)
                                                                                                                                               aTGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattg    <  1:709076/150‑1 (MQ=255)
                                                                                                                                                tGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattgc   >  2:702553/1‑150 (MQ=255)
                                                                                                                                                 gCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattgca  <  1:542322/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCA  >  NC_000913/1169236‑1169530

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: