Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,926 A→T 6.7% intergenic (‑113/+86) 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,9260AT 6.7% 250.9 / 7.9 75intergenic (‑113/+86)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 = 3.15e-03

AATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAACGAATCTGCTGATCAAA  >  NC_000913/3959656‑3960188
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                      
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                                                                                                                                                                                                                                                                    <  6:870585/141‑1 (MQ=255)
                                                                                                                                      cAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCg                                                                                                                                                                                                                                                                    >  5:870585/1‑141 (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                                                                                                                                                                                                                                                                >  6:874256/1‑100 (MQ=255)
                                                                                                                                                                                   ttGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATTTTCCGGGCt                                                                                                                                                                                                                                                                <  5:874256/100‑1 (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                                                                                                                                                                                                             >  6:733569/1‑107 (MQ=255)
                                                                                                                                                                                                                               ttGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                                                                                                                                                                                                             <  5:733569/107‑1 (MQ=255)
                                                                                                                                                                                                                                  tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggg                                                                                                                                                                                                                                             <  6:698471/72‑1 (MQ=255)
                                                                                                                                                                                                                                  tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggg                                                                                                                                                                                                                                             <  8:321671/72‑1 (MQ=255)
                                                                                                                                                                                                                                  tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggg                                                                                                                                                                                                                                             >  7:321671/1‑72 (MQ=255)
                                                                                                                                                                                                                                  tcagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgggg                                                                                                                                                                                                                                             >  5:698471/1‑72 (MQ=255)
                                                                                                                                                                                                                                      tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCg                                                                                                                                                                                                                                                                    >  8:482859/1‑45 (MQ=255)
                                                                                                                                                                                                                                      tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCg                                                                                                                                                                                                                                                                    <  7:482859/45‑1 (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)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                      
AATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCAACTCATGTGGCAAACATTGATTAAGATACCCTCTGTAAGTCCTTGCTCTGTTTAATCTGTTTCGACTGACTTGTTGTACCAGGTCGGCTTTGGCCCAACGAATCTGCTGATCAAA  >  NC_000913/3959656‑3960188

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