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

TCCGCGGCGACCATTCGACAAAAGAGGTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCA  >  NC_000913/3477607‑3478145
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                            
tCCGCGGCGACCATTCGCCAAAAGAGGTGAATCATGCAGGATTTATCAGTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTCCGTGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATAAAAGGGAGATGGCGAAACTGCGCGAGCATCTGCCTCTGCTGACGGAGATGTTAAAAACCAGCTAGCCGTCGAACATCCCGACCGATGTTCAAGACACGCCACCGCAACACACTTAAGGACTAAAAAAAACGGGGAtac                                                                                                                                                                                                                                                                            >  4:514551/1‑272 (MQ=255)
tCCGCGGCGACCATTCGACAAAAGAGGTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCAAAGTATTGAGGCGTAAAAAAAGCGGGGATTc                                                                                                                                                                                                                                                                            >  2:284720/1‑274 (MQ=255)
                                  tGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGGGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgca                                                                                                                                                                                                                                        <  2:85904/275‑1 (MQ=255)
                                                 tGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTTAGGAGTAAAAAAAGCGGGGGTTCCACGCTTTTTTGTCACTTTTTCGGTATTAGTGTCAACCCCAACCGCCAtt                                                                                                                                                                                                                         >  1:105781/1‑275 (MQ=255)
                                                 tGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGtt                                                                                                                                                                                                                         >  2:366256/1‑275 (MQ=255)
                                                 tGGAAGCACGCCTGGCAGAGCTGGAGAGCCGAATGGCTTTTCAGGAAATCACCATTGAATAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGAACGAGAAGTNAAAAGCCAGCAACCCGTCGAACATCGACTAGCAGGCTGAAGAAACGCCACCGCCAACATAATGAGGAGTAAAAAAAGCGGGGATTCACCGATCTTTTGCCACTTTTNc                                                                                                                                                                                                                                                     >  3:110784/1‑247 (MQ=255)
                                                     aGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCCCGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTACAAGCCAGCCAGCCGTCGAACATCCCGTCGCAGGCCACATAAACCCCACCCCCACATTATTCAGGCCTACAAAACCCCGGGACTCCCCGCTTTTTTGGTCCTTTTTCTGTATTAGTGGCAACCCCACCCCCCGTTGccc                                                                                                                                                                                                                     >  4:516335/1‑274 (MQ=255)
                                                            tgtgaagagCGGGATGGCCGCTTGCCTTTTCAGGAAATCCCCCTTTAATAACTGAACGTCCCGGTGACCGGCCATGAAATGGAGCTGGGGAACCTGCTCCATCATCTGCGTCGGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGc                                                                                                                                                                                                              <  4:404958/270‑1 (MQ=255)
                                                               gCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCca                                                                                                                                                                                                            >  2:472254/1‑274 (MQ=255)
                                                                    gCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAACAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCCACCGCAAACCCCTTTGCCTTTTCCGcccccgccc                                                                                                                                                                                                      >  4:105008/1‑268 (MQ=255)
                                                                                aCTGGCTTTTCAGGAAATCACCATTGAAGAACTGCACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTCAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTATGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGGTCCCTTTACCGCCACAGCAGCCTTAGcacccg                                                                                                                                                                                          >  2:439645/1‑270 (MQ=255)
                                                                                 cgggtatttCAGTAATACGCCATTGGANNCATGAACGTCACGGTGCGTGATCATGAAAGNGAGAGGGCGAAACTGCGCGTTCATCTGCGTCTGCTGACCGAGAAGTAACAAGCAAGCCAACCGTCAACAATCGANTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCAcc                                                                                                                                                                                           <  3:57955/267‑1 (MQ=255)
                                                                                                 taaCCATTGAAGAACGGAACGTCACGGTGACCGCTCATGAAATGGAGAGGGCGAAACTGCGCGATCGTCTGCGTCGGCTGACCGGGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATc                                                                                                                                                                         <  2:4103/273‑1 (MQ=255)
                                                                                                                    cGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCCGGACCGCCgcag                                                                                                                                                      >  1:9147/1‑275 (MQ=255)
                                                                                                                            tGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACcgtc                                                                                                                                              >  2:272381/1‑275 (MQ=255)
                                                                                                                             taCCGCTCATGAAAGGGAAATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACcgtcg                                                                                                                                             <  2:105009/274‑1 (MQ=255)
                                                                                                                             gccgtcccttGAAATTGGCATGGCGAAACTGCCTGACCATCTGCGTCCGCTAGCCGAGAAGTTACAACCCCGCCTGCGGTCGAACTCCGCTTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGAATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACcgtcg                                                                                                                                             <  4:284719/266‑1 (MQ=255)
                                                                                                                                     aTGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTAACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAACCGCCACCGCAACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggtcgt                                                                                                                                      <  2:271561/274‑1 (MQ=255)
                                                                                                                                              agaTGGCGAAACTGCGCGTTCATCTGCGGCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCGGAAGAAACGCCACCGCCACATTATTGATGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtg                                                                                                                            <  1:157600/275‑1 (MQ=255)
                                                                                                                                                  ggCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACCGCCGCCTTCGCCACCCTGTTCATGCCCGTGACCGTGGCCATTGCCGCCGCGCCAACACTCCTGTGCGTGATAgtgtgtgtga                                                                                                                        >  4:532336/1‑268 (MQ=255)
                                                                                                                                                  ggCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGCGATCGTCGCCATGACCGCCGCCGACACCGTCGGTGTCGGGATCGTGGTGGTGa                                                                                                                        >  4:70426/1‑275 (MQ=255)
                                                                                                                                                                    tcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgcaacc                                                                                                                                                                                                                                     >  2:560314/1‑148 (MQ=255)
                                                                                                                                                                        cGTCGGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgcaaccgcaa                                                                                                                                                                                                                                 <  4:560313/148‑1 (MQ=255)
                                                                                                                                                                          tCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGc                                                                                                                                                                                                                                                                       >  6:1045971/1‑109 (MQ=255)
                                                                                                                                                                          tCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGc                                                                                                                                                                                                                                                                       <  5:1045971/109‑1 (MQ=255)
                                                                                                                                                                          tCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACAGCCGTTGCCTTTACCCCCACCGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATTACCGCCGCACACACCGACCGGGTCGTGATCGTTGTGGGTATCGTGCGCGCCGTGAACGTGAACa                                                                                                >  4:116046/1‑275 (MQ=255)
                                                                                                                                                                                tacgggaaaGTCAAAAGCCAGCCAGCCGCCGAACATCGCGTCGCAGGCGGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGGGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGcc                                                                                          <  2:532337/268‑1 (MQ=255)
                                                                                                                                                                                gACCGGAAAGTTAAAGGCCAGCGGGCCGTGGACAAGAGGGTCGAAGGCTACGGAAACGCCACCGCCACATAATTGGGGCGTAAAAAAGGCGGGGGTTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCACCCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGcc                                                                                          <  4:131703/275‑1 (MQ=255)
                                                                                                                                                                                gACCCAGAAGGTACAAGCCAGCCAGCCGTCGAACAGCGCGTCGCAGGCTGAAGAACCGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGcc                                                                                          <  1:142952/275‑1 (MQ=255)
                                                                                                                                                                                 ccgggaaagggaaagGCCGGCCNNCCGTGNGACAACGCGCCGCACCGATGAGAAACCGCACCCGCACCTTTATTGTGGGCGAAAAAAAGCGGGTGTCCCCCGCTTTTTTTTCACTTTTTGGGTAGTGTGGGAAACCGCAACCGCCGTTGCCTTTCCCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGcc                                                                                          <  3:9147/264‑1 (MQ=255)
                                                                                                                                                                                      gaAGTAAAAAGCCATCCAGCCGTCGAACATCGCGTCGCAGGCGGAAGAAACGCAACCGCAACATTATTGGGGCGTAAAAAAGGCGGGGATTCCCCGCTTTTTTGTCACTTTTGCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGttc                                                                                     <  2:70427/274‑1 (MQ=255)
                                                                                                                                                                                                gccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCtt                                                                                                                                                                                                                                                                     <  8:770350/89‑1 (MQ=255)
                                                                                                                                                                                                gccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCtt                                                                                                                                                                                                                                                                     >  7:770350/1‑89 (MQ=255)
                                                                                                                                                                                                 ccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAt                                                                                                                                                                                                                                                                              >  8:654102/1‑79 (MQ=255)
                                                                                                                                                                                                 ccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAt                                                                                                                                                                                                                                                                              <  7:654102/79‑1 (MQ=255)
                                                                                                                                                                                                             cGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAACAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGc                                                             <  1:110784/275‑1 (MQ=255)
                                                                                                                                                                                                                         cGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTAtt                                                                                                                                                                                                                                               >  2:8774/1‑85 (MQ=255)
                                                                                                                                                                                                                              ggCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTg                                                                                                                                                                                                                                           <  4:8774/85‑1 (MQ=255)
                                                                                                                                                                                                                                  aaaaaaaCGCACCGCCAAATAAGTGAGGCGAAAAAAAAGCGGGTATTCCCCCCTTTTTTGTGACGTTGTAGGTATCAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTGCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAAttt                                        <  4:340873/271‑1 (MQ=255)
                                                                                                                                                                                                                                   aagaaACCCCACCCCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCTGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAAttt                                        <  2:116047/275‑1 (MQ=255)
                                                                                                                                                                                                                                         cgccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATAACCGTAATCGTGGCCAGGACCGCCGCAGCAACCATCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAAGGTGACCATGAGCCAGTTCTTCTTCGTTCGCCTCGCGCATCGCCAACACTTCAACGTTGAATTTCCGGtt                                  >  1:103054/1‑275 (MQ=255)
                                                                                                                                                                                                                                                    caTTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTCTTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCGTGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCCCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGATTCCCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGCCCGGCCa                       >  2:495774/1‑275 (MQ=255)
                                                                                                                                                                                                                                                        attGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCAGCGTGTTCATGACCGGGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCAt                   <  2:50158/275‑1 (MQ=255)
                                                                                                                                                                                                                                                         tgggggCGTAAAAAAAACGGGTGTTCCCTGCTTTGTGTGCAGTTTTTAGGTATTAGTGGAAACCCAAACCGCCGTTGCTCTGACGGCCACAGCATCCTCCGCCCCCGTGTTCATGACGGTGATGGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTTATCGTGGTGGTGATCGGGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCAtg                  <  4:120806/271‑1 (MQ=255)
                                                                                                                                                                                                                                                             ggtgTAAAAAAACCGGGCATCTCCCGTTTTTTTGTACTTTTTCCGGTATTATTGGCAACCGGAACCGCCGTTGCCTTTCCCGCCACAGCAGCCTTCGCCACCGGGGTCATGACCGGGGTCGTGGCCTGGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGt              <  4:259043/272‑1 (MQ=255)
                                                                                                                                                                                                                                                              gCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTAGCAACCGCAACCGCCGTTGCCTTTACCTCCACAGCATCCTTCGCCACCGCTTTCATGACCGTGCACGTGGCCATGACCGCCGCAGCAACAGTCGTGGCCGTCATCGTGGTGGTCAGCGCGCGCCCCGTGAACCTGACCATGCGCCAGTTCTCCTTCAGTCGCCTCTCGATTCCCCACAACTTCAACGTTACCTTCCAGGTTCTGCCccgcccgccactcgt              >  4:441502/1‑267 (MQ=255)
                                                                                                                                                                                                                                                               cGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGAACGGCCAGCCTGGTGTTa            >  2:300325/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                 tAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGAGCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCCCAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGGGGTTAcc          >  2:324293/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                        aGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggt                                                                                                                                         >  5:948366/1‑141 (MQ=255)
                                                                                                                                                                                                                                                                        aGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggt                                                                                                                                         <  6:948366/141‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCa  <  1:41349/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gggATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCt                                                                                                                                                                                                   >  6:956173/1‑79 (MQ=255)
                                                                                                                                                                                                                                                                            gggATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGCTGCCTTTACCGCCACAGCACCCt                                                                                                                                                                                                   <  5:956173/79‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ggATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCaa                                                                                                                                                   <  5:1097255/126‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ggATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCaa                                                                                                                                                   >  6:1097255/1‑126 (MQ=255)
                                                                                                                                                                                                                                                                               aTTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCa                                                                                                                                                                  >  5:123878/1‑109 (MQ=255)
                                                                                                                                                                                                                                                                               aTTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCa                                                                                                                                                                  <  6:123878/109‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                            
TCCGCGGCGACCATTCGACAAAAGAGGTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCA  >  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: