Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 422,247 A→T 15.3% Y421F (TAT→TTT)  proY → proline‑specific permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913422,2470AT15.3% 155.2 / 3.0 66Y421F (TAT→TTT) proYproline‑specific permease
Reads supporting (aligned to +/- strand):  ref base A (13/40);  new base T (8/3);  total (23/43)
Fisher's exact test for biased strand distribution p-value = 3.76e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCA‑AATATCCTTCAGATGATGAGTGATCCTGCATTAGG  >  NC_000913/422017‑422476
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                       
ttcGGTGTTTTTTAACTACATCCTGGCGGAAAACGTTTTCCTGGTGATTCCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCgg                                                                                                                                                                                                                                 <  1:323909/237‑1 (MQ=255)
 gCGGTGTATCTGAACTCCATCTTGCCGTAAACCGTCTTCCGGGTTATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGa                                                                                                                                                                                                                                <  1:89895/239‑1 (MQ=255)
 gCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGa                                                                                                                                                                                                                                <  2:658671/239‑1 (MQ=255)
     tgtATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCCACGACCATCGGCGGGCTGCGTTTCATC‑GTC‑TTTATTTTTCGGGGGAT‑TGGGTTTCACCCCGATAcg                                                                                                                                                                                                                            >  2:323960/1‑239 (MQ=255)
       tATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGCTTTTTCTG‑CTC‑TTTTTTATCTGGTGTGT‑TGGTATTCACCCGGAAAcacc                                                                                                                                                                                                                          >  2:103112/1‑236 (MQ=255)
       tATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCAc                                                                                                                                                                                                                                     <  2:271766/228‑1 (MQ=255)
              actacattcatccGGAAAACGTTCTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCg                                                                                                                                                                                                                   <  2:414678/228‑1 (MQ=255)
               cTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGCGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGCCGGGCAGCGTTTCCTG‑CTC‑CTTATTATCCGGGTGTG‑TTGGTCTCACCCGCATAagcggaatat                                                                                                                                                                                                                    >  1:604839/1‑231 (MQ=255)
                 aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGGCGGTTTTTCTT‑CTCTTTTTTTTTGGGGTTGTT‑GGGTTATCCCCCgaaacgccgcgctccttt                                                                                                                                                                                                                 >  1:271766/1‑222 (MQ=255)
                 aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTTCTC‑CCC‑TTTATTTTCCGGGTGCT‑TGTGTTTTACCCGGCTaacccgttctctc                                                                                                                                                                                                                  >  2:104925/1‑225 (MQ=255)
                  catcatGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTC‑CTC‑TTTTTTATCCGGTTGAT‑TGGTTATCACCCGGATACGCCTGtcttcccgg                                                                                                                                                                                                               >  2:651805/1‑233 (MQ=255)
                       tGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGGTTTTCCTGTCCT‑TTTTTTATC‑GGTTGAT‑TGGTTTTAACCCgggataccggattctcctgcatttc                                                                                                                                                                                                          >  2:408068/1‑216 (MQ=255)
                       tGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGGTTTTCCTG‑CTC‑TTTATTATCTGGGTGAT‑TGGGTATCACCCCGGTACGCGGATTTCGCTGTATTTc                                                                                                                                                                                                          >  2:105511/1‑239 (MQ=255)
                       tGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTT‑TTTTTTATTGGGGTGATGTGGGT‑TTACCCGGGTACGCGGTTTACGCTGTATGtg                                                                                                                                                                                                          >  2:90515/1‑238 (MQ=255)
                               agcGTTTTCCTGGTGATCGCTTCGCTGGGAACCTTTGCCACGGTTGGGGTGGGGATTATGATCCTGCTGGCGCAAATTGCCTTCCGTCGCCGTTTGGCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTcgc                                                                                                                                                                                                  <  1:651805/237‑1 (MQ=255)
                               agcGTCTTCCTGGGTATCGCTTCGCTGGCACCCGTGGCCAGGGTGTGGGTGTGGATGATGAGCCGGAGTTCGCAACTTGCCTTCCGGTGCGGTTGGCCGCCAGAAGAAGTTAAGGCGCTGAAATTGAAAGGTCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTcgc                                                                                                                                                                                                  <  1:704429/237‑1 (MQ=255)
                                    cTTCCTGTTGTTCCCTTCGCTGGCAACCTTCGCCACGGGTGGGGTGTGGTTTATGATCCTGGTGTCGCAAATTGCCTTCCGTCGCCGGTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGa                                                                                                                                                                                             <  1:608940/239‑1 (MQ=255)
                                             ggtcgcttcgctGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTTCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGc                                                                                                                                                                                    <  2:34142/237‑1 (MQ=255)
                                              aTCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGGTTTTTCTG‑CTC‑TTTTTTATCGGGGTGGT‑TGGGTTTTACCCGGGTACGCGGATCTCCCGCTATGTCGGgtttccgggggttttttttgg                                                                                                                                                                                   >  2:608940/1‑221 (MQ=255)
                                              aTCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGc                                                                                                                                                                                    <  2:708256/238‑1 (MQ=255)
                                              aTCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGGTGAT‑TTGTTTTTACCCCGGTACGCCTATCTCGCTGTATTTCGGGTTCGCGTGGGtttgtgt                                                                                                                                                                                      >  2:650064/1‑231 (MQ=255)
                                               tcgcttcgctGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGGTGTCGCAAATTGCCTTCCGGCGCCGTTTGCCGGCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCt                                                                                                                                                                                   <  2:294866/238‑1 (MQ=255)
                                                  cttcgctGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGCTTTTTCTC‑CTC‑TTTTTTTTTGGGTTGAT‑TGGTTTTCACCCGGATACGCGTAACTCGCTGTATGTCGGGTTTCCGGGGGTTGTtgtggggggg                                                                                                                                                                               >  1:470628/1‑233 (MQ=255)
                                                              ccTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTAtc                                                                                                                                                                                                                      <  2:375730/188‑1 (MQ=255)
                                                                            tgggtgtggATTATGATCCGTCTTGCGCAAATTGCCTTCCGGCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGGTCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCa                                                                                                                                                     <  2:127985/239‑1 (MQ=255)
                                                                             gggtgtggATTAAGATCCTGCTGTCTCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCa                                                                                                                                                     <  2:345247/238‑1 (MQ=255)
                                                                                    gATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGGTGAGTTTTCTT‑CTC‑TTTTTTAATTGGGTGAT‑TTGTTTTCACCCGGGTACCGGGATCTCCCTTGTTTTCGTGTTTGCGGGGGTTGTTGTTCTGTTGTTTGGGGTGGGGTTTAAAACCCGCCaacggcggc                                                                                                                                             >  1:288739/1‑231 (MQ=255)
                                                                                      ttATGATCCTGCTTTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCa                                                                                                                                            <  2:844082/238‑1 (MQ=255)
                                                                                       taagATCCTGGTTTCGAAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAgc                                                                                                                                          <  2:682212/236‑1 (MQ=255)
                                                                                        aTGAGCCTGCTGTCGCAAATTGCCTTCCGTGGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAgct                                                                                                                                         <  1:127996/239‑1 (MQ=255)
                                                                                               ggctgTCGAAAATTGCCTTCCGGCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGaa                                                                                                                                 <  2:816419/238‑1 (MQ=255)
                                                                                                      gctaatgtccTTCCGTGGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAATTGCCGGGTGGGGTAGCAACGCCGATCGGCGGGGTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTAGCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGGGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGt                                                                                                                          <  1:166860/231‑1 (MQ=255)
                                                                                                      cgaaaTTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGt                                                                                                                          <  2:167633/237‑1 (MQ=255)
                                                                                                                tttctcggCCTTTTGCCCCCAGAAAAAGTTAAGGGGGTGAAATTTAAAGTGCCGGGTGGGGTGACAACGTCCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCg                                                                                                                <  1:323960/232‑1 (MQ=255)
                                                                                                                                 ccatgaaaaaTTTAAGGGGCTGAATTTTAAAGGGCCGGGTGGTGTGGCAACGCCCATCGGCGGGGTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGATGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCa                                                                                               <  1:650064/235‑1 (MQ=255)
                                                                                                                                      aagaagTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCaggagg                                                                                          <  2:288739/239‑1 (MQ=255)
                                                                                                                                                     tGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTCCAGGAGGAAATCCTGATCCGGc                                                                           <  2:470628/239‑1 (MQ=255)
                                                                                                                                                         aTTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGc                                                                                                                                                                                    <  2:269800/131‑1 (MQ=255)
                                                                                                                                                                         tgggTAGCAACGGCTATTGGGTGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTTAC‑TGGTTATCCCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGc                                                        <  1:278808/237‑1 (MQ=255)
                                                                                                                                                                                      ccttCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctc                                           <  2:521723/235‑1 (MQ=255)
                                                                                                                                                                                      ccgtCGGCGGGCTGTTTTTCCTG‑CTC‑TTTTTTAGTGGGGTGGT‑TGGTTATCTCCCGGATACGCGTAGCTCGTTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctcc                                          <  1:105511/236‑1 (MQ=255)
                                                                                                                                                                                      ccATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctcc                                          <  2:889699/239‑1 (MQ=255)
                                                                                                                                                                                                gCTGATTTTTCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCCCCCGGGTACGCGTATCTCGCTGTATGTCGGGTTTGCGGGGGTTGTTTGTCTGTTTGTTTGGTGGGTGTTTAAAACCCCCCCCGCGTCTCAGCTGGCTTGAAACCCAG‑AAACCCTGCGGTCTTGTAGCCGGTTCAGGAGGCAATCCTCTTCCGGCCTTTTCGAAAATTTTCCCCCCCCCCCCCCCTCCCCCc‑aaaaa                                >  1:889699/1‑235 (MQ=255)
                                                                                                                                                                                                gCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGa                                                                                 <  2:558571/190‑1 (MQ=255)
                                                                                                                                                                                                   taTTTTCCTG‑CTT‑TTTATTATCGGGTTTAT‑TGGTTACCCCCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCA‑AATATCCtt                            <  2:20417/238‑1 (MQ=255)
                                                                                                                                                                                                     ttttCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCCCCCGGATACGCGTATCTCGCTGTATGGCGGTTTCGCGGGGATTGTTGGGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCA‑AATATCCTTCa                          <  2:403482/239‑1 (MQ=255)
                                                                                                                                                                                                      tttCCTG‑CTT‑TTTATTATCGGGTTGTT‑TGGTTATCCCCCGGATACGCGTATCTCGCTGTATGGCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCA‑AATATCCTTCAg                         <  2:600586/239‑1 (MQ=255)
                                                                                                                                                                                                        tCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATTCGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAGCGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcct                                            <  2:691290/219‑1 (MQ=255)
                                                                                                                                                                                                               tc‑tTTTTTTTCTGGTTGCT‑TGGTTTTCGCCCTGATCGGCGTTTCTCGGTGTATGTCGGTTGTGCGGTGATTGTTGGGCTGTTGATTGGCTGGATGCTTAAACGCCGTCACGATCGGCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCTGCTATTCCGAAAGTTATCCGCCCCCTTCCTCCTCCCCCA‑AATATCCTTCAGATGATGAg                 <  1:866069/239‑1 (MQ=255)
                                                                                                                                                                                                               tc‑tTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGccctcctcctgtctcttttacccatctccgag                   >  2:883737/1‑215 (MQ=255)
                                                                                                                                                                                                                    tattatCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAGCGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGGGGAAATCCTGATCCGGCCATTCCGAAAGTTATCCGCCCCCGTCTCCCTCCCCCA‑AAGATctgtctctttaacaaaacg             >  1:691290/1‑221 (MQ=255)
                                                                                                                                                                                                                      ttatCGGGTTGAT‑TGGTTATCACCGGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCA‑AATATCCtt                            <  2:646734/222‑1 (MQ=255)
                                                                                                                                                                                                                      ttatCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTAGGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCA‑AATATCCTTCAGATGATGAg                 <  2:780540/233‑1 (MQ=255)
                                                                                                                                                                                                                           ggggTGAT‑TGGGTATCACCCGGAATCGCGTATCTTGCTGTATGTCGGTTGCGCGGGGATTGTTGGGCGGGTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCA‑AATATCCTTCAGATGATGAGTGATCCTGCAt      <  2:605668/239‑1 (MQ=255)
                                                                                                                                                                                                                           gggTTGTT‑TGGTTTTCCCCCGGATCCGCGTATCTCGCTTTATTTCGGTTTCGCGTGGATTGGTTTGCTTTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCA‑AATATACTTCAGATGATGAGTGATCCTGCAt      <  1:785674/239‑1 (MQ=255)
                                                                                                                                                                                                                           gggTTGTT‑TGGTTATCCCCCGGTTGCGCGTCTTTGGCTGTATGTCGGTTTCGCGTGGATTGTTTGGCTGTTGATTGGCTGGATGTTTAAAGGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCA‑AATATCCTTCAGATGATGAGTGATCCTGCAt      <  2:242674/239‑1 (MQ=255)
                                                                                                                                                                                                                           gggTTGAT‑TGGTTTTCCCCCGGATACGCGTATCTCGGTGTATGTCGGTTTCGCGTGGATTGTTGGGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCA‑AATATCCTTCAGATGATGAGTGATCCTGCAt      <  2:269073/239‑1 (MQ=255)
                                                                                                                                                                                                                            ggTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGc                                                                           <  2:661064/170‑1 (MQ=255)
                                                                                                                                                                                                                            ggTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGc                                                                           >  1:661064/1‑170 (MQ=255)
                                                                                                                                                                                                                            ggTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCC‑AATATCCTTCAGATGTTGAGTGGTCCTGCAtt     >  2:828248/1‑239 (MQ=255)
                                                                                                                                                                                                                            ggTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCA‑CATATCCTTCAGATGATAAGTGATCCTCCAta     >  1:14402/1‑238 (MQ=255)
                                                                                                                                                                                                                            ggTTGAT‑TGGCTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTAGGGAGAGAAGCCGGAGCCGGCTCTTCTGAAAGTTACCCGCGCCCG‑CCCCCTCCCCCACAAAATCCTTCaatggattagggaaccggaatt     >  1:357952/1‑218 (MQ=255)
                                                                                                                                                                                                                                gAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTTTCAGCCGGTTCAGTCGGAAATCCTGGTCCGGCTATTCCGAAAGTTATCCCCCCCCGTCCCCCTCCCCCC‑ACTATCCTTCAgaatctacggtgtcccacacctgt  >  1:780540/1‑216 (MQ=255)
                                                                                                                                                                                                                                gAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGGGGAAATCCTGATCCGCCTATTCCGAAAGTTATCCGCCCCCGTACTCCCCCCCCA‑CATATCCTTCACATGTTGAGttgtctcttcaacaa  >  1:646734/1‑224 (MQ=255)
                                                                                                                                                                                                                                gAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCA‑AATATACTTCAGATGATGAGTGATCCTGCATTAg   >  2:785674/1‑237 (MQ=255)
                                                                                                                                                                                                                                       ttATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCA‑AATATCCTTCAGATGATGAg                 <  2:437532/217‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                       
TGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTAGCAACGACCATCGGCGGGCTGATTTTCCTG‑CTC‑TTTATTATCGGGTTGAT‑TGGTTATCACCCGGATACGCGTATCTCGCTGTATGTCGGTTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGC‑TGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCA‑AATATCCTTCAGATGATGAGTGATCCTGCATTAGG  >  NC_000913/422017‑422476

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑