Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,084,420 T→G 13.0% intergenic (‑231/+153) sdaC ← / ← amaC serine:H+ symport permease/D‑lysine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,084,4200TG13.0% 536.6 / 6.5 203intergenic (‑231/+153)sdaC/amaCserine:H+ symport permease/D‑lysine aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (48/115);  new base G (29/3);  total (85/118)
Fisher's exact test for biased strand distribution p-value = 7.49e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 209 total aligned reads displayed.

GGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGC  >  NC_002947/4084271‑4084569
                                                                                                                                                     |                                                                                                                                                     
ggggCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCt                                                                                                                                                       <  2:6264125/150‑1 (MQ=255)
 gggCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCtt                                                                                                                                                      <  2:2669274/150‑1 (MQ=255)
 gggCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCtt                                                                                                                                                      <  1:2310295/150‑1 (MQ=255)
 gggCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCtt                                                                                                                                                      <  2:556852/150‑1 (MQ=255)
   gcgcGATTGTGCACGGATAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGc                                                                                                                                                    <  1:1930609/150‑1 (MQ=255)
   gacagcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGc                                                                                                                                                    >  2:6617128‑M1/134‑150 (MQ=255)
    cgcgATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTATTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCg                                                                                                                                                   >  2:3909321/1‑150 (MQ=255)
     cagcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGt                                                                                                                                                  >  1:461743‑M1/132‑150 (MQ=255)
      agcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTc                                                                                                                                                 >  2:3348665‑M1/131‑150 (MQ=255)
      agcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTc                                                                                                                                                 <  2:4456487‑M1/20‑1 (MQ=255)
          tgtgCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGTGCGTcgcgc                                                                                                                                             >  2:685817/1‑149 (MQ=255)
          tgtgCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTGGCGGCGCGa                                                                                                                                             >  2:487428/1‑150 (MQ=255)
           gtgCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGAt                                                                                                                                            <  1:2751477/150‑1 (MQ=255)
            tgCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATa                                                                                                                                           <  2:6796818/150‑1 (MQ=255)
             gCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAg                                                                                                                                          <  1:802936/150‑1 (MQ=255)
             gCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAg                                                                                                                                          <  2:5243074/150‑1 (MQ=255)
              cACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAgg                                                                                                                                         <  2:4599735/150‑1 (MQ=255)
                cGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGc                                                                                                                                       <  1:5417654/150‑1 (MQ=255)
                cGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGc                                                                                                                                       <  2:3171705/150‑1 (MQ=255)
                 ggAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGcc                                                                                                                                      <  2:1573424/150‑1 (MQ=255)
                       gTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAg                                                                                                                                <  1:2537605/150‑1 (MQ=255)
                        ttGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGc                                                                                                                               <  2:5353080/150‑1 (MQ=255)
                         tGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCg                                                                                                                              <  2:6323335/150‑1 (MQ=255)
                             cTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTTGCGGccc                                                                                                                          <  2:3970455/150‑1 (MQ=255)
                             cTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGccc                                                                                                                          <  2:5551672/150‑1 (MQ=255)
                               tGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCa                                                                                                                        <  2:4291747/150‑1 (MQ=255)
                                gCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAg                                                                                                                       <  1:3063416/150‑1 (MQ=255)
                                 cGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTATTTACAAGCGTAAGGCATTGCTGGTGTCGATCGGTGGATTGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGc                                                                                                                      <  1:3632149/150‑1 (MQ=255)
                                   aaGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGTGGGGCGGGCGGGGGGCGCCAAGCCGcccccccaa                                                                                                                    >  1:51988/1‑145 (MQ=255)
                                   aaGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATGGGGCCGCTTTGCGGccctatcgc                                                                                                                    <  2:14904/150‑7 (MQ=255)
                                   aaGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCaa                                                                                                                    <  1:3989482/150‑1 (MQ=255)
                                    aGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCaaa                                                                                                                                 <  1:1076597/136‑1 (MQ=255)
                                        cTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCGGGGGGGGCGCCAAACGGCCCCACCaaaaaaa                                                                                                               >  1:5559941/1‑145 (MQ=255)
                                          tttCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGc                                                                                                             <  1:2947989/150‑1 (MQ=255)
                                                   agCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCt                                                                                                    <  1:6844869/150‑1 (MQ=255)
                                                    gCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTg                                                                                                   <  1:1065173/150‑1 (MQ=255)
                                                     cGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGc                                                                                                  <  2:4044169/150‑1 (MQ=255)
                                                         tGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATc                                                                                              <  2:2746413/150‑1 (MQ=255)
                                                         tGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCGCTAGGGCCGCAAAGCGGGCCCCGCACTTTCAGCGGCCAAGCTggactc                                                                                              >  2:4587823/1‑145 (MQ=255)
                                                          gCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCg                                                                                             <  1:2350421/150‑1 (MQ=255)
                                                           cACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGGGCGTCGCGCTAGGGGCGCCACGCGGCCCCCGCCAATGCAGCGGCCAAGCAGCCATAGt                                                                                            >  1:3399772/1‑150 (MQ=255)
                                                            aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCCTTCGGCGCCATAGGGCCGCCAAACCGCCCCCCCCAAGTAAGCGGCCAAGCGGCaaaagtt                                                                                           >  2:1691966/1‑145 (MQ=255)
                                                            aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTa                                                                                           <  1:2842288/150‑1 (MQ=255)
                                                            aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTa                                                                                           <  2:4162993/150‑1 (MQ=255)
                                                            aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTa                                                                                           <  1:2066898/150‑1 (MQ=255)
                                                            aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCCAAACGGCCCCACCAAAATAAACCGCAAAGATGCaaaagta                                                                                           >  2:1893360/1‑145 (MQ=37)
                                                             cTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAg                                                                                          <  1:1481015/150‑1 (MQ=255)
                                                              ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCACTCTCAGCGGCAAAGCTGCAATCGTAgg                                                                                         >  2:6848219/1‑150 (MQ=255)
                                                              ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGTGCGGCGCGATAGGGCCGCAAAGCGGCCCCCGCAAATACAGCGGCAAAGCAGCCAAAGTAgg                                                                                         >  1:4462883/1‑150 (MQ=38)
                                                              ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCCAAGCAGCCATCGTAgg                                                                                         >  2:5806915/1‑150 (MQ=255)
                                                              ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAgg                                                                                         <  1:4816265/150‑1 (MQ=255)
                                                              gtGTCTTTACAAGCGGAAGGCATTGCTGGTGGCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCCCCCAACGGCCGCAAACCGGCCCCCCCAATTTCAGCGGCAAACCGGCAAGCGTAgg                                                                                         >  2:5425472/2‑150 (MQ=255)
                                                                gTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAgggg                                                                                       <  2:3244821/150‑1 (MQ=255)
                                                                   tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGgcgc                                                                                    <  1:2347932/150‑1 (MQ=255)
                                                                   tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGgcgc                                                                                    <  2:2851223/150‑1 (MQ=255)
                                                                     tACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTc                                                                                  <  2:742410/150‑1 (MQ=255)
                                                                     tACAAGCGGAAGGCAGTGCGGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTc                                                                                  <  2:1331872/150‑1 (MQ=255)
                                                                     gctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTc                                                                                  <  1:6445279‑M1/83‑1 (MQ=255)
                                                                      aCAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGTGGGTCGGGATAGGGGCGCAAAGCGGCCCCAGCAAAGTGAGCGGCCAAGCAGCAAACGGAGGGGCGCGcc                                                                                 >  1:6623905/1‑150 (MQ=255)
                                                                      aCAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTcgc                                                                                                                                               <  1:16579/88‑1 (MQ=255)
                                                                      aCAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTcgc                                                                                                                                               >  2:16484/1‑88 (MQ=255)
                                                                      aCAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCCGCCCCCGCAAATTAAGCGGCAAAGCTGCCAACGGAGGGGCGCGcc                                                                                 >  2:3466395/1‑150 (MQ=37)
                                                                        aaGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGc                                                                               <  1:3219230/150‑1 (MQ=255)
                                                                        aaGCGGAAGGCATTGCTGGTGTCGATCGGTGGAAGGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGGAGGGGCGCTCCGc                                                                               >  2:2091781/1‑150 (MQ=255)
                                                                         ccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCa                                                                              >  1:4518094‑M1/64‑150 (MQ=255)
                                                                         aGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCCCCCTAGGGCCGCAAACCCGCCCCACCCAATTCACCCGCAAACCAGCAAACCGAGGGGCGCGccccc                                                                              >  2:2294281/1‑147 (MQ=255)
                                                                            ggAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCAccc                                                                           <  1:5936722/150‑1 (MQ=255)
                                                                            ggAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCAccc                                                                           <  1:3782727/150‑1 (MQ=255)
                                                                               aGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCttt                                                                        <  2:2396743/150‑1 (MQ=255)
                                                                                 gCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTc                                                                      <  2:197106/150‑1 (MQ=255)
                                                                                 gCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTc                                                                      <  2:3358840/150‑1 (MQ=255)
                                                                                 gCATTGCGGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTc                                                                      <  2:5695926/150‑1 (MQ=255)
                                                                                  cATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTcg                                                                     <  1:620508/150‑1 (MQ=255)
                                                                                  agtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTcg                                                                     >  1:95496‑M1/55‑150 (MQ=255)
                                                                                   aTTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCCCAGCCGCCCCAGCAAATTAAACCGCCAAGCTGCAAAAGTAGGGGCGCGCCCCACCCGTTTcgc                                                                    >  2:3160178/1‑150 (MQ=37)
                                                                                   aTTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAAATCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCAcccaaatcgc                                                                    >  2:4394605/1‑143 (MQ=255)
                                                                                       cTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGc                                                                                                                      <  1:4086897/96‑1 (MQ=255)
                                                                                          gtgtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTAATAGGCAAAGCTGCAATAGTAGGGGCGCTCCGCACCCTTTTAGCGACGCaa                                                             >  2:6825982/1‑150 (MQ=38)
                                                                                          gtgtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATATCAGCCGCCAAGCTGCAAACGAAGGGGCGCTCCGCGCCCATTTCGCGCCGCac                                                             >  2:6223292/1‑149 (MQ=38)
                                                                                          gccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCaa                                                             >  2:6488992‑M1/47‑150 (MQ=255)
                                                                                           tgtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAg                                                            <  1:5234599/150‑1 (MQ=255)
                                                                                            caacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAgg                                                           >  2:5119656‑M1/45‑150 (MQ=255)
                                                                                             tCGATCGGTGGAATGTAGCGGCCTGGGCGCCGGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGc                                                          <  1:1683983/150‑1 (MQ=255)
                                                                                                  cGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCCGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTc                                                     <  2:5158777/150‑1 (MQ=255)
                                                                                                    gTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGCTATGGCCCCAAAGCGGCCCCACCCAATTCAGCCGCCCAGCTGCACTCGGAGGGGCGCGCCGCACCCCTTTCGCCACGCAAGGCCGCTccc                                                   >  2:4326553/1‑150 (MQ=255)
                                                                                                         aTGTAGGGGCCTGGGCGCCTGGGGGGCGGCGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGAAACGGGGCCCCACCAATTTCAGCGCCAACGCTGCAATGGTCGGGGGGCTCCGCACCATTTTGGCGACGCACGGCCGCTCCCACaaa                                              <  2:1227496/150‑1 (MQ=25)
                                                                                                         aTGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGTGCGGCCTTGGGTCGCGATTGGGCCGCAAAGGGGCCCCTGCAATTTCAGCGGCAATGCTGCAATCGTTGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACaaa                                              <  2:4986994/150‑1 (MQ=37)
                                                                                                         aTGTAGCGGCCTGGGCGCCTGTGGGGCGCCCGGTGGGAGCGGCCTGGCGTGGCGATAGGGCCGAAAAGCGGCCCCAGCAATTTACGCGGCAAAGGTGCAATCTTAGGGGCGCTCCGCACCGTTTTCGCGAGTAAAGGCCGCTCCCACaaa                                              <  2:5065141/150‑1 (MQ=255)
                                                                                                            tAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGTGCCTCGCGCAAGGGCCGCCAAGCCGCCCCAGCAAATTAAGCGGCCAAGCAGGAAACGTAGGGGGGCGCCGCCCCCCTGTCGCGACGCAAGGCCGCGCCCCCAAACTc                                           >  2:3460043/1‑150 (MQ=19)
                                                                                                                 gCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGGGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCGAAGCTGCAATCGTTGGGGGGCTGCGCCCCCTTTTGGCGTCGCCAGGCCGCTCCCCCAAAGTCCAAcc                                      >  2:5040721/1‑150 (MQ=25)
                                                                                                                    tGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCATGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAATTCATCGGCAAACCCTCAATCGTAGGGGCGCGCCGCACCCGGGTGCCGCCGCAACGCCGCTCCCCCAAAGCCCAACcccg                                   >  2:3590439/1‑150 (MQ=25)
                                                                                                                           cTGTGGGGCGGCCTGTGGGTGCGGCCTTGCGTCGCGATAGGGCCGAAAAGCGGCCCCAGCAATTTCAGCGGCAATGCTGCAATCGTCGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCcaca                            <  2:2819149/150‑1 (MQ=255)
                                                                                                                           cTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCCAAGCTGCAAACGGAGGGGCGCGCCGCACCCATTTCGCGACGCAAGGCCGCTCCCCCCAAGTCCAAACGCGCCCcacc                            >  2:5502364/1‑149 (MQ=25)
                                                                                                                            tgtgGGGCGGCCTGTGGGAGCGGGCCTTCGTCGCGATAGGGCCCCAAACCCGCCCCCCCAAATTCAGCCGCAAAAAAGCAAAAGTAGGGGCGCGCCCCACCCATTTAGCCACCCAACGCCCCACCCAAAAAAAACAAACCCCCCCcacac                           >  2:2490186/1‑149 (MQ=255)
                                                                                                                                gggCGGCCTGTGGGAGCGGCCTTGCGTCGCGATGGGGCCGCTTTGCGGcccta                                                                                                                        <  2:71187/53‑3 (MQ=17)
                                                                                                                                gggCGGCCTGTGGGAGCGGCCTGGCGTCGCGATAGGGCCGCACCGCGGCCCCCGCAATTTCAGCGGCACAGCTGCAAGCGTCGGGGCGCCCCGCACCCTTTCCGCGACGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACaggagg                       <  2:3814325/150‑1 (MQ=25)
                                                                                                                                gggCGGCCTGTGGGAGCGGCCGTGCGTCGCGCTAGGGCCGCAAAGCGGCCCCCCCCATTTCAGCGGCAAAGCTGCAATCGTAGGGGGGCTCCCCACCCATTTAGCGCCGCAAGGCCGCTCCCACCAACTCCAACCGCGCCCcacacgcgg                       >  2:5991930/1‑145 (MQ=25)
                                                                                                                                gggCGGCCGGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAACGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACaggagg                       <  2:3814066/150‑1 (MQ=255)
                                                                                                                                          gGGGCGCGGCCTGGCGGCGCGATCGGGCCGCAAAGCGGCCCCAGCAAGTTCCGCGGAAAAGCTGCAATCGTCGGGGCGCGCCGCACCATTTCCGCGCCGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGtt             >  2:3730366‑M1/2‑150 (MQ=255)
                                                                                                                                                     TTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGc  <  1:1862572‑M1/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGC  >  NC_002947/4084271‑4084569

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: