Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,924 A→T 6.1% 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.1% 350.8 / 9.5 99intergenic (‑111/+88)ppiC/yifNpeptidyl‑prolyl cis‑trans isomerase C (rotamase C)/PemK toxin family pseudogene
Reads supporting (aligned to +/- strand):  ref base A (46/47);  new base T (2/4);  total (48/51)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

AAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGT  >  NC_000913/3959781‑3960064
                                                                                                                                               |                                                                                                                                            
aaGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATcc                                                                                                                                         >  4:1137595/1‑149 (MQ=255)
  ggATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCgg                                                                                                                                       <  4:978596/149‑1 (MQ=255)
  ggATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCgg                                                                                                                                       <  2:1445032/149‑1 (MQ=255)
    atatGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATTTTCCGGGc                                                                                                                                     <  1:967726/149‑1 (MQ=255)
     tatGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATTTTCCGGGCt                                                                                                                                    <  3:360513/149‑1 (MQ=255)
      atGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCtt                                                                                                                                   >  1:1164291/1‑149 (MQ=255)
        gCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATTTTCCGGGCtttt                                                                                                                                 <  4:70910/149‑3 (MQ=255)
        gCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCg                                                                                                                                 >  2:380076/1‑149 (MQ=255)
          aGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCa                                                                                                                               >  3:514407/1‑149 (MQ=255)
                  tgctgTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTc                                                                                                                       >  2:1438342/1‑149 (MQ=255)
                  tgctgTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTc                                                                                                                       <  3:1137595/149‑1 (MQ=255)
                     tgTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAg                                                                                                                    <  1:380076/149‑1 (MQ=255)
                      gTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAg                                                                                                                    <  4:817933/148‑1 (MQ=255)
                      gTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAg                                                                                                                    >  3:817933/1‑148 (MQ=255)
                       tttttGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                                                                                  >  3:677020/1‑149 (MQ=255)
                         tttGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGga                                                                                                                <  4:775187/149‑1 (MQ=255)
                             ccATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaa                                                                                                            >  2:1161445/1‑149 (MQ=255)
                              cATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg                                                                                                                   <  3:372245/141‑1 (MQ=255)
                              cATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg                                                                                                                   >  4:372245/1‑141 (MQ=255)
                              cATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaaa                                                                                                           <  4:514407/149‑1 (MQ=255)
                                  aTCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                                                       >  4:1446516/1‑149 (MQ=255)
                                  aTCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                                                       <  2:1164291/149‑1 (MQ=255)
                                   tCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCgg                                                                                                                                       <  2:91445/116‑1 (MQ=255)
                                   tCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCgg                                                                                                                                       >  1:91445/1‑116 (MQ=255)
                                   tCGTGCCTTGAGTGAATAAGTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                                                                      <  1:1161445/149‑1 (MQ=255)
                                        ccTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGATTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTc                                                                                                                                  <  1:875488/116‑1 (MQ=255)
                                        ccTTGAGTGAATAATTGTCAGGCGGCGATTTTAGAATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTc                                                                                                                                  >  2:875488/1‑116 (MQ=255)
                                              gTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATgg                                                                                           <  4:677020/149‑1 (MQ=255)
                                               tGAATAATTGTCAGGCGGCGATTTTAGTATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                                                                                  >  3:975996/1‑125 (MQ=255)
                                               tGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                                                                                  <  4:975996/125‑1 (MQ=255)
                                                  ataatTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                                                       <  3:697641/133‑1 (MQ=255)
                                                  ataatTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                                                       >  4:697641/1‑133 (MQ=255)
                                                   taatTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTGACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGga                                                                                                                >  1:484429/1‑123 (MQ=255)
                                                   taatTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTGACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGga                                                                                                                <  2:484429/123‑1 (MQ=255)
                                                   taatTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATTTTCCGGGCttt                                                                                                                                  <  4:507926/105‑2 (MQ=255)
                                                   taatTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATTTTCCGGGCttt                                                                                                                                  >  3:507926/1‑104 (MQ=255)
                                                   taatTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                                                                                     <  1:804748/134‑1 (MQ=255)
                                                   taatTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                                                                                     >  2:804748/1‑134 (MQ=255)
                                                      ttGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                                                                                  >  3:1116383/1‑118 (MQ=255)
                                                      ttGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                                                                                  <  4:1116383/118‑1 (MQ=255)
                                                         tCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                                                                                <  3:1446516/149‑1 (MQ=255)
                                                               ggcgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCgg                                                                                                                                       >  4:236631/1‑88 (MQ=255)
                                                               ggcgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCgg                                                                                                                                       <  3:236631/88‑1 (MQ=255)
                                                               ggcgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                                                                          >  1:599668/1‑149 (MQ=255)
                                                               ggcgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                                                                          <  1:1438342/149‑1 (MQ=255)
                                                                 cgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgg                                                                        <  2:599668/149‑1 (MQ=255)
                                                                         gCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtt                                                                >  1:702824/1‑149 (MQ=255)
                                                                          cATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCTCTCAAGGGGAGAAAACTTa                                                                                                       <  3:864698/109‑1 (MQ=255)
                                                                          cATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                                                       >  4:864698/1‑109 (MQ=255)
                                                                           aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccc                                                                                                                         <  2:772831/90‑1 (MQ=255)
                                                                           aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccc                                                                                                                         >  1:772831/1‑90 (MQ=255)
                                                                              aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATgg                                                                                           <  3:571202/117‑1 (MQ=255)
                                                                              aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATgg                                                                                           >  4:571202/1‑117 (MQ=255)
                                                                              aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttt                                                           >  3:1135317/1‑149 (MQ=255)
                                                                               tCCGGAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttta                                                          >  2:497960/1‑149 (MQ=255)
                                                                                  gCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                                                                                  >  4:529543/1‑106 (MQ=255)
                                                                                  gCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                                                                                  <  3:529543/106‑1 (MQ=255)
                                                                                        gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAAtt                                                 <  4:1135317/149‑1 (MQ=255)
                                                                                         cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAt                                                                    <  1:1249789/129‑1 (MQ=255)
                                                                                         cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAt                                                                    >  2:1249789/1‑129 (MQ=255)
                                                                                         cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAAttt                                                <  2:702824/149‑1 (MQ=255)
                                                                                          tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                                                                                <  3:208930/116‑1 (MQ=255)
                                                                                          tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                                                                                >  4:208930/1‑116 (MQ=255)
                                                                                          tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                                                                                <  1:681544/116‑1 (MQ=255)
                                                                                          tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTc                                                                                >  2:681544/1‑116 (MQ=255)
                                                                                           aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTAt                                              <  1:497960/149‑1 (MQ=255)
                                                                                            aCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCgg                                                                                                                                       <  3:1127432/59‑1 (MQ=255)
                                                                                            aCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCgg                                                                                                                                       >  4:1127432/1‑59 (MQ=255)
                                                                                            aCGTAATTGTCAGCCAGGGAGTGGCTGTTAGCTGTATGGCAAAAGCCCGGAAAACCCGGGCATCGCATGTCCCTCAAGGGGAGAAAACTATGGGCGCTTAGCGGCTAAAAAGCCATTAACAGGAGTGGTTTTTTTATTGTCAATTTCTg                                             >  2:1043953/1‑149 (MQ=255)
                                                                                                  ttGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                                                                       >  4:13113/1‑117 (MQ=255)
                                                                                                  ttGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                                                                       <  3:13113/117‑1 (MQ=255)
                                                                                                        gtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTATGGCCCCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTAt                                              >  3:1383635/1‑136 (MQ=255)
                                                                                                        gtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTAt                                              <  4:1383635/136‑1 (MQ=255)
                                                                                                           agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTAt                                              >  2:957191/1‑133 (MQ=255)
                                                                                                           agGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTAt                                              <  1:957191/133‑1 (MQ=255)
                                                                                                            gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAg                                                                                                                    <  4:9979/62‑1 (MQ=255)
                                                                                                            gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAg                                                                                                                    >  3:9979/1‑62 (MQ=255)
                                                                                                            gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                                                                      >  2:719086/1‑76 (MQ=255)
                                                                                                            gggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                                                                      <  1:719086/76‑1 (MQ=255)
                                                                                                             ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACtt                                                                                                        <  2:1166537/73‑1 (MQ=255)
                                                                                                             ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACtt                                                                                                        >  1:1166537/1‑73 (MQ=255)
                                                                                                                    cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                                                       <  2:19021/67‑1 (MQ=255)
                                                                                                                    cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                                                       >  1:19021/1‑67 (MQ=255)
                                                                                                                             atatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                                                       <  3:171916/58‑1 (MQ=255)
                                                                                                                             atatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                                                       >  4:171916/1‑58 (MQ=255)
                                                                                                                                tGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATgg                                                                  >  3:1427520/1‑92 (MQ=255)
                                                                                                                                tGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATgg                                                                  <  4:1427520/92‑1 (MQ=255)
                                                                                                                                tGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCctc         >  4:395814/1‑149 (MQ=255)
                                                                                                                                    aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                                                                     <  4:428805/85‑1 (MQ=255)
                                                                                                                                    aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGa                                                                     >  3:428805/1‑85 (MQ=255)
                                                                                                                                    aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGa                                            >  1:706037/1‑110 (MQ=255)
                                                                                                                                    aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGa                                            <  2:706037/110‑1 (MQ=255)
                                                                                                                                     aaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                                                                                  >  1:1387415/1‑39 (MQ=255)
                                                                                                                                     aaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                                                                                  <  2:1387415/39‑1 (MQ=255)
                                                                                                                                       aGCCCGGATTTTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGt  >  1:525263/1‑149 (MQ=255)
                                                                                                                                        gCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGc                                  <  2:1003101/116‑1 (MQ=255)
                                                                                                                                        gCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGc                                  >  1:1003101/1‑116 (MQ=255)
                                                                                                                                             gAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCattgat                            <  3:158096/117‑1 (MQ=255)
                                                                                                                                             gAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCattgat                            >  4:158096/1‑117 (MQ=255)
                                                                                                                                               |                                                                                                                                            
AAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGT  >  NC_000913/3959781‑3960064

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