Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,487,046 A→T 6.3% 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,487,0460AT 6.3% 651.3 / 16.1 217intergenic (+55/+75)yhdZ/degSamino acid ABC transporter ATP‑binding protein/serine endoprotease DegS
Reads supporting (aligned to +/- strand):  ref base A (89/105);  new base T (13/2);  total (107/110)
Fisher's exact test for biased strand distribution p-value = 2.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

AGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCT  >  NC_002947/1486897‑1487194
                                                                                                                                                     |                                                                                                                                                    
ccatcgaggtgatgcgcaacgggcagcagttgaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcaGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTa                                                                                                                                                      >  2:8923542‑M1/117‑150 (MQ=255)
  ggCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAgc                                                                                                                                                    >  2:7154273/1‑150 (MQ=255)
  ggCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAgc                                                                                                                                                    <  1:8295899/150‑1 (MQ=255)
    cGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAgcgc                                                                                                                                                  <  1:3108736/150‑1 (MQ=255)
      agaTTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTg                                                                                                                                                <  1:5010950/150‑1 (MQ=255)
           gtgGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCttt                                                                                                                                           <  1:1151717/150‑1 (MQ=255)
           gtgGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCttt                                                                                                                                           <  2:6545367/150‑1 (MQ=255)
            tgGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTg                                                                                                                                          <  2:7297190/150‑1 (MQ=255)
              gAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGa                                                                                                                                        <  2:4959914/150‑1 (MQ=255)
              gAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGa                                                                                                                                        <  1:5650644/150‑1 (MQ=255)
              gAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCACCCACAGGGGTTGCGCGGGCGTTggg                                                                                                                                        >  1:5762496/1‑149 (MQ=255)
              cgcaacgggcagcagttgaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGa                                                                                                                                        >  1:6110542‑M1/103‑150 (MQ=255)
                  aGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCt                                                                                                                                    <  1:651145/150‑1 (MQ=255)
                    gCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGTTGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGc                                                                                                                                  <  1:254330/150‑1 (MQ=255)
                    gCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTATCGCGGGCGTTGGTGCCTGc                                                                                                                                  >  1:4420283/1‑150 (MQ=255)
                     ccGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCg                                                                                                                                 >  1:5209844/1‑150 (MQ=255)
                        agcagttgaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTaa                                                                                                                              >  2:253576‑M1/93‑150 (MQ=255)
                         gcagttgaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATTACGCTGGCTTTGGAGCCTGCGTAAt                                                                                                                             >  2:5671347‑M1/92‑150 (MQ=255)
                          ccGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAAtt                                                                                                                            <  2:3412555/150‑1 (MQ=255)
                           cGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTc                                                                                                                           >  1:4718245‑M2/1‑149 (MQ=255)
                           cGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTc                                                                                                                           >  1:4718487‑M2/1‑149 (MQ=255)
                           cGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTc                                                                                                                           >  2:8383702‑M2/1‑149 (MQ=255)
                            gttgaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCTGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcc                                                                                                                          >  2:1113506‑M1/89‑150 (MQ=255)
                            gttgaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcc                                                                                                                          >  2:4365714‑M1/89‑150 (MQ=255)
                             gATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcct                                                                                                                         >  2:5704446‑M2/1‑147 (MQ=255)
                             gATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcct                                                                                                                         >  1:8199128‑M2/1‑147 (MQ=255)
                               tGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgt                                                                                                                       >  2:1214675‑M2/1‑145 (MQ=255)
                               tGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgt                                                                                                                       >  2:3212039‑M2/1‑145 (MQ=255)
                               tGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgt                                                                                                                       >  2:8765861‑M2/1‑145 (MQ=255)
                               ggaCGTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGGTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgt                                                                                                                       >  2:7428680‑M2/2‑145 (MQ=255)
                               gaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggacgagaaatacacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCtgt                                                                                                                       >  2:8396898‑M1/86‑150 (MQ=255)
                               gaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCtgt                                                                                                                       >  1:8000926‑M1/86‑150 (MQ=255)
                               gaagctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATTGCGCTGGCTTTGGAGCCTGCGTAATTCCtgt                                                                                                                       >  1:2499399‑M1/86‑150 (MQ=255)
                                 agctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTgg                                                                                                                     >  2:2731665‑M1/84‑150 (MQ=255)
                                 agctgatcgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATTACGCTGGCTTTGGAGCCTGCGTAATTCCTGTgg                                                                                                                     >  1:4223672‑M1/84‑150 (MQ=255)
                                 aCTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCAGTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgg                                                                                                                     >  2:1637500‑M2/1‑143 (MQ=255)
                                   ttcttcGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCATATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCCCGCCCGCACCCCCCGGTTTATCGCTGGCGGTGGAGGCAGCGTAAATcctgggggg                                                                                                                   <  1:7833798‑M2/150‑10 (MQ=255)
                                   ttcttcGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtggga                                                                                                                   <  1:5289377‑M2/150‑10 (MQ=255)
                                   ctgatcgccgaggtggggctgcggccgccaccggcgccgcctgccagccaggaagagaaataaacccagggccgcgttgcgGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGa                                                                                                                   >  2:7292634‑M1/82‑150 (MQ=255)
                                      ttcGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGGACGCCCGCTCCCACAGGTTTATCGCTTGGTTTTGTGCCTGCGTAATTgctttgggttcg                                                                                                                <  1:5780905‑M2/150‑13 (MQ=255)
                                      ttcGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcg                                                                                                                >  1:6486766‑M2/1‑138 (MQ=255)
                                        cgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgcggcaGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGa                                                                                                                   <  2:2178647‑M1/69‑1 (MQ=255)
                                        cgccgaggtggggctgcggccgccaccggcgccggctgccagccaggaagagaaataaacccagggccgcgcggcaGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGa                                                                                                                   >  1:2178647‑M1/77‑145 (MQ=255)
                                        cGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCCCGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTccggtgggagcggg                                                                                                              <  2:3381271‑M2/150‑15 (MQ=255)
                                        cGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATTACGCAGGCTCCAAAGCCAGCGCTAAAcctgtgggagcggg                                                                                                              >  2:1922608‑M2/1‑136 (MQ=255)
                                         gACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggc                                                                                                             >  2:5289377‑M2/1‑135 (MQ=255)
                                             ggCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggcgtgc                                                                                                         <  2:3008253‑M2/150‑20 (MQ=255)
                                             ggCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCGGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCCCTCCCCCAGGGTTAGCGCTGGCGTTGGGGCCGGCGTCATTcctgtgggggcgggcgtgc                                                                                                         <  2:2490078‑M2/150‑20 (MQ=255)
                                                 gcgcAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggcgtgcccgc                                                                                                     >  2:5789745‑M2/1‑127 (MQ=255)
                                                    cAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATTACGCAGGCTCCAAAGCCAGCGCTAAAcctgtgggagcgggcgtgcccgcgaa                                                                                                  >  1:7583186‑M2/1‑124 (MQ=255)
                                                      gtgccgggACAATGTTGTTCCTCAGCCAGATCCTGATTTGTTGCTGTTTGGGCTGCGCGGCTGCCCTTCGCGGGCACGCCCGCTCCCACTGTTTTAGCGCTGGCTTTGGAGCCTGCTTAATTCCTTTGGGCGCGGGCGTGCCCGCGAAgg                                                                                                <  1:6053437/144‑1 (MQ=255)
                                                      gTGACCGGACCACGTTGTTCCTCCGCCAGATCCTGCATTGATGCTGTTGGGGCTGCGCGGCAGCCCTGCGCGGGCACGCCCGCTCCCACCGGTTTAGCGCTGGCTTTGGAGCCGGCGTAATTCCTGTGGGCGCGGGCGTGCCCGCGAAgg                                                                                                <  2:4174769/150‑1 (MQ=38)
                                                                cAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgt                                                                                                                       <  1:7805610‑M2/117‑6 (MQ=255)
                                                                cAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgt                                                                                                                       >  2:7804669‑M2/1‑112 (MQ=255)
                                                                        tCCTCAGCCCGATCCTGCATTGAGGCGGTTGGGGCGGCGCGGCAGCCCTTCGCGGGCACGCCCGCGCCCACAGGTTTAGCGCTGGCTTTGGCGCCTGCGTAATGCCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTgg                                                                              <  2:4979335/150‑1 (MQ=38)
                                                                         ccTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCGGCGTAAATCCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCgggg                                                                             >  2:1106821/1‑150 (MQ=255)
                                                                                cAGATCCTGCAGTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTAttt                                                                      <  2:2481905/150‑1 (MQ=255)
                                                                                            tGATGCTGTTTGGGCTGCGCGGCCGCCCTTCGCGGGCACGCCCGCTCCCACCGGTTTAGCGGTGGCTTTGGAGCCTGCGTAATGCGTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCctggct                                                          <  2:6207383/150‑1 (MQ=39)
                                                                                             gATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGTCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCctggctg                                                         >  1:5563654/1‑150 (MQ=255)
                                                                                              aTGCTGTTTGGGCTGCGCGGCTGCCCTTCGGGGGCACGCCCGCTCCCACCGTTTTTGCGCTGGGTTTGGTGCCTGCGTATTTCGTGTGGGTGGGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCctggctgg                                                        <  1:7370843/150‑1 (MQ=255)
                                                                                              aTGCTGTTTGGGCTGCGCGGCCGCCCTTCGCGGGCACGCCCGCTCCCCCAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCctggctgg                                                        <  2:5126332/150‑1 (MQ=255)
                                                                                              aTGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCCACGCGGCCCTGGGTTTATTTCTCTTCctggctgg                                                        >  2:4442752/1‑150 (MQ=255)
                                                                                                gcggTTGGGGCTGCGCGGCAGCCCTTGGGGGGCCCGCCCGCTCCCACTGTTTTCGGGGGGGCTTTGGAGCCGGCGTAATTCGTGGGGGTGGGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCa                                                      <  1:8833609/147‑1 (MQ=255)
                                                                                                gCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTTGCGCTGGCTTTGGAGCCTGCGTAAATCCTGTGGGAGCGGGCGGGCCCGCGAAGGGCCGCAACCCGGCCCAGGGTTTATTTCTCTTCCAGGCTGGCa                                                      >  1:2981169/1‑150 (MQ=38)
                                                                                                    tttGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCCCGCCCACAGGTTTTGCGCTGGGGTTGGGGCCGGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCCACCCGGCCCGGGGTTTTTTTGTCGTCCTGGCTGGCAGCCg                                                  >  1:3652479/1‑150 (MQ=37)
                                                                                                      ggggCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGgc                                                <  2:7176807/149‑1 (MQ=255)
                                                                                                                  ctGCCGTTCGGGGGCACGCCCGCTCCCACTGTTTTTGCGGGGGCTTTGGTGCCGGCTTAATTCCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTggcggc                                    >  1:1059375‑M1/3‑150 (MQ=255)
                                                                                                                    GCCCGTCGCGGGCACGCCCGCTCCCACAGGTTTCGCGCGGGCTTTGGAGCCTGCGTATTTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCg                                  >  2:4580600‑M1/1‑150 (MQ=255)
                                                                                                                        tcCGCGGGCCCGCCCGCCCCACAGTTTTAGCGCTGGCTTGGGAGCCTGCGTAATGCCGGGGGGCGCGGGCGCGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGcc                             >  1:1585032‑M1/3‑150 (MQ=255)
                                                                                                                         TCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGcc                             >  2:6478018‑M1/1‑150 (MQ=255)
                                                                                                                         TCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCGGCGTAATTCCTGTGGGAGCGGGCGGGCCCGCGAAGGGCCGCAACGCGGCCCAGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGcc                             <  1:2825915‑M1/150‑1 (MQ=255)
                                                                                                                            CGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCa                          <  2:7577124‑M1/150‑1 (MQ=255)
                                                                                                                                 ACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCGGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCg                     <  1:1747185‑M1/150‑1 (MQ=255)
                                                                                                                                     CCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGa                 <  1:1872316‑M1/150‑1 (MQ=255)
                                                                                                                                        ccCCCACTGTTTTCGCGGGGGATTTGGCGCCTGCTTATTTCGTGTGGGCGGGGGCGTCCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCa              >  1:2248148‑M1/3‑150 (MQ=255)
                                                                                                                                         TCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAg             <  2:6328031‑M1/150‑1 (MQ=255)
                                                                                                                                           ccACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggcgtgcccgcgaagggctgccgcgcagcccaaacagcatcaatgcaggatctggctgaggaacaacttggtccggtcactgcgcggccggtcgaagaagt           <  2:8053502‑M2/150‑114 (MQ=255)
                                                                                                                                           ccACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggcgtgcccgcgaagggctgccgcgcagcccaaacagcatcaatgcaggatctggctgaggaacaacttggtccggtcactgcgcggccggtcgaagaagt           <  1:9070188‑M2/150‑114 (MQ=255)
                                                                                                                                           ccACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggcgtgcccgcgaagggctgccgcgcagcccaaacagcatcaatgcaggatctggctgaggaacaacttggtccggtcactgcgcggccggtcgaagaagt           <  2:1464122‑M2/150‑114 (MQ=255)
                                                                                                                                           ccACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggcgtgcccgcgaagggccgccgcgcagcccaaacagcatcaatgcaggatctggctgaggaacaacttggtccggtcactgcgcggccggtcgaagaagt           <  1:4056884‑M2/150‑114 (MQ=255)
                                                                                                                                           CCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCt           <  1:1325908‑M1/150‑1 (MQ=255)
                                                                                                                                           CCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCt           <  1:550619‑M1/150‑1 (MQ=255)
                                                                                                                                            cacaGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggcgtgcccgcgaagggctgccgcgcagcccaaacagcatcaatgcaggatctggctgaggaacaacttggtccggtcactgcgcggccggtcgaagaagtc          <  1:4589725‑M2/150‑115 (MQ=255)
                                                                                                                                            CACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGctg          <  1:6958061‑M1/150‑2 (MQ=255)
                                                                                                                                            CACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCtt          <  2:5890219‑M1/150‑1 (MQ=255)
                                                                                                                                            CACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCtt          <  1:4840921‑M1/150‑1 (MQ=255)
                                                                                                                                            CACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCtt          <  2:141304‑M1/150‑1 (MQ=255)
                                                                                                                                             ACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGTTCAGCTTc         <  2:4036959‑M1/150‑1 (MQ=255)
                                                                                                                                             ACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTc         <  1:8717135‑M1/150‑1 (MQ=255)
                                                                                                                                             ACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTc         <  1:2133304‑M1/150‑1 (MQ=255)
                                                                                                                                             ACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTc         <  2:125600‑M1/150‑1 (MQ=255)
                                                                                                                                             ACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTc         <  1:2606114‑M1/150‑1 (MQ=255)
                                                                                                                                             ACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTc         <  1:2301780‑M1/150‑1 (MQ=255)
                                                                                                                                               AGGTTTCGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCaa       >  1:4485252‑M1/1‑150 (MQ=255)
                                                                                                                                               AGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCctggctg                                                         >  1:3290790‑M1/1‑100 (MQ=255)
                                                                                                                                               AGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCaa       <  2:5441482‑M1/150‑1 (MQ=255)
                                                                                                                                                 gTTTAGCGCTGGCTTTGGAGCCTGCGTAATTcctgtgggagcgggcgtgcccgcgaagggccgcaacgcggcccaaacagcatcaatgcaggatctggctgaggaacaacttggtccggtcactgcgcggccggtcgaagaagtcatccg     <  1:3154674‑M2/150‑120 (MQ=255)
                                                                                                                                                  tttAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAActg    >  2:736422/1‑150 (MQ=255)
                                                                                                                                                    tAGCGCTGGCTTTGGAGCCTGCGCTAAAcctgtgggagcgggcgtgcccgcgaagggctgccgcgcagcccaaacagcatcaatgcaggatctggctgaggaacaacttggtccggtcactgcgcggccggtcgaagaagtcatccggcg  <  1:8333444‑M2/150‑123 (MQ=255)
                                                                                                                                                     aGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGc                   <  2:6102926/132‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCT  >  NC_002947/1486897‑1487194

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: