Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,451,083 Δ1 bp 10.8% coding (25/1503 nt) ytfR → putative sugar ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,451,0830A.10.8% 117.8 / 7.5 37coding (25/1503 nt)ytfRputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (16/17);  new base . (2/2);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-01

GCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCAATATGTCGGTCGCTGAT  >  NC_000913/4450797‑4451367
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                            
gcgcTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATGGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGAtc                                                                                                                                                                                                                                                                                            <  1:217310/289‑1 (MQ=255)
                                                  aaCCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGa                                                                                                                                                                                                                                          <  2:220137/289‑1 (MQ=255)
                                                                        cctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTgcg                                                                                                                                                                                                                     >  2:12226/1‑288 (MQ=255)
                                                                          tgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCg                                                                                                                                                                                                                  >  1:151567/1‑289 (MQ=255)
                                                                            aTACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTg                                                                                                                                                                                                                >  2:191174/1‑289 (MQ=255)
                                                                                cGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGa                                                                                                                                                                                                            <  2:151567/289‑1 (MQ=255)
                                                                                          gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCa                                                                                                                                                                                                                                                    <  2:243912/239‑1 (MQ=255)
                                                                                          gaaTTAGAGAAGAAGAAAAATATGGGGAATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCa                                                                                                                                                                                                                                                    >  1:243912/1‑239 (MQ=255)
                                                                                            aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGgcg                                                                                                                                                                                                  <  2:56482/287‑1 (MQ=255)
                                                                                                          aaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGaaaa                                                                                                                                                                                    <  2:74537/287‑1 (MQ=255)
                                                                                                                   ggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCCGCCTGCGCCGTGGCGAAATCATGGCCCTGCGCGGTGAAAACGGGGCGGGaa                                                                                                                                                                         >  2:188950/1‑289 (MQ=255)
                                                                                                                            ggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCa                                                                                                                                                                                                                                                    >  1:145642/1‑205 (MQ=255)
                                                                                                                            ggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCa                                                                                                                                                                                                                                                    <  2:145642/205‑1 (MQ=255)
                                                                                                                               tGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTaa                                                                                                                                                             >  1:263716/1‑289 (MQ=255)
                                                                                                                                 cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCCCTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGACAACGGGGCGGGAAAATCAACGCTAATc                                                                                                                                                           >  1:309779/1‑289 (MQ=255)
                                                                                                                                          aTGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTa                                                                                                                                                  >  2:181392/1‑289 (MQ=255)
                                                                                                                                           tGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTaa                                                                                                                                                 <  1:229879/289‑1 (MQ=255)
                                                                                                                                           tGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTa                                                                                                                                                  >  2:229879/1‑288 (MQ=255)
                                                                                                                                            gCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTaa                                                                                                                                                 <  1:124024/288‑1 (MQ=255)
                                                                                                                                            gCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTaa                                                                                                                                                 <  2:263716/288‑1 (MQ=255)
                                                                                                                                             ccGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAAc                                                                                                                                                <  1:12226/288‑1 (MQ=255)
                                                                                                                                                            gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTAAGCCTTCCCCGTGGCGAAATCATGGCGCTGCCCGGTGAAACCGGGGCGGGAAAATCCACGCAAATCAAAGCATTAACTGgtgtcacacacgtc                                                                                                                                >  2:171037/1‑279 (MQ=255)
                                                                                                                                                                     cTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGc                                                                                                                       <  1:191174/289‑1 (MQ=255)
                                                                                                                                                                      tCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGc                                                                                                                       >  1:55680/1‑288 (MQ=255)
                                                                                                                                                                                                  aaGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGTCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTaa                                                                                                                                                 >  1:33752/1‑233 (MQ=255)
                                                                                                                                                                                                  aaGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGTCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTaa                                                                                                                                                 <  2:33752/233‑1 (MQ=255)
                                                                                                                                                                                                               cGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGTCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGt                                                                                                                                                                                                                                    >  1:103186/1‑137 (MQ=255)
                                                                                                                                                                                                               cGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGTCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGt                                                                                                                                                                                                                                    <  2:103186/137‑1 (MQ=255)
                                                                                                                                                                                                                                 aTGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTc                                                           <  2:309779/289‑1 (MQ=255)
                                                                                                                                                                                                                                            gAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCa                                                 <  2:55680/288‑1 (MQ=255)
                                                                                                                                                                                                                                              gCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCa                                                                                                                                                          <  2:165810/181‑1 (MQ=255)
                                                                                                                                                                                                                                              gCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCa                                                                                                                                                          >  1:165810/1‑181 (MQ=255)
                                                                                                                                                                                                                                              gCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCg                                              >  1:171247/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                  aaCCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTg                          <  2:171247/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                              cgaccAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCaa                  >  1:177641/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                                   aaCACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCAATATgtcggt         <  1:181392/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          ggAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGACCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGCCACCGTCTATCATGAAGTCAACCTTCTACCCAATATGTCGGTCGCTGAt  >  1:66354/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                            agaTCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCAATATGTCGGTCGCTGAt  >  2:34871/1‑287 (MQ=255)
                                                                                                                                                                                                                                                                                            agaTCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCAATATGTCGGTCGCTGAt  <  1:34871/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                            
GCGCTGGAGAAATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCAATATGTCGGTCGCTGAT  >  NC_000913/4450797‑4451367

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