Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,924 A→T 6.7% intergenic (‑111/+88) ppiC ← / ← yifN peptidyl‑prolyl cis‑trans isomerase C (rotamase C)/PemK toxin family pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,959,9240AT 6.7% 260.6 / 7.6 75intergenic (‑111/+88)ppiC/yifNpeptidyl‑prolyl cis‑trans isomerase C (rotamase C)/PemK toxin family pseudogene
Reads supporting (aligned to +/- strand):  ref base A (33/37);  new base T (2/3);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

GGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAACGAATCTGCTGATCAAA  >  NC_000913/3959651‑3960188
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                        
ggCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGACGATGCTAACGTAATTGTCAGCCAGGGAGTGGCTGGTACCATTATGGCAAAACCCCGGaaa                                                                                                                                                                                                                                                                         >  4:25952/1‑275 (MQ=255)
     aaTTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATcc                                                                                                                                                                                                                                                                     >  1:124020/1‑274 (MQ=255)
      aTTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGGGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCGGGGAGTGGCTGGTACCTAATGGGCAAAAGCCCGGAAAATCCg                                                                                                                                                                                                                                                                    >  4:286307/1‑274 (MQ=255)
       ttCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCggg                                                                                                                                                                                                                                                                  <  2:25953/275‑1 (MQ=255)
                  cACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGt                                                                                                                                                                                                                                                        <  2:301348/274‑1 (MQ=255)
                    ctgccttgTTTGCCGGATGGTCACTTNNNTTGTTTCTTCGCCAGCTTGCCGAATTCGTCCCCGTGTCTATGTTGCTCGAGAAGATCCAGGGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTTTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGt                                                                                                                                                                                                                                                        <  3:159515/265‑1 (MQ=255)
                                ccTTATGGGCAAATGGGGTGTTTCTTCGCCAACTTGCCGAAATCGGCCCCGTTCTTAATCTGCGCGAGAAGATTCAGAGCCAGTTTCTCTTCTTTTACAAGGTTATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaa                                                                                                                                                                                                                                         <  4:188168/275‑1 (MQ=255)
                                     tGGTCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACtt                                                                                                                                                                                                                                    <  2:286308/275‑1 (MQ=255)
                                                  ggtTGCTGCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCGTTGACAAGGATATGCAGTGCTGCTGCTGTTTTTGCTAGGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGTATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGACGGGAAAATCCGGGCTTCGCAGGTCCATCAAGGGGAGAAAACTTAGGGCCTCTATgg                                                                                                                                                                                                                       >  1:85376/2‑275 (MQ=255)
                                                       cttcGCCAGCTTTCCGGAATCGGCCCCGTTCTTAATCTGCTCGGGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTa                                                                                                                                                                                                                   <  4:216297/274‑1 (MQ=255)
                                                            ccAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                                                                                                                                                                                                             <  4:497504/275‑1 (MQ=255)
                                                                 tttccGAAAGCGGCCCCGTTCTTAATCTGCTCGAGAATATCCAGAGCCAGTTTCGCTTCTTTTACAAGGATATGCAGTGCTGCGGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGGCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAtt                                                                                                                                                                                                         <  1:127306/271‑1 (MQ=255)
                                                                        aaTCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                                                                                                                                                                                                 <  2:164494/275‑1 (MQ=255)
                                                                                                        tCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCCAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATCTATGGCAAAAGCCCGGTAAATCCGGGCTTCGCACGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACCACGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat                                                                                                                                                                 >  2:187268/1‑275 (MQ=255)
                                                                                                             agCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                                                                                                                                                                                                  >  4:8343/1‑205 (MQ=255)
                                                                                                             agCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                                                                                                                                                                                                  <  2:8343/205‑1 (MQ=255)
                                                                                                                 aGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGCACATCCTGGCTTCGCAGGTCCCTCAAGGGGATAAAACTTAGGGCCTCTATGGGATACAACGTCATTATCGGGGCTGGTTATTTTATTGTACTTTTCTGAGCGATATGGCattga                                                                                                                                                         >  4:136858/1‑274 (MQ=255)
                                                                                                                             tttACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGTTCGTGCCTTGAGTGAAAAATTGTCAGGCTGCGATGTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGAc                                                                                                                                            <  2:390727/275‑1 (MQ=255)
                                                                                                                                 caagaattaggtgtggtgTTGCTGTTTTGGCCGGGATCGTGCCTTGAGTTAAAAATTGTCAGGCGGGGATTTTAGCATGTTCCGGAGATGCTAACGTAATTGTCAGTCCGGGGGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGcc                                                                                                                                        <  4:187267/263‑1 (MQ=255)
                                                                                                                                  aaGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCg                                                                                                                                          <  2:290094/272‑1 (MQ=255)
                                                                                                                                   aGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGAGGCTGGTACATATATGGCAAAAGCCCGGAAAATcc                                                                                                                                                                                                                                                                     >  7:778613/1‑148 (MQ=255)
                                                                                                                                           cAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCg                                                                                                                                                                                                                                                                    >  5:870585/1‑141 (MQ=255)
                                                                                                                                           cAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCg                                                                                                                                                                                                                                                                    <  6:870585/141‑1 (MQ=255)
                                                                                                                                              tgctgctgcTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATTCCCCCCCtctgtctgtctg                                                                                                                           >  4:348058/1‑275 (MQ=255)
                                                                                                                                                ctgctgctgTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGGTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGATTTGATGCCCGCCCTCTGTATGTCTGca                                                                                                                         >  2:136083/1‑275 (MQ=255)
                                                                                                                                                       tgTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCaa                                                                                                                                                                                                                                                 <  8:778613/148‑1 (MQ=255)
                                                                                                                                                               ccaaggtCGTGCCTTAGGGTAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGGGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGc                                                                                                          <  4:136082/269‑1 (MQ=255)
                                                                                                                                                                    aTCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATcc                                                                                                                                                                                                                                                                     <  8:794279/115‑1 (MQ=255)
                                                                                                                                                                    aTCGTGCCGTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGACAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATcc                                                                                                                                                                                                                                                                     >  7:794279/1‑115 (MQ=255)
                                                                                                                                                                     tcttGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGCCTGGTACATATATGGCAAAAGCCCGGAAAATCCg                                                                                                                                                                                                                                                                    >  8:721962/4‑115 (MQ=255)
                                                                                                                                                                     tCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCg                                                                                                                                                                                                                                                                    <  7:721962/115‑1 (MQ=255)
                                                                                                                                                                            ttGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATTTTCCGGGc                                                                                                                                                                                                                                                                 <  8:1007217/111‑1 (MQ=255)
                                                                                                                                                                               aGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGAGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATTACCGCCCTCTGTCTGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGa                                                                                           >  2:65398/1‑274 (MQ=255)
                                                                                                                                                                               aGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCAGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGACTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGCTGGCCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGCTTAAGa                                                                                           >  2:397147/1‑274 (MQ=255)
                                                                                                                                                                                   acagatTGTCAGGCGGCGATTTTAGCATGATCCGAAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATAcc                                                                                       <  2:172113/270‑1 (MQ=255)
                                                                                                                                                                                   aataatTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATAcc                                                                                       <  2:453877/274‑1 (MQ=255)
                                                                                                                                                                                        ttGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATTTTCCGGGCt                                                                                                                                                                                                                                                                <  5:874256/100‑1 (MQ=255)
                                                                                                                                                                                        ttGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATTTTCCGGGCt                                                                                                                                                                                                                                                                >  6:874256/1‑100 (MQ=255)
                                                                                                                                                                                             aGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTc                                                                            >  2:187021/1‑275 (MQ=255)
                                                                                                                                                                                             aGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTc                                                                            >  2:341326/1‑275 (MQ=255)
                                                                                                                                                                                               gcggcgATCTTGGCAGGATCCGCAGTAGTTAACGTAATTGTCAGTCAGGGAGTGGCTGTTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTcc                                                                           <  2:252047/274‑1 (MQ=255)
                                                                                                                                                                                                  gcgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATTTGGCAAACCTTGATTAAgc                                                                                           >  4:172112/1‑254 (MQ=255)
                                                                                                                                                                                                       tttAGCAAGATCCGCAGAAGCTAAAGTAATTGTAAGTCAGGTAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGGCCCTCAAGGGGAGAAAACTTATGTCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGt                                                                  <  4:397146/275‑1 (MQ=255)
                                                                                                                                                                                                             aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                                                                                                                                                                                                 >  6:556166/1‑142 (MQ=255)
                                                                                                                                                                                                             aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                                                                                                                                                                                                 <  5:556166/142‑1 (MQ=255)
                                                                                                                                                                                                               gATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                                                                                                                                                                                   <  8:569799/106‑1 (MQ=255)
                                                                                                                                                                                                               gATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                                                                                                                                                                                   >  7:569799/1‑106 (MQ=255)
                                                                                                                                                                                                                   cGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATTTTCCGGGc                                                                                                                                                                                                                                                                 <  8:69338/72‑1 (MQ=255)
                                                                                                                                                                                                                     cAGATTCTAACGTAATTGTCAGGCAGGGGGTTGCTGGTACAGATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGGCCCGCAAGGGGAGAAAACTTAGGGCCTCTATGAGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCga                                                    <  4:65397/275‑1 (MQ=255)
                                                                                                                                                                                                                           cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCa                                                                                                                                                                                                                                                           >  7:579584/1‑70 (MQ=255)
                                                                                                                                                                                                                           cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCa                                                                                                                                                                                                                                                           <  8:579584/70‑1 (MQ=255)
                                                                                                                                                                                                                           cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGcaaccaac                                                                                                                   >  4:207195/1‑206 (MQ=255)
                                                                                                                                                                                                                               cGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCAt                                                                                                               <  2:207196/206‑1 (MQ=255)
                                                                                                                                                                                                                                    ttGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                                                                                                                                                                                                             <  5:733569/107‑1 (MQ=255)
                                                                                                                                                                                                                                    ttGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                                                                                                                                                                                                             >  6:733569/1‑107 (MQ=255)
                                                                                                                                                                                                                                       tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggg                                                                                                                                                                                                                                             <  8:321671/72‑1 (MQ=255)
                                                                                                                                                                                                                                       tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggg                                                                                                                                                                                                                                             >  7:321671/1‑72 (MQ=255)
                                                                                                                                                                                                                                       tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggg                                                                                                                                                                                                                                             <  6:698471/72‑1 (MQ=255)
                                                                                                                                                                                                                                       tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggg                                                                                                                                                                                                                                             >  5:698471/1‑72 (MQ=255)
                                                                                                                                                                                                                                           tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCg                                                                                                                                                                                                                                                                    <  7:482859/45‑1 (MQ=255)
                                                                                                                                                                                                                                           tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCg                                                                                                                                                                                                                                                                    >  8:482859/1‑45 (MQ=255)
                                                                                                                                                                                                                                                 aGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                                                                                                                                                                                                      >  4:555433/1‑101 (MQ=255)
                                                                                                                                                                                                                                                     gCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgggg                                                                                                                                                                                                  <  2:555434/101‑1 (MQ=255)
                                                                                                                                                                                                                                                     gCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTATGTAAGTACTTGCACTGTTTAATCTGTTTCgactga                                                >  1:12422/1‑247 (MQ=255)
                                                                                                                                                                                                                                                         gTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGCGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATTTGGCAAAACTTGATTAAGATACCCTATGTAAGTACTTGCTCTTTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAACGa                >  4:216558/1‑275 (MQ=255)
                                                                                                                                                                                                                                                           aCATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAACGAAt              >  2:89095/1‑275 (MQ=255)
                                                                                                                                                                                                                                                             atatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                                                                                                                                                                                                      >  1:84221/1‑89 (MQ=255)
                                                                                                                                                                                                                                                                 atGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgggg                                                                                                                                                                                                  <  3:84221/89‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aaGCCGGGAAAATCCGTGCTGCGCAGGGCCCTCAAGGGGAGAAAACTTAGGGCCTTTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGGTTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAACGAATCTGCTGATCaag  <  2:183991/274‑2 (MQ=255)
                                                                                                                                                                                                                                                                          gCCCGGATTTTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTtgatg                                                                                                                                              >  7:1096223/1‑132 (MQ=255)
                                                                                                                                                                                                                                                                             cGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                                                                                                                                                                                                      <  7:998918/73‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                                                                                                                                                                                                      >  8:998918/1‑73 (MQ=255)
                                                                                                                                                                                                                                                                             cGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATgg                                                                                                                                                                                              >  7:478989/1‑81 (MQ=255)
                                                                                                                                                                                                                                                                             cGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATgg                                                                                                                                                                                              <  8:478989/81‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATgg                                                                                                                                                                                              <  6:537739/81‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATgg                                                                                                                                                                                              >  5:537739/1‑81 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                        
GGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAACGAATCTGCTGATCAAA  >  NC_000913/3959651‑3960188

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