Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,356,852 A→T W87R (TGG→AGG)  sapA ← antimicrobial peptide transport ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,8520AT97.6% 145.1 / NA 42W87R (TGG→AGG) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (22/19);  minor base G (0/1);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAA  >  NC_000913/1356561‑1357132
                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                        
cGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGc                                                                                                                                                                                                                                                                                 <  3:49443/301‑1 (MQ=255)
                 ggtgtaTCCGGCGGGGCGAGGCGGAACTCAACGGTATGAGTAACCAGGTTGCGGACGCTTTGCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGc                                                                                                                                                                                                                                                                <  1:279048/296‑1 (MQ=255)
                      aTCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGc                                                                                                                                                                                                                                                                >  4:279051/1‑296 (MQ=255)
                                             cAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCcagcag                                                                                                                                                                                                                                    >  3:86039/1‑301 (MQ=255)
                                             cAACGCTCTGTTTATCCACGATGCGCACGCTTTTCACGTTATCTGGTAATCGCAGGATGTCGAAGTAGGGGAAGTTGCTGCCATTGACGGTATGCCACGGGTTGATGCGGTCCAACGTGCGCTGGAATGTAAAAACCACATCGTCGGCACTCATTTGACGAGTGGGAGTAAACCAGGCGGTTTGTGAAAACGGAACATCGCGAGGCAGGTGGAAGCGATAGGTCGCGGCGTTGGCGAGTACTTCCATGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGAAGAGACATCcagcag                                                                                                                                                                                                                                    <  1:86661/301‑1 (MQ=255)
                                                        cgtgtgctcttccgatctcTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATaaa                                                                                                                                                                                                                         <  3:10293/283‑1 (MQ=255)
                                                         tATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATacaa                                                                                                                                                                                                                        >  3:69449/1‑301 (MQ=255)
                                                              aGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAATCAGTCGGTTTTTTGAAACGGAACATCGAGACGCAGGTGGAATCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCAGACATCAgcc                                                                                                                                                                                                                                                                >  1:104087/1‑254 (MQ=255)
                                                              aGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGc                                                                                                                                                                                                                                                                <  4:104088/256‑1 (MQ=255)
                                                                        cGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAggg                                                                                                                                                                                                          <  2:85943/300‑1 (MQ=255)
                                                                          cTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATaaaagatcggaagagcgtcgg                                                                                                                                                                                                       >  2:10278/1‑284 (MQ=255)
                                                                           ttttCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATc                                                                                                                                                                                                      >  1:171027/1‑301 (MQ=255)
                                                                                  ggtATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGGTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACGTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAAtt                                                                                                                                                                                               <  1:16/299‑1 (MQ=255)
                                                                                                         cGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATgg                                                                                                                                                                        >  3:111031/1‑301 (MQ=255)
                                                                                                                                             aCGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGAc                                                                                                                                    <  4:171028/301‑1 (MQ=255)
                                                                                                                                                 gttgttGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACaaa                                                                                                                                 >  4:329472/1‑300 (MQ=255)
                                                                                                                                                    tttgCGGTAAAAAAAGCGCGGGAAAGTAAACACCACAGCGGCGGCATTCATTTTACGCGTGGGAGTAAACCAGTAGGTTTGTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAATACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCg                                                                                                                             <  2:69368/300‑1 (MQ=255)
                                                                                                                                                          gTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGcc                                                                                                                                                    <  4:152811/272‑1 (MQ=255)
                                                                                                                                                          gTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGcc                                                                                                                                                    >  1:152810/1‑272 (MQ=255)
                                                                                                                                                          gTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTcg                                                                                                                       >  4:50862/1‑301 (MQ=255)
                                                                                                                                                             aaaaaTTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATggg                                                                                                                                                                       <  1:119261/250‑1 (MQ=255)
                                                                                                                                                             aaaaaTTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATggg                                                                                                                                                                       >  4:119262/1‑250 (MQ=255)
                                                                                                                                                             aaaaaTTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCgg                                                                                                                    >  4:183037/1‑301 (MQ=255)
                                                                                                                                                                            ggTAAACACCACATCGTCGGCATTCATTTTACGGGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCggg                                                                                                     <  1:183036/301‑1 (MQ=255)
                                                                                                                                                                             gTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGACCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCTCGCTGACGTAATCGACAAAACCGCTGTCGCGGATATCAGCATGCgggg                                                                                                    >  4:151211/1‑301 (MQ=255)
                                                                                                                                                                                 acacCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGAc                                                                                                                                                                                                                                             <  1:50405/160‑1 (MQ=255)
                                                                                                                                                                                 acacCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGAc                                                                                                                                                                                                                                             >  4:50405/1‑160 (MQ=255)
                                                                                                                                                                                                              gTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGCTGGGTTAAAGGTGTTGACTTGCCCGCTGCCGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTca                                                                   >  3:34728/1‑301 (MQ=255)
                                                                                                                                                                                                               tGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGc                                                                                                                                           <  4:251349/228‑1 (MQ=255)
                                                                                                                                                                                                               tGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGc                                                                                                                                           >  1:251346/1‑228 (MQ=255)
                                                                                                                                                                                                                gggAGTAAGCCAGTCGGTTTGTTGAAACGGAACATCGGTAAGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTGCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAAGTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTcac                                                                  <  2:110896/300‑1 (MQ=255)
                                                                                                                                                                                                                                  tttttGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCaa                                                >  1:257695/1‑300 (MQ=255)
                                                                                                                                                                                                                                       gAAACAGAACAACGCGACGCAGGTAGAAGCGATAGGTCGAGCCGTTGTTGAGTACTTCCCTGCTGTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAagag                                           <  1:197498/300‑1 (MQ=255)
                                                                                                                                                                                                                                            ggAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGata                                     >  3:77697/1‑301 (MQ=255)
                                                                                                                                                                                                                                                   cgcgACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTg                              >  1:98573/1‑301 (MQ=255)
                                                                                                                                                                                                                                                      gACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCGGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGgcg                           <  1:50862/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGAcc                                                                        >  4:241920/1‑242 (MQ=255)
                                                                                                                                                                                                                                                                    cGATAGGTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGAcc                                                                        <  1:241917/242‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ggTCGCGCCGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGAtttt        <  4:98574/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                ccGTTGTCGAGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGCTAGGTATAGGGATCGACATCCAGCAGCCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTCACCCACTGCTCGCTTTGGACTGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAACACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGCTTGCCTGCCCACTCACAAGTTCAGCAATCACCAAAAGAGACGGTAATACCTTGCCCATAATCATAAGGGATTTTATGTaa  >  1:12627/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                         aGTACTTCCCTGCTTTTGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATc                                                                                                               >  4:36229/1‑182 (MQ=255)
                                                                                                                                                                                                                                                                                         aGTACTTCCCTGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATc                                                                                                               <  1:36229/182‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                        
CGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAA  >  NC_000913/1356561‑1357132

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