Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,477,878 T→A 10.6% 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,8780TA10.6% 222.2 / 15.5 75intergenic (+20/+29)slyX/slyDphi X174 lysis protein/FKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base T (37/30);  new base A (5/3);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.19e-01

AATCCGCGGCGACCATTCGACAAAAGAGGTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCA  >  NC_000913/3477605‑3478145
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                             
aaTCCGCGGCGACCATTCGACAAAAGAGGTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGACACCGCCACATTATTGAAGCGTAAAAAAAGCGGGGAt                                                                                                                                                                                                                                                                               >  2:231262/1‑274 (MQ=255)
                tCGACAAAAGAGGTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCAACGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGATTTTTTGtac                                                                                                                                                                                                                                                             >  4:53939/1‑273 (MQ=255)
                       aaGAGGTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACttttt                                                                                                                                                                                                                                                       <  2:53940/274‑1 (MQ=255)
                           ggTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAAAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTa                                                                                                                                                                                                                                                  <  2:434461/275‑1 (MQ=255)
                            gTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGGTTCCCCGCCTTTTTGTCACTTTTTCGGTAt                                                                                                                                                                                                                                                 >  1:14689/1‑275 (MQ=255)
                            gTGAATCATGCAGGATTTATCATTGGAACCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGATCGTCACGGTGACCGCTCATGAAATGGAGATGGAGAAACTGCGCGATCTTCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAACAAAAGCGGGTATTCCCCGCTTTTTTGTCACTTTTTCGGTAt                                                                                                                                                                                                                                                 >  2:360702/1‑275 (MQ=255)
                                      cAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTCTTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGGATTATTGgcaa                                                                                                                                                                                                                                        >  1:135862/1‑274 (MQ=255)
                                                                                 gACTGGCTTTTCAGGAAATACCCATGGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGTCGAAACTGCGCAATCCTCTGCGTCTGCTGACCGAGAAGTTAAGAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCAcc                                                                                                                                                                                           <  2:431010/275‑1 (MQ=255)
                                                                                          ttCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCAt                                                                                                                                                                                   <  2:303170/274‑1 (MQ=255)
                                                                                                aaaTCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCCCTTTTTCGGCATTAATTGCCACCGCAACCCCCGTTGCCTTTACCGCCACCGCATCCTTCGCCCCCCGTTCCATGACCGt                                                                                                                                                                             >  4:53156/1‑274 (MQ=255)
                                                                                                              gaaCTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCAACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTTTTAATGACCGAGATCGTGGCCATGa                                                                                                                                                               >  1:8720/1‑274 (MQ=255)
                                                                                                               aaCGGAACGTCACGGTGACCGCTCATGAAATGGGGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACCCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGAcc                                                                                                                                                             <  2:53157/275‑1 (MQ=255)
                                                                                                                      cGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCAGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGCACATCGCGACGCAGGCTGAAGAAACGCCACCGCCAACTTATTGAGGAGTAAAAAACGCGGGGATTACCCGCTTTTTTGTATCTCTTTCGGTATTAGTGGCCACTGCCACCGTCGTTACCTTTACCCCCCAATCCGCCTTCGCCACCGTGTTCAGGCCCGTAATCTTGGcccttacagtccccg                                                                                                                                                      >  4:175942/1‑262 (MQ=255)
                                                                                                                                aCCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTcc                                                                                                                                                                                                                                                                            <  6:137241/149‑1 (MQ=255)
                                                                                                                                aCCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCACTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCACACGCCGTTGCCTTTAACGCCCCAGCAGCCTTCCACACCGTGTTCATTACCGTGATCGTGGCGATGCCCGCCGAAACAACagtcat                                                                                                                                            >  4:212566/1‑273 (MQ=255)
                                                                                                                                  cGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTcccc                                                                                                                                                                                                                                                                          <  5:273187/149‑1 (MQ=255)
                                                                                                                                  cGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTcccc                                                                                                                                                                                                                                                                          <  5:116476/149‑1 (MQ=255)
                                                                                                                                    cTCATGAAATGGGGATGGCGAAACTGCCGGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggt                                                                                                                                         <  2:212569/274‑1 (MQ=255)
                                                                                                                                       aTGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCttt                                                                                                                                                                                                                                                                    >  6:151500/1‑149 (MQ=255)
                                                                                                                                       aTGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCttt                                                                                                                                                                                                                                                                    <  6:350401/149‑1 (MQ=255)
                                                                                                                                        tGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCtttt                                                                                                                                                                                                                                                                   >  6:355533/1‑149 (MQ=255)
                                                                                                                                           aaTGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAACCGCCACCGCCACATTATTGAGGAGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATc                                                                                                                                  <  2:136460/274‑1 (MQ=255)
                                                                                                                                             tGGAGATGGCGAAACTGCGCGACCATCTGCGTCTGCTGACCAAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTGTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgt                                                                                                                                <  2:175944/274‑1 (MQ=255)
                                                                                                                                                    ggCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGa                                                                                                                                    >  2:166015/1‑264 (MQ=255)
                                                                                                                                                    ggCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgtgt                                                                                                                                                                                       >  2:127306/1‑212 (MQ=255)
                                                                                                                                                         aaCTGAGCGACAATCTGAGTATGCTGACCGAGAAGTTAAAAGCCAGCGAGCCGTCGAACATCGCGTCGCAGGCGGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTg                                                                                                                   <  2:263613/275‑1 (MQ=255)
                                                                                                                                                          aCTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTAt                                                                                                                                                                                                                                                 <  5:175476/149‑1 (MQ=255)
                                                                                                                                                              cgcgATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTc                                                                                                                                                                                                 >  2:504366/1‑192 (MQ=255)
                                                                                                                                                              cgcgATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGGGGGCGTGATCGTGGTGGTGATCGTgggcg                                                                                                               >  2:501319/1‑274 (MQ=255)
                                                                                                                                                                       cTGCGTCTGCTGACAGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACg                                                                                                      <  1:139010/274‑1 (MQ=255)
                                                                                                                                                                       cTGCGTCTCCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGt                                                                                                     <  2:48325/275‑1 (MQ=255)
                                                                                                                                                                              tgctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAAccgcc                                                                                                                                                                                                                             >  6:359475/1‑149 (MQ=255)
                                                                                                                                                                              tgctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAAccgcc                                                                                                                                                                                                                             >  8:433296/1‑149 (MQ=255)
                                                                                                                                                                               cgtgACCGAGAAGTAAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCTCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATg                                                                                              <  2:195215/272‑1 (MQ=255)
                                                                                                                                                                                  gACCGAGAATTTAAAAGCCAGCCAGCCGTCGACCATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGcc                                                                                          <  2:252084/275‑1 (MQ=255)
                                                                                                                                                                                     cGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAtt                                                                                                                                                                                                                                                                              >  8:72021/1‑94 (MQ=255)
                                                                                                                                                                                     cGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAtt                                                                                                                                                                                                                                                                              <  7:72021/94‑1 (MQ=255)
                                                                                                                                                                                      gagaAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTa                                                                                                                                                                                                                                               >  7:150302/1‑123 (MQ=255)
                                                                                                                                                                                      gagaAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTa                                                                                                                                                                                                                                               <  8:150302/123‑1 (MQ=255)
                                                                                                                                                                                           gTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCg                                                                                                                                                                                                               >  8:322587/1‑149 (MQ=255)
                                                                                                                                                                                             tAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGttcttcttc                                                                               >  2:325302/1‑275 (MQ=255)
                                                                                                                                                                                                 agccagccagccGTCGAACATCGCATCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCAcag                                                                                                                                                                                                         >  8:335512/1‑149 (MQ=255)
                                                                                                                                                                                                      gccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCtttttt                                                                                                                                                                                                                                                                 >  2:343913/1‑89 (MQ=255)
                                                                                                                                                                                                         agccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGt                                                                                                                                                                                                                           <  7:191869/124‑1 (MQ=255)
                                                                                                                                                                                                         agccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGt                                                                                                                                                                                                                           >  8:191869/1‑124 (MQ=255)
                                                                                                                                                                                                          gccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCa                                                                                                                                                                                                                                                             <  4:343909/89‑1 (MQ=255)
                                                                                                                                                                                                           ccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCttttt                                                                                                                                                                                                                                                                  <  6:90883/83‑1 (MQ=255)
                                                                                                                                                                                                           ccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCttttt                                                                                                                                                                                                                                                                  >  5:90883/1‑83 (MQ=255)
                                                                                                                                                                                                                   cATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCAGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgtgt                                                                                                                                                                                       <  8:312288/149‑1 (MQ=255)
                                                                                                                                                                                                                        cgTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCCTGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTc                                                    >  2:230913/1‑275 (MQ=255)
                                                                                                                                                                                                                                gCTGAAGAAACGCCACCGCCCCATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCAAACGTGTTCNTGACCGTTATCGTGGCCATGCCCGCCGCACCAACCCTCGTGTTCGGGCACCTGGTTGTGATCGTGCCCCACCTGNACGTGACCCTGAGCCAGTTCCTCCTCTCTCNNCTCGCCGATCACCACAACtatcaatctga                                            >  3:67910/1‑265 (MQ=255)
                                                                                                                                                                                                                                  tGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCg                                                                                                                                                                        >  7:344003/1‑149 (MQ=255)
                                                                                                                                                                                                                                    aagaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTg                                                                                                                                                                      <  7:25997/149‑1 (MQ=255)
                                                                                                                                                                                                                                     agaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTgg                                                                                                                                                                     <  5:359475/149‑1 (MQ=255)
                                                                                                                                                                                                                                     agaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTgg                                                                                                                                                                     <  5:151500/149‑1 (MQ=255)
                                                                                                                                                                                                                                      gaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGc                                                                                                                                                                    >  8:238568/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         tttggggctgaaaaaaaGGAGGGATGGCCGGCTTTTTTGTCACGTTTTCGGTATGAGTGGAAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCAt                   <  2:59483/265‑1 (MQ=255)
                                                                                                                                                                                                                                                                cGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGt                                                                                                                                                                             >  4:48459/1‑114 (MQ=255)
                                                                                                                                                                                                                                                                cGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGtt             >  2:132011/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                cGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTa            >  1:139349/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                  tAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGGTCGTGGTGGTGATCGTGCGCGCCGTGACGGTGACCATGAGCCAGGTTTTCTTCATTCGCTTCGCGAACCGCCACACCTTCCACGTTGAATTTCCGTTTCTGACCCCGCCGCATTTGGTTAcc          >  4:307631/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                  tAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTAcc          >  2:361947/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                  tAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGGGACCATGAGCCCGTTCTTCTTCAGTCGCTTCGCGAAGCGCACCAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCCGCATGTGGTTAcc          >  2:329631/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                  tAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAAACGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTGGTGATCGTGGTGGTCCACGTGCGCGCCGGGAACGTGGCCAGGATCCGGTTCTTCTTCAGTCGCTTCGCGAATCGCCACACCATCAACGTTGAATTTCAGGTTCTAGCCGGCTAGcctagtgttccc          >  4:76762/1‑264 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaagGCGGTGTTTCCCCCTGTTTGTTGCAGTTTGTGGGTATGAGCGTCAACCGTAGCCCTCGTTGTCTCTGCCGCCACGCGAGCGTGCGCCACCGTGTTCATGACCGTGAGCGTGGCCATGTCCGCCGAAGCACCCGTGGTGTTCGTGATCGTGGTGGTGATCGTGGGCGCCGGGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCAt        <  4:286858/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaaaGCTGGGATTCCCCGTTTGTTGTTCCCTTTTTCGTTTTTGGTGCTACCCCGAACACCCGTTGCATTTCCCACCACGGCAGCCTCCGCCACGTTGGTCATGACCGTGATAGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTAATCGTGGTGGTGAGCGTGCGCGCCGTGAACGTGACCATGTGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCGGGTTCTGACCGGCCAGCATGTGGTTACCAt        <  4:325298/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaaaGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATc                                                                                                                                                                         <  2:48459/114‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaaaGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggtcgt                                                                                                                                      >  5:411464/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaaaGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggtcgt                                                                                                                                      >  5:296609/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaaaGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggtcgt                                                                                                                                      >  8:433653/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaaaGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggtcat                                                                                                                                      >  6:38245/1‑147 (MQ=255)
                                                                                                                                                                                                                                                                     aaaaaGCGGGGAATCCCCGCTTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggtcgt                                                                                                                                      >  5:84922/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      aaaaGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGa                                                                                                                                    >  7:2996/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        aaGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATc                                                                                                                                  <  8:138048/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCg                                                                                                                                 >  5:455305/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          gCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTTACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCCAAGCGCCACAACTTCAACGTTTAATTTCAGGTTCTGACCGGCCAGCATTTCGTTACCATCACCCa  >  2:159318/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                             gggAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCATGATCgtggt                                                                                                                             <  5:38245/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ggAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtg                                                                                                                            <  6:296609/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                             
AATCCGCGGCGACCATTCGACAAAAGAGGTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCA  >  NC_000913/3477605‑3478145

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