Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,477,878 T→A 7.7% intergenic (+20/+29) slyX → / ← slyD phi X174 lysis protein/FKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,477,8780TA 7.7% 156.6 / 6.1 54intergenic (+20/+29)slyX/slyDphi X174 lysis protein/FKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base T (23/26);  new base A (2/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-04

TCCGCGGCGACCATTCGACAAAAGAGGTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCA  >  NC_000913/3477607‑3478145
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                            
tCCGCGGCGACCATTCGCCAAAAGAGGTGAATCATGCAGGATTTATCAGTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTCCGTGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATAAAAGGGAGATGGCGAAACTGCGCGAGCATCTGCCTCTGCTGACGGAGATGTTAAAAACCAGCTAGCCGTCGAACATCCCGACCGA‑TGTTCAAGACACGCCACCGCAACACACTTAAGGACTAAAAAAAACGGGGAtac                                                                                                                                                                                                                                                                            >  4:514551/1‑272 (MQ=255)
tCCGCGGCGACCATTCGACAAAAGAGGTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCAAAGTATTGAGGCGTAAAAAAAGCGGGGATTc                                                                                                                                                                                                                                                                            >  2:284720/1‑274 (MQ=255)
                                  tGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGGGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGgca                                                                                                                                                                                                                                        <  2:85904/275‑1 (MQ=255)
                                                 tGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTTAGGAGTAAAAAAAGCGGGGGTTCCACGCTTTTTTGTCACTTTTT‑CGGTATTAGTGTCAACCCCAACCGCCAtt                                                                                                                                                                                                                         >  1:105781/1‑275 (MQ=255)
                                                 tGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGtt                                                                                                                                                                                                                         >  2:366256/1‑275 (MQ=255)
                                                 tGGAAGCACGCCTGGCAGAGCTGGAGAGCCGAATGGCTTTTCAGGAAATCACCATTGAATAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGAACGAGAAGTNAAAAGCCAGCAACCCGTCGAACATCGACTAGCA‑GGCTGAAGAAACGCCACCGCCAACATAATGAGGAGTAAAAAAAGCGGGGATTCACCGATCTTTTGCCACTTTTN‑c                                                                                                                                                                                                                                                     >  3:110784/1‑247 (MQ=255)
                                                     aGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCCCGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTACAAGCCAGCCAGCCGTCGAACATCCCGTCGCA‑GGCCACATAAACCCCACCCCCACATTATTCAGGCCTACAAAACCCCGGGACTCCCCGCTTTTTTGGTCCTTTTT‑CTGTATTAGTGGCAACCCCACCCCCCGTTGccc                                                                                                                                                                                                                     >  4:516335/1‑274 (MQ=255)
                                                            tgtgaagagCGGGATGGCCGCTTGCCTTTTCAGGAAATCCCCCTTTAATAACTGAACGTCCCGGTGACCGGCCATGAAATGGAGCTGGGGAACCTGCTCCATCATCTGCGTCGGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGc                                                                                                                                                                                                              <  4:404958/270‑1 (MQ=255)
                                                               gCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCca                                                                                                                                                                                                            >  2:472254/1‑274 (MQ=255)
                                                                    gCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAACAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCCACCGCAAACCCCTTTGCCTTTTCCGcccccgccc                                                                                                                                                                                                      >  4:105008/1‑268 (MQ=255)
                                                                                aCTGGCTTTTCAGGAAATCACCATTGAAGAACTGCACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTCAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTATGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGGTCCCTTTACCGCCACAGCAGCCTTAGcacccg                                                                                                                                                                                          >  2:439645/1‑270 (MQ=255)
                                                                                 cgggtatttCAGTAATACGCCATTGGANNCATGAACGTCACGGTGCGTGATCATGAAAGNGAGAGGGCGAAACTGCGCGTTCATCTGCGTCTGCTGACCGAGAAGTAACAAGCAAGCCAACCGTCAACAATCGANTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCAcc                                                                                                                                                                                           <  3:57955/267‑1 (MQ=255)
                                                                                                 taaCCATTGAAGAACGGAACGTCACGGTGACCGCTCATGAAATGGAGAGGGCGAAACTGCGCGATCGTCTGCGTCGGCTGACCGGGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATc                                                                                                                                                                         <  2:4103/273‑1 (MQ=255)
                                                                                                                    cGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCCGGACCGCCgcag                                                                                                                                                      >  1:9147/1‑275 (MQ=255)
                                                                                                                            tGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACcgtc                                                                                                                                              >  2:272381/1‑275 (MQ=255)
                                                                                                                             taCCGCTCATGAAAGGGAAATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACcgtcg                                                                                                                                             <  2:105009/274‑1 (MQ=255)
                                                                                                                             gccgtcccttGAAATTGGCATGGCGAAACTGCCTGACCATCTGCGTCCGCTAGCCGAGAAGTTACAACCCCGCCTGCGGTCGAACTCCGCTTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGAATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACcgtcg                                                                                                                                             <  4:284719/266‑1 (MQ=255)
                                                                                                                                     aTGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTAACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAACCGCCACCGCAACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggtcgt                                                                                                                                      <  2:271561/274‑1 (MQ=255)
                                                                                                                                              agaTGGCGAAACTGCGCGTTCATCTGCGGCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCGGAAGAAACGCCACCGCCACATTATTGATGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtg                                                                                                                            <  1:157600/275‑1 (MQ=255)
                                                                                                                                                  ggCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACCGCCGCCTTCGCCACCCTGTTCATGCCCGTGACCGTGGCCATTGCCGCCGCGCCAACACTCCTGTGCGTGATAgtgtgtgtga                                                                                                                        >  4:532336/1‑268 (MQ=255)
                                                                                                                                                  ggCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGCGATCGTCGCCATGACCGCCGCCGACACCGTCGGTGTCGGGATCGTGGTGGTGa                                                                                                                        >  4:70426/1‑275 (MQ=255)
                                                                                                                                                                    tcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGgcaacc                                                                                                                                                                                                                                     >  2:560314/1‑148 (MQ=255)
                                                                                                                                                                        cGTCGGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGgcaaccgcaa                                                                                                                                                                                                                                 <  4:560313/148‑1 (MQ=255)
                                                                                                                                                                          tCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGc                                                                                                                                                                                                                                                                       >  6:1045971/1‑109 (MQ=255)
                                                                                                                                                                          tCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGc                                                                                                                                                                                                                                                                       <  5:1045971/109‑1 (MQ=255)
                                                                                                                                                                          tCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACAGCCGTTGCCTTTACCCCCACCGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATTACCGCCGCACACACCGACCGGGTCGTGATCGTTGTGGGTATCGTGCGCGCCGTGAACGTGAACa                                                                                                >  4:116046/1‑275 (MQ=255)
                                                                                                                                                                                tacgggaaaGTCAAAAGCCAGCCAGCCGCCGAACATCGCGTCGCA‑GGCGGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGGGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGcc                                                                                          <  2:532337/268‑1 (MQ=255)
                                                                                                                                                                                gACCGGAAAGTTAAAGGCCAGCGGGCCGTGGACAAGAGGGTCGAA‑GGCTACGGAAACGCCACCGCCACATAATTGGGGCGTAAAAAAGGCGGGGGTTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCACCCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGcc                                                                                          <  4:131703/275‑1 (MQ=255)
                                                                                                                                                                                gACCCAGAAGGTACAAGCCAGCCAGCCGTCGAACAGCGCGTCGCA‑GGCTGAAGAACCGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGcc                                                                                          <  1:142952/275‑1 (MQ=255)
                                                                                                                                                                                 ccgggaaagggaaagGCCGGCCNNCCGTGNGACAACGCGCCGCACCGATGAGAAACCGCACCCGCACCTTTATTGTGGGCGAAAAAAAGCGGGTGTCCCCCGCTTTTTTTTCACTTTTT‑GGGTAGTGTGGGAAACCGCAACCGCCGTTGCCTTTCCCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGcc                                                                                          <  3:9147/264‑1 (MQ=255)
                                                                                                                                                                                      gaAGTAAAAAGCCATCCAGCCGTCGAACATCGCGTCGCA‑GGCGGAAGAAACGCAACCGCAACATTATTGGGGCGTAAAAAAGGCGGGGATTCCCCGCTTTTTTGTCACTTTTG‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGttc                                                                                     <  2:70427/274‑1 (MQ=255)
                                                                                                                                                                                                gccagccagccGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCtt                                                                                                                                                                                                                                                                     <  8:770350/89‑1 (MQ=255)
                                                                                                                                                                                                gccagccagccGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCtt                                                                                                                                                                                                                                                                     >  7:770350/1‑89 (MQ=255)
                                                                                                                                                                                                 ccagccagccGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAt                                                                                                                                                                                                                                                                              >  8:654102/1‑79 (MQ=255)
                                                                                                                                                                                                 ccagccagccGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAt                                                                                                                                                                                                                                                                              <  7:654102/79‑1 (MQ=255)
                                                                                                                                                                                                             cGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAACAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGc                                                             <  1:110784/275‑1 (MQ=255)
                                                                                                                                                                                                                         cGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTAtt                                                                                                                                                                                                                                               >  2:8774/1‑85 (MQ=255)
                                                                                                                                                                                                                              ggCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTg                                                                                                                                                                                                                                           <  4:8774/85‑1 (MQ=255)
                                                                                                                                                                                                                                  aaaaaaaCG‑CACCGCCAAATAAGTGAGGCGAAAAAAAAGCGGGTATTCCCCCCTTTTTTGTGACGTTGT‑AGGTATCAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTGCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAAttt                                        <  4:340873/271‑1 (MQ=255)
                                                                                                                                                                                                                                   aagaaACCCCACCCCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCTGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAAttt                                        <  2:116047/275‑1 (MQ=255)
                                                                                                                                                                                                                                         cgccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATAACCGTAATCGTGGCCAGGACCGCCGCAGCAACCATCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAAGGTGACCATGAGCCAGTTCTTCTTCGTTCGCCTCGCGCATCGCCAACACTTCAACGTTGAATTTCCGGtt                                  >  1:103054/1‑275 (MQ=255)
                                                                                                                                                                                                                                                    caTTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTCTTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCGTGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCCCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGATTCCCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGCCCGGCCa                       >  2:495774/1‑275 (MQ=255)
                                                                                                                                                                                                                                                        attGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCAGCGTGTTCATGACCGGGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCAt                   <  2:50158/275‑1 (MQ=255)
                                                                                                                                                                                                                                                         tgggggCGTAAAAAAAACGGGTGTTCCCTGCTTTGTGTGCAGTTTTT‑AGGTATTAGTGGAAACCCAAACCGCCGTTGCTCTGACGGCCACAGCATCCTCCGCCCCCGTGTTCATGACGGTGATGGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTTATCGTGGTGGTGATCGGGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCAtg                  <  4:120806/271‑1 (MQ=255)
                                                                                                                                                                                                                                                             ggtgTAAAAAAACCGGGCATCTCCCGTTTTTTTGT‑ACTTTTTCCGGTATTATTGGCAACCGGAACCGCCGTTGCCTTTCCCGCCACAGCAGCCTTCGCCACCGGGGTCATGACCGGGGTCGTGGCCTGGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGt              <  4:259043/272‑1 (MQ=255)
                                                                                                                                                                                                                                                              gCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTAGCAACCGCAACCGCCGTTGCCTTTACCTCCACAGCATCCTTCGCCACCGCTTTCATGACCGTGCACGTGGCCATGACCGCCGCAGCAACAGTCGTGGCCGTCATCGTGGTGGTCAGCGCGCGCCCCGTGAACCTGACCATGCGCCAGTTCTCCTTCAGTCGCCTCTCGATTCCCCACAACTTCAACGTTACCTTCCAGGTTCTGCCccgcccgccactcgt              >  4:441502/1‑267 (MQ=255)
                                                                                                                                                                                                                                                               cGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGAACGGCCAGCCTGGTGTTa            >  2:300325/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                 tAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGAGCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCCCAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGGGGTTAcc          >  2:324293/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                        aGCGGGGAATCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggt                                                                                                                                         >  5:948366/1‑141 (MQ=255)
                                                                                                                                                                                                                                                                        aGCGGGGAATCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggt                                                                                                                                         <  6:948366/141‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCa  <  1:41349/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gggATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCt                                                                                                                                                                                                   >  6:956173/1‑79 (MQ=255)
                                                                                                                                                                                                                                                                            gggATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGCTGCCTTTACCGCCACAGCACCCt                                                                                                                                                                                                   <  5:956173/79‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ggATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCaa                                                                                                                                                   <  5:1097255/126‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ggATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCaa                                                                                                                                                   >  6:1097255/1‑126 (MQ=255)
                                                                                                                                                                                                                                                                               aTTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCa                                                                                                                                                                  >  5:123878/1‑109 (MQ=255)
                                                                                                                                                                                                                                                                               aTTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCa                                                                                                                                                                  <  6:123878/109‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                            
TCCGCGGCGACCATTCGACAAAAGAGGTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCA‑GGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTT‑CGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCA  >  NC_000913/3477607‑3478145

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑