Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,796,792 A→G 12.7% D272G (GAT→GGT)  waaC → ADP‑heptose:LPS heptosyl transferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,796,7920AG12.7% 187.5 / 3.0 72D272G (GAT→GGT) waaCADP‑heptose:LPS heptosyl transferase I
Reads supporting (aligned to +/- strand):  ref base A (10/52);  new base G (8/2);  total (18/54)
Fisher's exact test for biased strand distribution p-value = 1.22e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTTAAACTTCATTCATTGAAACCTTACACTCTGAAATCATCAATGATTTTAGAGATAAT  >  NC_000913/3796554‑3797017
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    
cccTGGCCGGAAGAACACTGGCGAGAATTGAGTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGa                                                                                                                                                                                                                                      <  2:201291/237‑1 (MQ=255)
 ccttgCCGGAAGACCACTGGCGAGAATTGTTTGGGTTACTGGCTGATTCAGGAATGCGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGAt                                                                                                                                                                                                                                     <  1:359761/235‑1 (MQ=255)
 acTGGCCGGAAGAACCCTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGAt                                                                                                                                                                                                                                     <  2:985406/239‑1 (MQ=255)
  cTGGCCGGAAGAACGCTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTTGGGCGCGCCGGATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATa                                                                                                                                                                                                                                    <  2:966926/239‑1 (MQ=255)
       cGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGAcc                                                                                                                                                                                                                               <  2:872352/238‑1 (MQ=255)
              acacTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATc                                                                                                                                                                                                                        <  2:869811/238‑1 (MQ=255)
               cccTGGCGAGAATGTATTGGTTTACTGGCTGATTCAGGAATACGGATTAACCTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCAc                                                                                                                                                                                                                      <  2:628794/237‑1 (MQ=255)
               cacTGGGGGTAATTTGTTGGTTTTCTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCAGGAGGAAGAACGGGCGAAACGACGGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCAc                                                                                                                                                                                                                      <  1:161556/239‑1 (MQ=255)
               cacTGGCGAGAATTGATTGGTTTTCTGGCTTATTCCGGAATACGGATTAACCTTCCGTGGGGCCCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCAc                                                                                                                                                                                                                      <  2:924480/239‑1 (MQ=255)
               cacTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGGTTTGCTTATGTTGAAGTATTGCCGAAGATGAGGCTGGAAGGCGTTGGCCGCGTGTTGGCCGGGGCTAAATTTGTAGGGTCGGGGGGTACGGGGGTAAGCCATTTTAAGGCGGGCCCGGGTGGACCCAAAATCAc                                                                                                                                                                                                                      >  1:59840/1‑239 (MQ=255)
               actggggggaaattttttgGTTTGCTGCGTGATTCAGGGAAACGGATTAAATTCCCGTGGGGGGGGCCGCATGAGGAAGTATGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTGAAAGTATTGCCGAAGATGAGTCTGCAAGGCGTTGCCCGCGCCCTGGCCGGGGCTAAAGTTGTAGTGGCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACGCAATATCAc                                                                                                                                                                                                                      <  1:528298/223‑1 (MQ=255)
                    gCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGGGGGTACGGGGTTAAGCCATTTAACGGCGGCCATGGGAAGAACCCAAAATCAGGGtt                                                                                                                                                                                                                  >  2:334485/1‑239 (MQ=255)
                              attgGTTTACTGGGGTGATCAGGGATCCGGGTTAAACTTGCGGGGGGCCGGCCGTATGGGGAATAAGGGGCGAAACGGCTGGCTGAAGGATTTTCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAAGTGGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAAc                                                                                                                                                                                                      <  1:120620/239‑1 (MQ=255)
                              attgGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAAc                                                                                                                                                                                                      <  2:59840/239‑1 (MQ=255)
                                   tttACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGGTACGGGGTTAAGCCATTTAACGGCGGCCACGGGTGAAACCATAATCACGGGTTATGGACCaaaccat                                                                                                                                                                                                   >  2:1019734/1‑233 (MQ=255)
                                                ttCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGAccc                                                                                                                                                                                                                              <  2:854997/199‑1 (MQ=255)
                                                  cgggAATCGGGATTAACCTTCCGGGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATtggt                                                                                                                                                                                   <  1:440552/237‑1 (MQ=255)
                                                  caggatttgggtttAACCTCCCGTGGGGGGCGCCGAAGGAGGAAGACCTGGCGAAACGGATGGCGGAAGGATTTGCTTATGTGGAAGTATTGCCGAAGATGAGTCGGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATtggt                                                                                                                                                                                   <  1:334485/227‑1 (MQ=255)
                                                  cAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGGTAGACCCAATATTCCGGTGTATTGACCAAACGGTCCGGGGTTAATttgt                                                                                                                                                                                   >  2:361466/1‑239 (MQ=255)
                                                             aTTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGGTAGACCCAAAATTAAGGGTTTTGGACCAACCGGACCGGGGTTAATTTGTggggtttggga                                                                                                                                                                        >  2:424765/1‑231 (MQ=255)
                                                             aTTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGAt                                                                                                                                                                                                                                     <  2:514101/180‑1 (MQ=255)
                                                                 aaCTTCCGGGGGGCGCGGCGCAGGAGGAAGAACGGGCGAAACGACGGGGGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAAt                                                                                                                                                                   <  1:424765/239‑1 (MQ=255)
                                                                  aCTTGCGTGGGGCGGGCCGAATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGGCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATGTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAAGATCACGGTTTATGGACCAACCGATCCTGGATTAATTGGTGGGTATGGGAAGAATc                                                                                                                                                                  <  2:847712/239‑1 (MQ=255)
                                                                       cGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACGGGATAGAACCCAAATCACGGGTTTTTGGCACACCGGTACTGGATTTATTGGTAGGTATTGGgaagataaaatg                                                                                                                                                              >  2:379993/1‑232 (MQ=255)
                                                                       cGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGGGTCGGTGGGTACGGGGGTAAGCCATTTAACGGCGGCCCCGGGTTGGACCAAAATAACGGTTTTTGGGCCAACCCGACCCGGGTTAATTGGGGGGGCTGGGAAGAATACGATg                                                                                                                                                             >  2:822106/1‑239 (MQ=255)
                                                                                cgccgcATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGGACTTGATAGACCCAATATACCGGATTATGGACCAACCCGTCCGGGGGTAATTGTGGGGGATGGGGAGAAtaggaagggattgagg                                                                                                                                                    >  2:621583/1‑224 (MQ=255)
                                                                                cgccgcATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCCCTGGGTAGACCCCAAATCAAGGGTTTTGGGCACAACGATCCGGGGGTCATTTGTGGGTATGGGGAAAAAAAAAGGGGATGTAgg                                                                                                                                                    >  1:467078/1‑239 (MQ=255)
                                                                                ccccgtgtggggAAGAACGGGCGAAACGGTTGGCGGAAGGATTTTGTTAGTTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAg                                                                                                                                                     <  1:361466/228‑1 (MQ=255)
                                                                                            aaaaaCGGGCGAAACGACTGGGGGGTGGATTTGCTTATGGTGAAGGTTTGACGAAGATGGGGTTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAgaga                                                                                                                                         <  1:858215/235‑1 (MQ=255)
                                                                                                aCGTGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCGTAAGATGAGTCTTGAAGGCGTTGCCCGCGTGCTGCCCGGGCCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGaaaa                                                                                                                                      <  2:334501/237‑1 (MQ=255)
                                                                                                aCGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGaaaa                                                                                                                                      <  2:467078/237‑1 (MQ=255)
                                                                                                    gCGAAACGACTGGCTGAAGGATTTTCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATttaa                                                                                                                                 <  2:292793/238‑1 (MQ=255)
                                                                                                            aCTGGCGGAAGGATTTGCTTATGTTGAAGTAGTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAg                                                                                                                                                                <  2:213377/199‑1 (MQ=255)
                                                                                                                     aGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAg                                                                                                               <  2:680011/239‑1 (MQ=255)
                                                                                                                                       aGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGGCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTat                                                                                             <  2:395532/239‑1 (MQ=255)
                                                                                                                                         tATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATttaa                                                                                                                                 <  2:610091/201‑1 (MQ=255)
                                                                                                                                          tttGCCGAAGTTGGGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGTGGCTAAATTTGTAGTGGCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGGCCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTatcat                                                                                          <  1:717371/238‑1 (MQ=255)
                                                                                                                                           ttGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTa                                                                                                                                                                                                                <  2:159689/122‑1 (MQ=255)
                                                                                                                                           ttGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTa                                                                                                                                                                                                                >  1:159689/1‑122 (MQ=255)
                                                                                                                                           ttGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGGTTATGGACCAACCCATCCGGGATTAATTGGGGGGGATGGGAAGAATCAGATGGAATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGGTatcatc                                                                                         >  2:321242/1‑239 (MQ=255)
                                                                                                                                           ttGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATAACGGTTTATTGACCAACCGATCCGGGGTTAATTGGTGGGGATTGGAAGAATCAGAAGGTATGAAGGGCTCCAAGAGAAAATTTAATTAAACTCAACAGTCAAGCAGTTTTGGAAAAGGTatcaac                                                                                         >  2:634096/1‑237 (MQ=255)
                                                                                                                                            tGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGTCCAACCGATCCGGGGTTAATTGGTGGGTATTGGAAGAATCAGATGGGATGTAGGGCTCCCAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGGTTTTGAAAAGTTAtcatca                                                                                        >  2:226460/1‑239 (MQ=255)
                                                                                                                                                gAAGAGGAGTCGGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATtat                                                                                    <  2:738094/239‑1 (MQ=255)
                                                                                                                                                   caggggtttggaggtggTTTCCCGCGTGTTGGCCGGGGCTAAATGTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCAGTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTACGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATtata                                                                                   <  1:405355/221‑1 (MQ=255)
                                                                                                                                                        gTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGGTAGAACCAATATCACGGTTTATGGACCAACCGCTCCGGGGTTAATTGGTGGGTATGGGAAGAATCAGATGGGATGTAGGGGTCCAAGGGAAAATTTAATTAAACCCAACAAGCAAACAAGTTTGGGAAAATTaaaaataatataaaagaaa                                                                            >  1:373716/1‑221 (MQ=255)
                                                                                                                                                              aaGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGGTTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCTTCCTTTTaaccgtccctttta                                                                      >  2:1011160/1‑227 (MQ=255)
                                                                                                                                                                 gCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGGAGAATCAGATGGTATGTAGGGCTCCAAGGGAAAATTTAATTAACCTCAAAAGTCAAGCAGTTTTGGAAAAGGTATCATCCTTATAAAGGTaaaaaaccgt                                                                    >  2:405355/1‑232 (MQ=255)
                                                                                                                                                                   gTTTCCCGCGTGCGGGCCTGGCCTAAATTTGTAGTGTCGGTGAATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTAGGGACCAACCGATCCGGGATTAATTGGTGGTTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAAc                                                                 <  1:379993/239‑1 (MQ=255)
                                                                                                                                                                   gTTGCCCGCGTTTTGGCCGGGGCTAAATTTGTGGTGTCTGGGGAAACGGGGTTAAGCCATTTAACGGCGGCACTGGGTAGACCCAATATCCCGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAAc                                                                 <  1:822106/239‑1 (MQ=255)
                                                                                                                                                                    tttccctggttcaggCCGGGAGTAATGTTGTGGTGTCGGTTGATACGGGGTTAAGCCATTTAACGGCGGCTCTGGATAGACCCAATATCGCGGTTTATGGACCATCCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACa                                                                <  1:634096/226‑1 (MQ=255)
                                                                                                                                                                    ttGCCCGCGTGCTGGCCGGGGTTAAATTGGTAGTGTCGGTGGATACGGGGTCAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAAc                                                                 <  2:1720/238‑1 (MQ=255)
                                                                                                                                                                            gTGCTGGCCGGGGCTAAATTTGTAGTGGCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTTa                                                         <  2:1002535/238‑1 (MQ=255)
                                                                                                                                                                                   ccGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTTAAACttca                                                  <  2:373716/238‑1 (MQ=255)
                                                                                                                                                                                    cGGGGCTAATTTTTGAGGGCCGGGGGATCAGGGGTTAAGCCATGTAACGGCGGCATTGGATGGACCTAATATCACGGTTTAGGGACCAACCGATCCGGGATTAATTGGTGTGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTTAAACttcat                                                 <  1:1019734/238‑1 (MQ=255)
                                                                                                                                                                                          ttatttttTAGTGGCGGTGGATTCGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTTAAACTTCATTCATTg                                           <  1:289300/234‑1 (MQ=255)
                                                                                                                                                                                          tAAATTTTTAGTTGCGGTGGATACGGGGTTAAGCCATTTAACGGAGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTTAAACTTCATTCATTg                                           <  1:621583/238‑1 (MQ=255)
                                                                                                                                                                                          tAAATTTGTAGGGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTTAAACTTCATTCATTg                                           <  2:249528/238‑1 (MQ=255)
                                                                                                                                                                                          gatatttgttggttgggTGGATTCGGGGTTAAGGCATTTAACGGTGGCGCTGAATAGACCCAATTTCCCGGTTTATGGAGCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTACAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAGCATGCTGACATCCTTTGAAGTACATTCGTTg                                           <  1:76509/223‑1 (MQ=255)
                                                                                                                                                                                           atttttttggtttCGTTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTTAAACTTCATTCATTGaa                                         <  1:321242/227‑1 (MQ=255)
                                                                                                                                                                                                 gttggttCGGTGGATTCGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTTAAACTTCATTCATTGAAACCTTa                                   <  1:226460/233‑1 (MQ=255)
                                                                                                                                                                                                 ggaGTTTCGGGTGATTCGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTTAAACTTCATTCATTGAAACCTTa                                   <  2:769638/237‑1 (MQ=255)
                                                                                                                                                                                                       tCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATttaa                                                                                                                                 >  2:7450/1‑139 (MQ=255)
                                                                                                                                                                                                       tCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATttaa                                                                                                                                 <  1:7450/139‑1 (MQ=255)
                                                                                                                                                                                                                gacgggTTAAGCCAATTAACGGCGGGACTGGATAGACCCAAGATCACGGTGTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTTAAACTTCATTCATTGAAACCTTACACTCTGAAatcatc                    <  2:440912/236‑1 (MQ=255)
                                                                                                                                                                                                                aCGGGGGTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGAGTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTTAAACTTCATTCATTGAATCCTTa                                   <  2:32525/224‑1 (MQ=255)
                                                                                                                                                                                                                          gCCATTTAACGTGGGCCGTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCGGGGATTAATTGGGGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTAACATCATTATAAAGGTAAAACATGCTAACATCCTTTAAACTTCATTCAGTGAAACCTTACACTCTGAAATCATCAATGATTTTa          <  1:309967/239‑1 (MQ=255)
                                                                                                                                                                                                                            caattAACGGCGGGAGTGGGTAGACCCAAGATCACGGTTTATGGACCAACCGATCCGGGTGTAATTGGTGGGTATGGGAATAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTGAAACTTCATTCATTGAAACCTTACACTCTGAAATCATCAATGATTTTa          <  1:869177/234‑1 (MQ=255)
                                                                                                                                                                                                                            cATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTAGGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTTAAACTTCATTCATTGAAACCTTACACTCTGAAATCATCAATGATTTTaga        <  2:8616/239‑1 (MQ=255)
                                                                                                                                                                                                                              tttAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATc                                                                                                                                                                                                 <  2:506849/53‑1 (MQ=255)
                                                                                                                                                                                                                              tttAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATc                                                                                                                                                                                                 >  1:506849/1‑53 (MQ=255)
                                                                                                                                                                                                                              tttAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAg                                                                                                               <  2:304860/134‑1 (MQ=255)
                                                                                                                                                                                                                                  acggggggccTGGATGGACCGAATATTACGGTTTATGGACCAAGCGATCCGGGATTTATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATGCTCATTATAAAGGTAAAACATGCTAACATCCTTTAAACTTCATTCATTGAAACCTTACACTCTGAAATCATCAATGATTTTAGAGataat  <  2:352003/230‑1 (MQ=255)
                                                                                                                                                                                                                                  acggggggcatGGATAGACCCAATATCACGGTTTAAGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTTAAACTTCATTCATTGAAACCTTACACTCTGAAATCATCAATGATTTTAGAGataat  <  2:91095/229‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    
CACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTTGCCCGCGTGCTGGCCGGGGCTAAATTTGTAGTGTCGGTGGATACGGGGTTAAGCCATTTAACGGCGGCACTGGATAGACCCAATATCACGGTTTATGGACCAACCGATCCGGGATTAATTGGTGGGTATGGGAAGAATCAGATGGTATGTAGGGCTCCAAGAGAAAATTTAATTAACCTCAACAGTCAAGCAGTTTTGGAAAAGTTATCATCATTATAAAGGTAAAACATGCTAACATCCTTTAAACTTCATTCATTGAAACCTTACACTCTGAAATCATCAATGATTTTAGAGATAAT  >  NC_000913/3796554‑3797017

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