Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,796,150 G→A A58T (GCA→ACA)  waaC → ADP‑heptose:LPS heptosyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,796,1500GA98.0% 169.5 / NA 50A58T (GCA→ACA) waaCADP‑heptose:LPS heptosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/30);  minor base T (0/1);  total (19/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAG  >  NC_000913/3795909‑3796389
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
aaTCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGt                                                                                                                                                                                                                                         <  1:86939/250‑1 (MQ=255)
 aTCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCg                                                                                                                                                                                                                                       <  2:149181/251‑1 (MQ=255)
 aTCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCg                                                                                                                                                                                                                                       <  1:271741/251‑1 (MQ=255)
    gACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTg                                                                                                                                                                                                                                    <  2:49139/251‑1 (MQ=255)
      cATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGc                                                                                                                                                                                                                                  >  1:306292/1‑251 (MQ=255)
             ccccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATTCCCAGCAGCCAATCCCAGGGATTAAGTTTGGCTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATGCCTGTGACAATACGTCGCTGGCGTaaa                                                                                                                                                                                                                             <  1:295469/249‑1 (MQ=255)
                     cgTACTGGAAGAACTCGACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTtct                                                                                                                                                                                                                   <  2:87570/251‑1 (MQ=255)
                                aaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCAt                                                                                                                                                                                                         >  1:303345/1‑250 (MQ=255)
                                      aCGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCg                                                                                                                                                                                                  >  1:388916/1‑251 (MQ=255)
                                          gcTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAAc                                                                                                                                                                                              <  2:84196/251‑1 (MQ=255)
                                               tGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCaa                                                                                                                                                                                          >  1:5721/1‑250 (MQ=255)
                                                           aaGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCACAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGACCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGa                                                                                                                                                                              >  2:299054/1‑250 (MQ=255)
                                                            aGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGa                                                                                                                                                                              >  1:312666/1‑249 (MQ=255)
                                                                tGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCa                                                                                                                                                                                           >  2:5180/1‑232 (MQ=255)
                                                                tGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCa                                                                                                                                                                                           <  1:5180/232‑1 (MQ=255)
                                                                 gACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTCAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCt                                                                                                                                                                                                                                     >  2:238152/1‑189 (MQ=255)
                                                                 gACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGAGGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCt                                                                                                                                                                                                                                     <  1:238152/189‑1 (MQ=255)
                                                                       tGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGc                                                                                                                                                                  >  1:42243/1‑250 (MQ=255)
                                                                                  aTCGTTAAAACATCGTCGAGGGGCGAGCTTCTCCATCCGTTGCCCGCACTCCATGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATTCGTCACTGGCGTAAAGCCTGGATCTCGGCCCCCATAAACACGGAACGCAATGCGGTTCGTGAAGCGCTACACGCAGAGAACTATg                                                                                                                                                       >  2:222249/1‑250 (MQ=255)
                                                                                       tAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCa                                                                                                                                                  <  1:291086/250‑1 (MQ=255)
                                                                                         aaaCATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCTCTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATCCGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGtt                                                                                                                                               <  1:401146/251‑1 (MQ=255)
                                                                                                   gATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTc                                                                                                                                     <  1:37433/251‑1 (MQ=255)
                                                                                                      gggCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAgg                                                                                                                                  >  1:223476/1‑251 (MQ=255)
                                                                                                                           gCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGAGTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGTGCTGGTAAAAAGCgcggcg                                                                                                              <  2:159193/250‑1 (MQ=255)
                                                                                                                                    cactGATGCCCAGCAGGCAAGCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGAc                                                                                                      <  2:186870/249‑1 (MQ=255)
                                                                                                                                           gCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGgc                                                                                             <  2:96154/251‑1 (MQ=255)
                                                                                                                                             ccAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGCAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGgcgc                                                                                           >  1:329760/1‑251 (MQ=255)
                                                                                                                                               agcagGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCAt                                                                                         >  1:43853/1‑251 (MQ=255)
                                                                                                                                                gcagGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATg                                                                                        <  2:142298/251‑1 (MQ=255)
                                                                                                                                                  agGCAATCCCAGGGCTTAAGATTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTAAAAGCAGAGAACTATGACGCAGGTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCAGGAGACGCGTCTGGAGCATGGc                                                                                      >  1:104041/1‑251 (MQ=255)
                                                                                                                                                  agGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGc                                                                                      <  1:184304/251‑1 (MQ=255)
                                                                                                                                                                  tAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGcatggcatgg                                                                      >  2:205664/1‑251 (MQ=255)
                                                                                                                                                                               gtggtgGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAgcatggcat                                                                        >  1:219120/1‑236 (MQ=255)
                                                                                                                                                                               gtggtgGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAgcatggcat                                                                        <  2:219120/236‑1 (MQ=255)
                                                                                                                                                                                  gtggAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAGAGCATGGCATGGACTGGCAAACCGCTcg                                                      <  2:88161/251‑1 (MQ=255)
                                                                                                                                                                                       agaagGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAAc                                                 <  1:292134/251‑1 (MQ=255)
                                                                                                                                                                                       agaagGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAAc                                                 <  2:43853/251‑1 (MQ=255)
                                                                                                                                                                                          agGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTGGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACACAAAGCGTTTCGTGAAGCGCTGCAAGCAGAGAACTATGACGCGGTTATCGACTCTCAGTGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCTTGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCtt                                              <  1:69488/251‑1 (MQ=255)
                                                                                                                                                                                               tCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCt                                                                                                                              <  1:51462/166‑1 (MQ=255)
                                                                                                                                                                                               tCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCt                                                                                                                              >  2:51462/1‑166 (MQ=255)
                                                                                                                                                                                                 gccgtgATTCCTGCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGTCACTACGTCGCTGTCGTAAAGCTTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCTTTTCGTGAAGCGCTACAAGCAGATAACTATGACGCAGTTATCGACGCTCATGTGCTGGTAAAAATCGCGGCGCTGGTGACGCGGCTGGTGCATTGCGTAAAGCATGGCATGGACTGGCAACCCGATCGCGAACCTTTagccag                                       <  2:145656/246‑1 (MQ=255)
                                                                                                                                                                                                      gATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGt                                  <  2:312666/251‑1 (MQ=255)
                                                                                                                                                                                                      gATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGt                                  <  1:62497/251‑1 (MQ=255)
                                                                                                                                                                                                         tccttcgtGGCAAGCTGCCGCCGAACGCGCTATTCCTGTTACAAAACGCTGCTGGTGTAAAGCCTGGTTCTCGCCCCCGCTAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGCGGCGGGTAAAAAGGGCGGCGCTGGTGACGCATCTGACGCATGGAGTAAAGCATGGCATGGACTGGCAAACCGACCGCGAACCTTTAGCCAGCCTGtctt                               >  2:410310/1‑248 (MQ=255)
                                                                                                                                                                                                         tccttcctGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGttt                                <  1:7604/250‑1 (MQ=255)
                                                                                                                                                                                                                              tGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAATCGTTTCGTGAAGCGCTACAAGCAGAGAACTATCACGCAGTTATCGACGCTCAGGTGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTTTTTTACAATCGTAAGCATCATAtt           <  2:293874/250‑1 (MQ=255)
                                                                                                                                                                                                                              tGAGCGAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATAtt           <  1:299054/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gAGTTATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACATGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCaaa      >  1:363088/1‑250 (MQ=255)
                                                                                                                                                                                                                                      ttATTCCTGTGACAATTAGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTACGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCTCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGTGAGCCTGTTTTGCAATCGTAAGCTTCATATTGCaaatca   <  1:322563/250‑4 (MQ=255)
                                                                                                                                                                                                                                       tATTCCTGTGACAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAAcag  <  1:325422/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
AATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGTAAAGCCTGGTTCTCGGCCCCCATAAAAGCGGAACGCAAAGCGTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAG  >  NC_000913/3795909‑3796389

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: