Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,123,412 T→A 13.4% I215F (ATT→TTT)  grxB ← glutaredoxin 2 (Grx2)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,123,4120TA13.4% 155.0 / 2.4 71I215F (ATT→TTT) grxBglutaredoxin 2 (Grx2)
Reads supporting (aligned to +/- strand):  ref base T (11/48);  new base A (8/2);  total (21/50)
Fisher's exact test for biased strand distribution p-value = 2.86e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

TTAATTTCTTGTTCGGTGATGGTGTATTGGGTGAGTTGATTACAGCCAACGAGCAGGCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAG  >  NC_000913/1123176‑1123628
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                           
ttAATTTTTTTTGCGGTGTTGGTGTATTGGGTGAGTTGATTACAGCCAACGAGCAGGCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCg                                                                                                                                                                                                                           <  2:1319932/239‑1 (MQ=255)
    tttCTTGTTCGGTGATGGTGTATTGGGTGAGTTGATTACAGCCAACGAGCAGGCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGACAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCAGCGGCGGGA‑CC‑GGA‑TTTAATCG‑ccct                                                                                                                                                                                                                      >  1:1193359/1‑237 (MQ=255)
         tGTTCGGTGATGGTGTATTGGGTGAGTTGATTACAGCCAACGAGCAGGCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGCCGGCAACCCCCGAATACAGGCTGCCCCGGCGGGG‑CC‑GGA‑TTAAAACG‑CCCTtgatg                                                                                                                                                                                                                 >  1:1255141/1‑239 (MQ=255)
          tttgggTGTTGGGTTATTGGGTGTGTTGATATCAGCCAACGGGCCGGCCATTGACGATCAATGACGCGGCAAATAAAATCTTGTTCATGGTAGTCCTCGACATGAAATCTGAGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATtgatga                                                                                                                                                                                                                <  2:394907/235‑1 (MQ=255)
                  aTGGTGTATTGGGTGAGTTGATTAAAGCCAACGAGCAGGCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACaaa                                                                                                                                                                                                         <  1:895402/238‑1 (MQ=255)
                   tGGTGTATTGGGTGATTTGATTACAGCCAACGGGCAGGGCGCTGACGATCAATGCAGCGGCAAATAAAAACTTGGTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAAt                                                                                                                                                                                                        <  1:520428/238‑1 (MQ=255)
                          ttGGGTGAGTTGATTACAGCCAACGAGCAGGCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CC‑GGA‑TTTAAAAC‑CCCCAGATGATAACAAATTGAttttg                                                                                                                                                                                                >  1:331921/1‑237 (MQ=255)
                          ttGGGTGAGTTGATTACAGCCAACGAGCAGGCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGG‑CA‑GGG‑TTAAAATG‑CCCTCTATGATAACAAATTGAttttg                                                                                                                                                                                                >  1:1171523/1‑237 (MQ=255)
                          ttGGGTGAGTTGATTACAGCCAACGAGCAGGCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCATCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGTCGGCAACCCGCGAATAAAGGCTGCCGTGGCGGGG‑CG‑GGG‑TTAAAATG‑CCATTGTTGATAACAAAttacatttt                                                                                                                                                                                                >  1:1081610/1‑232 (MQ=255)
                               tGAGTTGATTACAGCCAACGAGCAGGCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCCGCCGCGGGGGGG‑CA‑CGAGTAAAACCC‑CCATTGATGAAAAAAAAAT‑ATTtgtgggtg                                                                                                                                                                                           >  2:707880/1‑235 (MQ=255)
                                  gTTGATTACAGCCAACGAGCAGGCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTCCCCCGGCGGGG‑CA‑GGA‑TGAAAACA‑CCCATGATGAAAACaaaaaaatttgggttgtttt                                                                                                                                                                                        >  2:1275418/1‑222 (MQ=255)
                                    ttaTTACACCCAACGGGCAGGCCACTGACGATCAATGCAGCGGAAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCg                                                                                                                                                                                      <  2:1274404/237‑1 (MQ=255)
                                    tGATTAAAGCCAACGAGCAGGCCACTGACGATCAATGCCGCGGAAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCGGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCg                                                                                                                                                                                      <  2:411703/239‑1 (MQ=255)
                                       tttcagcacaccggaaGGGCATTGGCGATCAATGCAGCGGAAAATAAAAACTTGTTCATGGTGCGCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCa                                                                                                                                                                                   <  2:649823/224‑1 (MQ=255)
                                       ttACAGCCAACGAACAGGCCAAGGACGATCAATGCAGCGGCAAATAAAAACTTGTTCAGGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCa                                                                                                                                                                                   <  1:148063/239‑1 (MQ=255)
                                        tACAGCCAACGAGCAGGCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACCGGGTGCCGCGGCGGGG‑GA‑GGG‑TGAAAAAGCCCATTGGTGAAAAAAAAtaatttttgttgtttttgccaa                                                                                                                                                                                   >  2:1168044/1‑218 (MQ=255)
                                             ccAACGAGCAGGCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAAAAG‑CCATTGATGATAACAAATTGATTTTTGTCTGTTTTCCCAcatta                                                                                                                                                                              >  1:567628/1‑238 (MQ=255)
                                               agaggccgggcccatGACGATAAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACAGAAAATCGGCGTCAATATCCTGACACAACGCAGCATGTGTCCCCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATcg                                                                                                                                                                           <  2:1343252/225‑1 (MQ=255)
                                               aaCGGGCAGGCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGCCATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATcg                                                                                                                                                                           <  1:159106/239‑1 (MQ=255)
                                                 cGAGCAGGCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATcgcg                                                                                                                                                                         <  1:335458/239‑1 (MQ=255)
                                                       ggCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCCCCAGCAGAAAAAACTGAAAACGGCGGCAACCCCCGCATACAGGCTGCCGCGGCGGGTCCA‑GAA‑TAAAAACC‑CCATT‑ATGATTACAAAATGTTTTTTTTCTTCTTCGCCCCATTTTTGCGGGaaa                                                                                                                                                                    >  2:685447/1‑237 (MQ=255)
                                                             tGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATAAAGGCTGCCGCGGCGGGG‑CAGGGT‑TTAAACCC‑CCATTTTGGAAAAAAAATT‑CTTTTTGTTTTTTTTGCCataattttccgggaataaaaaa                                                                                                                                                              >  2:516468/1‑219 (MQ=255)
                                                                  aTCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGG‑TTAAATCG‑CCATTTAGGTTAAAAAATTGATTTGTGTCTGGTTTGCCCAATTATTCCGGTAATCAGCAAagcg                                                                                                                                                          >  1:461814/1‑237 (MQ=255)
                                                                   tCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTGGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGc                                                                                                                                                        <  1:345761/238‑1 (MQ=255)
                                                                      aTGCAGCGGCAAATAAAAACTTGTTCATGGCAATCCTCGGCATGAAATCTGCGTCAACATCCTGACACAACGCCGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTgg                                                                                                                                                    <  1:87388/239‑1 (MQ=255)
                                                                        tcAGCGGCAAAAAAAAACTTGTTCCAGGTAGGCCTCGACATGAAATCGGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACCCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGcc                                                                                                                                                  <  2:1255141/238‑1 (MQ=255)
                                                                                        cgtttttaatggtgtgccTCGAAATGAAATTTGCGTCAATATCCTGACACAACGGAGAATGTGGCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTAc                                                                                                                                 <  2:567628/223‑1 (MQ=255)
                                                                                        aCTTGTTCTTGGTATTCTCCGACATGAAATCTTCCTCAAATTCCCGACACAACGCAGCATGGGTCACCAGCCATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATCAAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTAc                                                                                                                                 <  2:67523/239‑1 (MQ=255)
                                                                                             ttCATGGTAGTCCTCGACATGAAATCTTCTTCAAATTCCTGACAAAACGCAGCATGTGTCCCCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCg                                                                                                                            <  2:461814/239‑1 (MQ=255)
                                                                                             ttCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCg                                                                                                                            <  1:1275418/239‑1 (MQ=255)
                                                                                             ttCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACAAAACGCAGCATGGGGCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATcg                                                                                                                                                                           <  2:440274/193‑1 (MQ=255)
                                                                                                       tCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTTAATCG‑CCATTTATGATAACAAATTGATTTTTTTCTTGTTCTCCATATTATCGCGGTAATCAGCCAACCGGCTTTGCCCG‑TTAATTCCGGCCACCcgcgcccgctttc                                                                                                                   >  1:1065573/1‑229 (MQ=255)
                                                                                                       agacagACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACg                                                                                                                  <  1:118634/234‑1 (MQ=255)
                                                                                                        ccTAGAAATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACg                                                                                                                  <  1:940368/238‑1 (MQ=255)
                                                                                                         ccggACATGAAATCTGCGTCAATATCCTGCAACAACGCAGCATGGGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCa                                                                                                                <  1:1168044/236‑1 (MQ=255)
                                                                                                         cTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCa                                                                                                                <  1:516468/239‑1 (MQ=255)
                                                                                                              cATGAAATGTGGGTCAATACCCTGAAAAAACGAAGCATGTGTCACCAGCGATAAACCCGCCAGCAGAAAAAACGGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAg                                                                                                            <  1:685447/238‑1 (MQ=255)
                                                                                                                         cGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCg                                                                                                                                                                                                                           <  2:368250/118‑1 (MQ=255)
                                                                                                                         cGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCg                                                                                                                                                                                                                           >  1:368250/1‑118 (MQ=255)
                                                                                                                                      acacAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGGAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCCGATTACGCAGTAACGGGGATAGCC‑GAATATCCTCTTc                                                                                   >  1:316441/1‑238 (MQ=255)
                                                                                                                                            cGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTTGCCATATTATCCCGGTAATCAGCAACCCGGCTTGGCCCA‑TTAATTCCGGCTACCCGCGTCAGGTTACGGCGCAACGGGGAAAGCT‑GAATATCATCTTCCCAAAGt                                                                            >  1:1092644/1‑239 (MQ=255)
                                                                                                                                            cGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCA‑CCCTTGATGATAACAAATTGATTTGTGTCTGTTTCGCCCTATTATCGCGGGAATCAGCAACCCGGCTTGGGCAG‑TTAATTTCGGATACCAGCGTGAAATTACGCAGGAGCGGGGAAAGCT‑GAATATCATATTCCGcaaaa                                                                            >  1:152029/1‑237 (MQ=255)
                                                                                                                                            cGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGG‑CA‑GGA‑TTAAATCG‑CCATTTATGATAACAAATTTATTTGTGTCTGTTTTGCCATATTATCCCGGTAATCAGCAAAACGGCCTTGCCAGTTTAATTCC‑GCTACCACCGTCCAATTTACCAGTAGCGGGGAAAACTGGAATAT‑ATCTTCCCaaaaa                                                                            >  1:91282/1‑237 (MQ=255)
                                                                                                                                               agcaTGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTc                                                                          <  1:769775/238‑1 (MQ=255)
                                                                                                                                                 caTGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGc                                                                        <  1:122916/238‑1 (MQ=255)
                                                                                                                                                         aCCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGG‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAAATGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGGTTAAGCAGTAGCGGGAAAAGCT‑GAAAATCATCTTCCGAAAATTGCCCATTCCCgg                                                               >  1:1226679/1‑239 (MQ=255)
                                                                                                                                                             gCGATAAACCCGCCAACAGAAAAAACTGAAAACGGCGGCAACCCTTGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGtt                                                           <  1:1007384/239‑1 (MQ=255)
                                                                                                                                                                        ggcagcagAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGGGGCGGGT‑CA‑GGA‑TTAAAACG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGa                                                <  1:124050/237‑1 (MQ=255)
                                                                                                                                                                        cccAGCAGAAAAAACTGAAAACGGGGGCAACCCGCGAATACAGGCTGGCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGa                                                <  2:1171523/238‑1 (MQ=255)
                                                                                                                                                                        cccAGCAAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGa                                                <  1:225412/238‑1 (MQ=255)
                                                                                                                                                                             cagAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAAc                                                                                                                                                                                                            <  1:319117/80‑1 (MQ=255)
                                                                                                                                                                              agAAAAAACTGAAAACGGCGGCACCCCGGGAAATCAGGCTGCCGGGGGGGGT‑CC‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGGGTCTTTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGt                                           <  2:415585/238‑1 (MQ=255)
                                                                                                                                                                              aaaaaaaaCTGAAAACGGCGGCAACCCGGGAATACGGGCTGCCGCGGCGGGG‑CA‑GGA‑TTAAATCG‑CCATTGATTATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGt                                           <  2:1092644/236‑1 (MQ=255)
                                                                                                                                                                                     aCTTAAAACGGCGGCAACCCGGGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAg                                   <  1:1203779/239‑1 (MQ=255)
                                                                                                                                                                                     aCTGAAAACGGCGGCAACCCGGGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAg                                   <  1:707880/239‑1 (MQ=255)
                                                                                                                                                                                      cTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAgc                                  <  1:945020/239‑1 (MQ=255)
                                                                                                                                                                                       tGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTTAATCG‑CCATTGTTGATAAAAAATTGAttt                                                                                                                                                                                                  >  2:319117/1‑80 (MQ=255)
                                                                                                                                                                                                ccccacccccccattACCGGTTGGCGGGGGGGGT‑CA‑GGA‑TTAAATTG‑CCATTGATGATAACAAATTGTTTTGTGGTTGTTTCCGCAAATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTTCCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAat                        <  2:1158429/225‑1 (MQ=255)
                                                                                                                                                                                                    aaCCCGCGAAATAAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAatcat                     <  2:152029/238‑1 (MQ=255)
                                                                                                                                                                                                        ccctaatacgttgggtccgtgcggGT‑‑A‑GGA‑TTAAATCT‑CCATTGAGGATAACAAATTGATTTGTGTTTGTTGCGCTTGATTATCGCGGTAATCCGAGATGGTGCTTGGTCAG‑TTAAGTTCGGCTACCAACGTCAGATTTCGCAGTAGCGGGAATAGCT‑TAATATCATCTTCCGAAAGTTAGCTATTCACGGCGTTCGGTTTGACGATCAGTGTGTCCAGCGCACGTAAATGATCGCTGa               <  2:710080/219‑1 (MQ=255)
                                                                                                                                                                                                         gtcaatacgggtttgcgcgGCGGTT‑CA‑GGA‑TTAAATCG‑CCATTGATGAAAACAAATTGTTTTGTGTGTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGa               <  2:577481/224‑1 (MQ=255)
                                                                                                                                                                                                         gggAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGa               <  1:143073/237‑1 (MQ=255)
                                                                                                                                                                                                         cgaaaTAAAGGCGGCCGGGGGGGGT‑CA‑GGA‑TTAAATCC‑GCATTGATTCAAACAAATTGGTTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGa               <  2:1193359/236‑1 (MQ=255)
                                                                                                                                                                                                                    cttccGCGGGGGGT‑CA‑GTA‑TTTAATTG‑CCATTGGTTATAACAAATTGCGTTGTTTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAAt     <  2:470929/235‑1 (MQ=255)
                                                                                                                                                                                                                    cTGCCGGGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGGGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAAt     <  1:584205/238‑1 (MQ=255)
                                                                                                                                                                                                                    cTGCCCGGGCGGGT‑CA‑GGA‑TTAAAACG‑CCATTGATGATAACAAATTGATGTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAAt     <  1:1246041/238‑1 (MQ=255)
                                                                                                                                                                                                                      ggcgcgGCGGGT‑CA‑GGA‑ATAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCa   <  1:1335326/236‑1 (MQ=255)
                                                                                                                                                                                                                      gCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCa   <  1:474208/238‑1 (MQ=255)
                                                                                                                                                                                                                       ccGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAg  <  1:659353/238‑1 (MQ=255)
                                                                                                                                                                                                                                              aTCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGt                                                                            <  2:366693/144‑1 (MQ=255)
                                                                                                                                                                                                                                              aTCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTAGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGt                                                                            >  1:366693/1‑144 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                           
TTAATTTCTTGTTCGGTGATGGTGTATTGGGTGAGTTGATTACAGCCAACGAGCAGGCCACTGACGATCAATGCAGCGGCAAATAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGT‑CA‑GGA‑TTAAATCG‑CCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAG‑TTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCT‑GAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAG  >  NC_000913/1123176‑1123628

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