Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,491,995 A→T 6.8% intergenic (+55/+75) yhdZ → / ← degS amino acid ABC transporter ATP‑binding protein/serine endoprotease DegS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,491,9950AT 6.8% 602.1 / 12.4 190intergenic (+55/+75)yhdZ/degSamino acid ABC transporter ATP‑binding protein/serine endoprotease DegS
Reads supporting (aligned to +/- strand):  ref base A (93/78);  new base T (10/5);  total (105/85)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

GGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGC  >  NC_002947/1491848‑1492142
                                                                                                                                                   |                                                                                                                                                   
   gagaTTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCt                                                                                                                                                <  2:672183/150‑1 (MQ=255)
   gagaTTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCt                                                                                                                                                <  1:3856196/150‑1 (MQ=255)
   gagaTTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCt                                                                                                                                                <  1:4337882/150‑1 (MQ=255)
   gagaTTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCt                                                                                                                                                <  2:5265747/150‑1 (MQ=255)
    agaTTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTg                                                                                                                                               >  2:6945757/1‑150 (MQ=255)
    agaTTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTg                                                                                                                                               <  1:6263876/150‑1 (MQ=255)
     gaTTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTgg                                                                                                                                              >  1:6767971/1‑150 (MQ=255)
      aTTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGGTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGc                                                                                                                                             <  1:7070966/150‑1 (MQ=255)
         gtgGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCttt                                                                                                                                          <  1:1768293/150‑1 (MQ=255)
         gtgGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCttt                                                                                                                                          <  1:3199888/150‑1 (MQ=255)
         atgcgcaacgggcagcagttgaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCttt                                                                                                                                          >  2:4119789‑M1/106‑150 (MQ=255)
          tgGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTg                                                                                                                                         <  1:19998/150‑1 (MQ=255)
          tgGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTg                                                                                                                                         <  1:578723/150‑1 (MQ=255)
             gcaacgggcagcagttgaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggctgcgcggcaGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAg                                                                                                                                      <  1:2293210‑M1/49‑1 (MQ=255)
              gCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTTGCGCTGGCTTTGGAGc                                                                                                                                     >  2:5641504/1‑150 (MQ=255)
              gCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCGCCCACAGGTTTAGCGCTGGCGTTGGAGc                                                                                                                                     >  1:4798915/1‑150 (MQ=255)
                aGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCt                                                                                                                                   <  2:5510788/150‑1 (MQ=255)
                  gCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGc                                                                                                                                 <  2:7309759/150‑1 (MQ=255)
                    cgcgCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGc                                                                                                                                     <  1:4895316/144‑1 (MQ=255)
                      cgcCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTaa                                                                                                                             <  1:4103376/150‑1 (MQ=255)
                        ccGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAAtt                                                                                                                           <  1:2114526/150‑1 (MQ=255)
                         cGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTc                                                                                                                          >  1:5241998‑M2/1‑149 (MQ=255)
                          ggATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcc                                                                                                                         <  2:1112302‑M2/150‑3 (MQ=255)
                          ggATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCCCACCCAAAGGTTTATCGCTGGCGTTTGTGGCGGCGTAAtttc                                                                                                                         >  1:6330938/1‑148 (MQ=255)
                           gATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTg                                                                                                                                         >  1:6877030/1‑133 (MQ=255)
                           gATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcct                                                                                                                        >  2:1348991‑M2/1‑147 (MQ=255)
                           gATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcct                                                                                                                        >  1:383050‑M2/1‑147 (MQ=255)
                           gATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcct                                                                                                                        >  2:6438828‑M2/1‑147 (MQ=255)
                            tgaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCtg                                                                                                                       >  1:133539‑M1/87‑150 (MQ=255)
                             tGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgt                                                                                                                      >  1:3315119‑M2/1‑145 (MQ=255)
                             tGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgt                                                                                                                      >  2:3506044‑M2/1‑145 (MQ=255)
                             tGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgt                                                                                                                      >  1:228533‑M2/1‑145 (MQ=255)
                             tGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgt                                                                                                                      >  2:6219267‑M2/1‑145 (MQ=255)
                             tGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgt                                                                                                                      >  1:2464640‑M2/1‑145 (MQ=255)
                             tGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgt                                                                                                                      >  1:5744407‑M2/1‑145 (MQ=255)
                             tGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGGTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgt                                                                                                                      >  1:5815956‑M2/1‑145 (MQ=255)
                             gaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCtgt                                                                                                                      >  2:456550‑M1/86‑150 (MQ=255)
                             gaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCtgt                                                                                                                      >  1:4609981‑M1/86‑150 (MQ=255)
                             gaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCtgt                                                                                                                      >  2:681384‑M1/86‑150 (MQ=255)
                             gaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCtgt                                                                                                                      >  2:6108670‑M1/86‑150 (MQ=255)
                             gaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCtgt                                                                                                                      >  2:4999916‑M1/86‑150 (MQ=255)
                             gaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCtgt                                                                                                                      >  1:2831478‑M1/86‑150 (MQ=255)
                             gaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATTACGCTGGCTTTGGAGCCTGCGTAATTCCtgt                                                                                                                      >  1:4164485‑M1/86‑150 (MQ=255)
                              gACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtg                                                                                                                     >  1:5161938‑M2/1‑144 (MQ=255)
                              gACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtg                                                                                                                     >  2:6101233‑M2/1‑144 (MQ=255)
                              aagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCtgtg                                                                                                                     >  1:4969433‑M1/85‑150 (MQ=255)
                              aagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCtgtg                                                                                                                     >  2:2707734‑M1/85‑150 (MQ=255)
                              aagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcaGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCtgtg                                                                                                                     >  1:471158‑M1/85‑150 (MQ=255)
                               agctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTgg                                                                                                                    >  1:5959445‑M1/84‑150 (MQ=255)
                               agctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgtggcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTgg                                                                                                                    >  1:3296381‑M1/84‑150 (MQ=255)
                               aCTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGTCAGCCCTTCGCGTGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgg                                                                                                                    >  1:6439031‑M2/1‑143 (MQ=255)
                               aCTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgg                                                                                                                    >  2:3731120‑M2/1‑143 (MQ=255)
                               aCTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgg                                                                                                                    >  2:2807614‑M2/1‑143 (MQ=255)
                                gctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTggg                                                                                                                   >  1:441333‑M1/83‑150 (MQ=255)
                                 ttcttcGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtggga                                                                                                                  >  2:5698656‑M2/1‑141 (MQ=255)
                                 ttcttcGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtggga                                                                                                                  >  2:1699372‑M2/1‑141 (MQ=255)
                                 ctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGa                                                                                                                  >  2:5712464‑M1/82‑150 (MQ=255)
                                 ctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGa                                                                                                                  >  2:5133927‑M1/82‑150 (MQ=255)
                                 ctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGa                                                                                                                  >  1:5036895‑M1/82‑150 (MQ=255)
                                  tgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggctgcgcggcaGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAg                                                                                                                 >  2:2292957‑M1/81‑150 (MQ=255)
                                  tgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggctgcgcggcaGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAg                                                                                                                 >  2:6044033‑M1/81‑150 (MQ=255)
                                   cttcGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagc                                                                                                                >  1:2992651‑M2/1‑139 (MQ=255)
                                     tcGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgg                                                                                                              >  1:14765‑M2/1‑137 (MQ=255)
                                       gACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATTACGCAGGCTCCAAAGCCAGCGCTAAAcctgtgggagcgggc                                                                                                            >  1:1105096‑M2/1‑135 (MQ=255)
                                       gACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCACCCACCGGTTTTGCGCTGGGGTTTGTGCCGGCGTAATTgcggtgggggcgggc                                                                                                            <  1:3513786‑M2/150‑16 (MQ=255)
                                         cgtggtggggcggcggccgccaccggcgccggctgccagccagtaagataaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGTTTTAGCTCTGTCTTTGGAGCCTGATTAATTCCTGTGGGAGCGGGCGt                                                                                                          >  1:950342‑M1/74‑150 (MQ=255)
                                           ggCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggcgtgc                                                                                                        >  2:3010392‑M2/1‑131 (MQ=255)
                                               gcgcAGTGCCCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTGGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCGGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggcgtgcccgc                                                                                                    >  1:4310868‑M2/1‑127 (MQ=255)
                                                       aCCGGACCAAGTTGTGCCTCCGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGcc                                                                                            <  2:3229313/150‑1 (MQ=255)
                                                                     ccggctgccagccaggaagagaaataaacccagggctgcgcggcaGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCtgt                                                                                                                      <  2:5447665‑M1/65‑1 (MQ=255)
                                                                     ccggctgccagccaggaagagaaataaacccagggctgcgcggcaGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCtgt                                                                                                                      >  1:5447665‑M1/46‑110 (MQ=255)
                                                                         tCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGtt                                                                          >  2:2293283/1‑150 (MQ=255)
                                                                            gCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCGCCCACAGGGTTAGCGCAGGCTTTGGAGCCGGCGGAAATCCGGTGGGAGCGGGCGTGCCCCCGCAGGGCCGCCACGCGGCCCGGGGGTTAt                                                                       >  1:5378492/1‑150 (MQ=25)
                                                                             ccAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTCGCGCTGGCGTTGGAGCCTGCGTCGTTCCTGTGGGCGCGGGCGTGCGCGCGAAGGGCCGCACCGCGGCGCTGGGGTTCtt                                                                      >  2:4255351/1‑150 (MQ=38)
                                                                                    cgtctttgttgCTGTGGGGGCGGCGCGGCTGCCGTTCGGGGGCCCGCCGGCTCCCACTGTTTTCGCGGTGGCTTTGGAGCGTGCGTAATTCCGGTGGGCGGGGGGGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTc                                                               <  1:6500315/145‑1 (MQ=255)
                                                                                     tGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTATCGCTGGCTTTGGAGCCTGCGTAATTCCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTcc                                                              >  2:5399037/1‑150 (MQ=255)
                                                                                          aaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCg                                                                                                               >  1:4426272‑M1/25‑96 (MQ=255)
                                                                                          aaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCg                                                                                                               <  2:4426157‑M1/72‑1 (MQ=255)
                                                                                              aaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATTACGCTGGCTTTGGAGCCTGCGTAATTCCtgt                                                                                                                      <  2:2184489‑M1/65‑1 (MQ=255)
                                                                                                                           GGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCGTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCAc                        <  2:2171125‑M1/150‑1 (MQ=255)
                                                                                                                           GGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCGGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCAc                        <  2:6911319‑M1/150‑1 (MQ=255)
                                                                                                                           GGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCGTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCGGGGTTTATTTCTCTTCCTGGCTGGCCGCCGGCGCCGGGGGCGGCCGCAGCCCCAc                        <  2:4285659‑M1/150‑1 (MQ=255)
                                                                                                                                 GCCCGCTCCCACTGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGc                  >  2:5846015‑M1/1‑150 (MQ=255)
                                                                                                                                      CTCCCACCGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCa             >  2:5315510‑M1/1‑150 (MQ=255)
                                                                                                                                         ccACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggcgtgcccgcgaagggctgccgcgcagcccaaacagcatcaatgcaggatctggctgaggaacaacttggtccggtcactgcgcggccggtcgaagaagt          <  2:196685‑M2/150‑114 (MQ=255)
                                                                                                                                         ccACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggcgtgcccgcgaagggctgccgcgcagcccaaacagcatcaatgcaggatctggctgaggaacaacttggtccggtcactgcgcggccggtcgaagaagt          <  1:6786711‑M2/150‑114 (MQ=255)
                                                                                                                                         CCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCt          <  1:2877063‑M1/150‑1 (MQ=255)
                                                                                                                                         CCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCt          <  1:1799833‑M1/150‑1 (MQ=255)
                                                                                                                                         CCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCt          <  1:7330122‑M1/150‑1 (MQ=255)
                                                                                                                                         CCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGAGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCt          <  1:5302823‑M1/150‑1 (MQ=255)
                                                                                                                                          ccccGGTTTCGCGCGGGCTTTGGAGCCTGCGAAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCtt         >  1:5306432‑M1/5‑150 (MQ=255)
                                                                                                                                          cacaGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggcgtgcccgcgaagggctgccgcgcagcccaaacagcatcaatgcaggatctggctgaggaacaacttggtccggtcactgcgcggccggtcgaagaagtc         <  1:279869‑M2/150‑115 (MQ=255)
                                                                                                                                          cacaGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggcgtgcccgcgaagggccgccgcgcagcccaaacagcatcaatgcaggatctggctgaggaacaacttggtccggtcactgcgcggccggtcgaagaagtc         <  1:6174282‑M2/150‑115 (MQ=255)
                                                                                                                                          CACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCtt         <  1:2158171‑M1/150‑1 (MQ=255)
                                                                                                                                          CACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCtt         <  1:2206025‑M1/150‑1 (MQ=255)
                                                                                                                                           acaGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggcgtgcccgcgaagggccgccgcgcagcccaaacagcatcaat                                                                   >  2:4853224‑M2/1‑35 (MQ=255)
                                                                                                                                              GGTTTTGCGCTGGCTTTGGCGCCTGCTTATTGCCTGTGGGTGGGGGGGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAAc     >  1:6945757‑M1/1‑150 (MQ=255)
                                                                                                                                                tttAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGGTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAActg   >  1:3112387/1‑150 (MQ=255)
                                                                                                                                                 ttAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAActgc  >  2:3978814/1‑150 (MQ=255)
                                                                                                                                                   aGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTc                                                               <  1:7301099/87‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGC  >  NC_002947/1491848‑1492142

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