Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,132,787 A→C 100% I204S (ATC→AGC)  wcaA ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,132,7870AC98.1% 134.3 / ‑5.3 48I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (26/21);  minor base T (1/0);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGAC  >  NC_000913/2132577‑2133010
                                                                                                                                                                                                                  |                                                                                                                                                                                                                               
tCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCAGGAGAATTTTTTCGGAGACGAGGGGGTCTGCATTTCACCGTGATTGATATACAGGCCCTGCGGCGCCtcttct                                                                                                                                                                                                                >  2:159321/1‑228 (MQ=255)
                  ttaccgcccGAGAGGAGGGTCAGAAAAGTGCGCCAGGTCAGACGTTTAGTGCGGGATTGGGAGAGGGTAAACATGTGATATTTTTTGCTGGCGCGGTCGAATTTGTTTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTGCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                                                                                                                                                                              <  2:464653/221‑1 (MQ=255)
                               gcagcGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAATccaccacca                                                                                                                                                                                 >  1:98696/1‑228 (MQ=255)
                                    gTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGCATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGGTATGCAGGCTCTGCTTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGc                                                                                                                                                                            >  1:498661/1‑228 (MQ=255)
                                            aGTGCGCCCGGTCATACGTTTATTGCGGTTCTGGTATAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGa                                                                                                                                                                     <  1:159321/227‑1 (MQ=255)
                                                   cAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTAGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGGGATCTGCATTTCAGCGAGATTGATATGAAGGTCCTGCGTCGCCTCTTTTACTTTCCACGGTTCGCCGTCCTCACCCACCATCTCCAGGAAAATGTGGtca                                                                                                                                                             >  2:135292/1‑226 (MQ=255)
                                                          tACGTTTATTGCGGATCTGGTAGAGGGTACACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCCCAGGAAGATGTGGTAATCCtgtg                                                                                                                                                      >  2:657895/1‑228 (MQ=255)
                                                          tACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCCCAGGAAGATGTCGGAATCCttgg                                                                                                                                                      >  2:23741/1‑225 (MQ=255)
                                                                  ttGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                                                                                                               <  1:23741/227‑1 (MQ=255)
                                                                      ggagttGGTAGTGGGTAAACAGTTGATTTTTTTTGCTGGCGCGGTCGAATTTGGCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                                                                                                          <  2:236227/223‑1 (MQ=255)
                                                                      ggATCTGGTGGAGGGTAAAAAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTGTCTCGGTGACGAGGTGATCTGCATTGCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                                                                                                          <  1:252541/228‑1 (MQ=255)
                                                                        aTCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTAACCGTGATTGATGTGCAGGCTCTGGGTCGCCTGTTATACTTTCCACGGTTCCTCGTACTCCACCACCACCCTCAGTCAGATGACTGAATCCTGTGCCGCTTTCATTTCg                                                                                                                                        >  2:20851/1‑228 (MQ=255)
                                                                                  gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCCTTCAGTCCggaacccacca                                                                                                                              >  2:146507/1‑219 (MQ=255)
                                                                                  gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTTGCCCGCTTTCAGTTCGGCATCGaaca                                                                                                                              >  1:158080/1‑228 (MQ=255)
                                                                                     tAAACCGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                                                                                                            <  1:20851/227‑1 (MQ=255)
                                                                                     tAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCATGCTCTGCGTCGCCTCTTCTCCTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGACGAGGTCGCAATCCTGTGCCGCTTTCAGTTCGGTTTCGaacaaa                                                                                                                            >  2:358017/1‑227 (MQ=255)
                                                                                     tAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCtt                                                                                                                                                >  2:440559/1‑207 (MQ=255)
                                                                                      aaaCATGTGATTTTTTGTGCTGGCGCGGTCTAATTTGTTTTTGTGCTTGCGGTAAAAGGGGAAGTTCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaaca                                                                                                                          <  1:254517/228‑1 (MQ=255)
                                                                                      aaaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                                                                                                               <  1:440559/207‑1 (MQ=255)
                                                                                       aaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctct                                                                                                                                                                                                                   >  2:473803/1‑138 (MQ=255)
                                                                                        aCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                                                                                                                  <  1:473803/138‑1 (MQ=255)
                                                                                             tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCCCCGGAAGCTTTGCTAATCCTGTGCCGCTTTCGTTTCGGTTTCGAACAAAAACTtttt                                                                                                                    >  2:552073/1‑227 (MQ=255)
                                                                                               atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTCTCGAACAAACACCCTTTTaa                                                                                                                 >  2:553512/1‑228 (MQ=255)
                                                                                                     tttGCTGCCGCGGTCTAATTTGTTTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTGTCGGTGACGCGGTGATCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCa                                                                                                           <  2:174443/228‑1 (MQ=255)
                                                                                                        ggtgtcgcgGTCGGATTTGTCTTTGGTCGTGCGGTAAAATTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGcc                                                                                                            <  1:135292/219‑1 (MQ=255)
                                                                                                            gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCAGCTTTCAGTTCGGTATCGAGCAACAACTCTTTGAAACGCCAGGc                                                                                                        >  1:431221/1‑224 (MQ=255)
                                                                                                                     agtttGTCGTTGTGATTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGGCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAc                                                                                            <  1:587177/225‑1 (MQ=255)
                                                                                                                         tGTCTTTGTGCTTGCGCTAAAAGTGGAAGGACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAGCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTgg                                                                                        <  1:595104/227‑1 (MQ=255)
                                                                                                                           ttttttggCTTGGGGGAAAAGGGGAAGTCCCGGGATAATTTTTGCGGGGGCGAGGTGGTCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGtt                                                                                      <  2:609435/220‑1 (MQ=255)
                                                                                                                              ttGGGCTTGCTGTAAAAGTGTGACTCCCCGGAGAATTTTTTCGGTGACGAGGTGATCTTCATTTCACCGTGATTGCTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCTACGGTTCGCCGGTCTCCACCACCATCCTCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCa                                                                                  <  1:611152/228‑1 (MQ=255)
                                                                                                                                     tGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAGTTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACCCCCTTTGAAAAGCCATCCCCCGGAAAAGCCATCGTTGCCAAtaaaat                                                                           >  2:527136/1‑223 (MQ=255)
                                                                                                                                        ggTAAAAGTTGAATTTCCCGGGAAATTTTTTCGGTGACGAGGTGATCTGCATTTCCCCGTGATTGATAGCCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGc                                                                        <  2:117214/228‑1 (MQ=255)
                                                                                                                                        ggTAAAAGTTGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGc                                                                        <  1:553512/228‑1 (MQ=255)
                                                                                                                                                    aGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAACCGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACa                                                            >  2:500589/1‑228 (MQ=255)
                                                                                                                                                    aGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACa                                                            >  1:172692/1‑228 (MQ=255)
                                                                                                                                                        cccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCATTGATATTGCGTTTGTAGACCAggcg                                                        >  1:262420/1‑228 (MQ=255)
                                                                                                                                                                                ggTTATCTGCATTTCCCCTTGGTTGATATGCGGGCTCTGCGTCGCCTTTTCTACTTTCCACGTTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATa                                <  2:366939/228‑1 (MQ=255)
                                                                                                                                                                                    aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGACGTGAATACGGCGATTTGGGATACAGc                            >  1:171226/1‑228 (MQ=255)
                                                                                                                                                                                           tttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGc                     >  1:136688/1‑228 (MQ=255)
                                                                                                                                                                                            ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGc                                                                                                        >  2:24784/1‑144 (MQ=255)
                                                                                                                                                                                            ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGGAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGGGAATACGGGGCTTTGGGATACAGCGGCaaaca                    >  2:219658/1‑225 (MQ=255)
                                                                                                                                                                                             tCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGcc                                                                                                       <  1:24784/144‑1 (MQ=255)
                                                                                                                                                                                                 cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGGGAAACCGGCGATTTGGGAAACAGCGGAAAACCTGCCg               >  2:385665/1‑228 (MQ=255)
                                                                                                                                                                                                   tgattgatATGCAGGCTCTGGGGCGCCACTTCTTGTTTCCACGGTTCGCCGTACCCCCCCCCCATCCGCAGGAAGATGTCGTAAACCGGGCCCGCTTTCAGTTCGGTATCGAACAAACACGCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGt             <  2:54921/228‑1 (MQ=255)
                                                                                                                                                                                                   tgattgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACACCGGCAAGCTTGCCGGt             >  1:494828/1‑228 (MQ=255)
                                                                                                                                                                                                   tgattgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAACCACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGTTACACCGGCAAGCTTGCCg               >  1:322977/1‑226 (MQ=255)
                                                                                                                                                                                                   gttttgatATCCAGGCTCTGGGGCGCCTTTTCTACTTTCCCCGGTTCGCCGTACTCCACCCCCATCCGCAGGAAGAGTTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGt             <  2:158080/225‑1 (MQ=255)
                                                                                                                                                                                                              cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGAc  >  1:72751/1‑228 (MQ=255)
                                                                                                                                                                                                                  |                                                                                                                                                                                                                               
TCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGAC  >  NC_000913/2132577‑2133010

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