Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,939 G→A F7F (TTC→TTT yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,9390GA89.8% 141.1 / 7.0 59F7F (TTC→TTTyaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base A (23/30);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.05e-01

GTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCA  >  NC_000913/7695‑8183
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
gtcacccGGGGGTAATTGCCGGCGCACGGCCGGCCAGGTGGCAGCTCCCCCCCCGCGCTGCAAGACGGGTACACAATGAGGTAAATGAGTGTAACCCGCCGGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGTGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGgag                                                                                                                                                                                                                                                  <  2:283527/246‑1 (MQ=255)
         cgctactggCCAGCGCAAGGCCTGCCAGGTGGCCTCGACCCCAGCGCGCCGAAAGACTGGTAAACAATTACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTGTAAGACTTGTGCCAAACTGGCGGATGTAGCGAAACTGCACAAAACCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGGAAATCAGTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAGAAAATCt                                                                                                                                                                                                                                         <  2:274353/243‑1 (MQ=255)
            tAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGGGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGGGGAAAAATCTggg                                                                                                                                                                                                                                      >  2:663552/1‑250 (MQ=255)
             aattgCCAGCGCACCGCCGGCCAGGTGGCCGCTACCCACGCGCTCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGCCTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGgca                                                                                                                                                                                                                                     <  1:376219/247‑1 (MQ=255)
             aatggccaggcaaaaaccGGCAAGTTTGCCCCTACCCACGCGGGCCTCAAGATTGTTACACAATTACGTAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTATGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGgca                                                                                                                                                                                                                                     <  2:67150/235‑1 (MQ=255)
                   cAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCata                                                                                                                                                                                                                               <  1:508498/251‑1 (MQ=255)
                       gcAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATccc                                                                                                                                                                                                                           <  1:663552/251‑1 (MQ=255)
                                   aGGTTGCCGATACCCAAGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGCCTTTTGCCAAACTGGCGGATGGAGCGAAACTGCACAAATACGGTGCGAAAAGTGAACCAACAACCGGCGCAGCAGAGCAGGGAAATCATTACCGATACCCAAATGACGCTTGTAGTAAAGGAGAAAAAATCTGGCCAGCATATCCCTCTTATTGCCgg                                                                                                                                                                                                              >  2:519399/1‑251 (MQ=255)
                                            cTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATCAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGGGCGAAAAGTGAACCAACAACCTGGGCCGAAGAGCAGGGAAATCATTACCGATCCCCAAAGGACGCTGGTAATAAAGAAGAAAAAATCTTGCATGCTCATCCCTCTTATTGCCggcaccgatg                                                                                                                                                                                                      >  2:103408/1‑242 (MQ=255)
                                               cccACGCGCGCCGCAAGACTGTTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACtt                                                                                                                                                                                                 <  1:657629/251‑1 (MQ=255)
                                                ccACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACttt                                                                                                                                                                                                <  1:120823/251‑1 (MQ=255)
                                                     cgcgcCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGGGAACCAACAACCTGCGCCGGAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTACtgt                                                                                                                                                                                           >  2:266422/1‑251 (MQ=255)
                                                         ccgcAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTaa                                                                                                                                                                                       <  1:20626/251‑1 (MQ=255)
                                                            cAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGATGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGGTGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGCTCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACg                                                                                                                                                                                    >  1:262154/1‑251 (MQ=255)
                                                              aGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATTGTACGAAACTGCACAAATCCGGTGCGAAAAGGGAACCAACAAACTGCGCCGAAGAACAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAAAAAAGGAGGAAAAATCTGGCATGCAAATCCCTCTTATTTCCGGTCGCGATACATTTCCTGTGTAAAAGtt                                                                                                                                                                                  >  2:475894/1‑251 (MQ=255)
                                                                    gTACACAATGACGGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAAt                                                                                                                                                                            <  1:85672/251‑1 (MQ=255)
                                                                    gTACACAATGAATGAAAGGGGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAAt                                                                                                                                                                            <  2:522018/251‑1 (MQ=255)
                                                                     tACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGga                                                                                                                                                                                                                                                   >  2:220793/1‑181 (MQ=255)
                                                                     tACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGga                                                                                                                                                                                                                                                   <  1:220793/181‑1 (MQ=255)
                                                                        acaATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGtt                                                                                                                                                                        <  1:199050/251‑1 (MQ=255)
                                                                        acaATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGtt                                                                                                                                                                        <  1:516741/251‑1 (MQ=255)
                                                                         caagtaCTGAAAGTAGGTTAAACCGCATGGCTGTGTATGAATGCTATTTTTAAGAATTTTGCAAAACTGGCGGATGTAGCGAAACTGCAAAAATCCGGTGCGAAAAGTGAACCAACAGCCTGCGCCGAAGAGCAGGTAAATCATGACAGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAAtttttt                                                                                                                                                                       <  2:364508/246‑1 (MQ=255)
                                                                         caATGACTGAAAGGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGttt                                                                                                                                                                       <  1:195973/251‑1 (MQ=255)
                                                                           aTGACTGAAATGATGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaa                                                                                                                                                                     <  1:103408/251‑1 (MQ=255)
                                                                                        ggTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAACAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCGTATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAAAGTTACCAATTGTTTAAGAAGTATCTACGc                                                                                                                                                        >  1:646635/1‑251 (MQ=255)
                                                                                            aaaCCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGAGCGCAAAGTGAACCAACAACATGCGCCAAAGAGCAGGTCAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCAGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACg                                                                                                                                                    >  1:413949/1‑251 (MQ=255)
                                                                                                cGCCTGGCGGTGGATTAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTCGCGAAACTGTACAAATTCCGTGCGATAATTGACTCAACTACCTGCGCCGAAGCGCAGGTAAATCATTCCCGATCTCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTTATTCTGGTCTCGATGACCTTCCGGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGt                                                                                                                                                <  2:501607/251‑1 (MQ=255)
                                                                                                      gagggggATGAATGCTATTTTTAAGACGTTTGCCAATCTGGCGGATGTAGCGAAACTGCTCAAATCCGTTGCGAAAAGTGAACCAACAACCTGCGTCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTGAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGACGACTTTACTGTGTAAACGTGACCAATTGTTTAAGAAGTATATCCGCGACGATGTACTTGa                                                                                                                                          <  2:227743/246‑1 (MQ=255)
                                                                                                        tgtgGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGAt                                                                                                                                         >  1:337071/1‑250 (MQ=255)
                                                                                                         gtgGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATa                                                                                                                                        <  1:526969/250‑1 (MQ=255)
                                                                                                                aaTGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGc                                                                                                                                >  1:482912/1‑251 (MQ=255)
                                                                                                                                cTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGt                                                                                                                >  1:670600/1‑251 (MQ=255)
                                                                                                                                 ttttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGGGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGCACTTGATAACTTCTGCGTAGCATACATGAGGtt                                                                                                               >  1:170649/1‑251 (MQ=255)
                                                                                                                                     gCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGCAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACCTGCGGATTTGt                                                                                                           >  2:444413/1‑251 (MQ=255)
                                                                                                                                       cAAACTGGCGGATGTAGCGAAACTGCCCAAATCCGGTGCGAAAAGTGAGCCAACAACCTGCGCCGAAGGGCAGGTAAATCATTACCGATCCCCAAAAGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGCTTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGCGGTACTTGATAACTTCTACGAAGCATAcctgcgggtttggat                                                                                                         >  1:209078/1‑237 (MQ=255)
                                                                                                                                                 gATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATggcg                                                                                               >  1:464553/1‑251 (MQ=255)
                                                                                                                                                     gTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAGCAGCCTGCGCCTAAGAGCAGGTAAAACATTACCGAGCCCCAAAGGACGCTGTTAAGAAAGTAGAAAAAATCGGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTGCCAATTGGTTAAGAAGTATATACGCTACGAGTTACTTGATAACGTCTGCGTAGCATACATGAGGTTTTGTATATAAATggcgggc                                                                                            <  2:617306/251‑1 (MQ=255)
                                                                                                                                                       aGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGATTTGTATAAAAATggcgggcgg                                                                                          >  2:307438/1‑250 (MQ=255)
                                                                                                                                                           aaaCTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATc                                                                                      <  1:676572/251‑1 (MQ=255)
                                                                                                                                                              cTGCACAAATACGGAGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCACAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGGGatataaa                                                                                    >  2:428135/1‑247 (MQ=255)
                                                                                                                                                                 cacaAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGAACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATCCATGAGGTTTTGTATACAAATGGCGGGCGGTATccacac                                                                                 >  2:263319/1‑245 (MQ=255)
                                                                                                                                                                    aaaTCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAgtg                                                                             <  1:415830/251‑1 (MQ=255)
                                                                                                                                                                         cGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGa                                                                        <  2:464553/251‑1 (MQ=255)
                                                                                                                                                                         cGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTCCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAACATACATGAGGTTTTGTATAAAAATGGCGGGGGATGTAAACGCAgtgaca                                                                          >  2:173988/1‑246 (MQ=255)
                                                                                                                                                                           gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGCGAACACATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGAACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCCACGCCGTGTCAgga                                                                       >  1:273343/1‑248 (MQ=255)
                                                                                                                                                                              cGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGt                                                                                                                                                                <  1:488088/158‑1 (MQ=255)
                                                                                                                                                                              cGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGt                                                                                                                                                                >  2:488088/1‑158 (MQ=255)
                                                                                                                                                                                  aaGTGAACCTACAACCTGCGCCGAAGAGCAGGTAAATCATTAAAGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGAATACATGAGGTTTTGTATAAAAATGGCGGGGGATATCAACGCaaatgcacaaattcgaa                                                                >  2:444666/1‑234 (MQ=255)
                                                                                                                                                                                  aaGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGaaa                                                               <  1:468907/251‑1 (MQ=255)
                                                                                                                                                                                         ccAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGAAACAGtctc                                                        <  1:493752/251‑1 (MQ=255)
                                                                                                                                                                                                          agagCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGtt                                                                                                                                                                                  <  1:517324/112‑1 (MQ=255)
                                                                                                                                                                                                          agagCAAGAAAATCCTTACCGATCCCCACAGGACTCTGTTAATAAAGGAGGAAAAATAGGGCATGCATATCCATCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGtt                                                                                                                                                                                  >  2:517324/1‑112 (MQ=255)
                                                                                                                                                                                                               aGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTc                                  <  1:225550/251‑1 (MQ=255)
                                                                                                                                                                                                                      tCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGaa                                                                <  2:128125/214‑1 (MQ=255)
                                                                                                                                                                                                                      tCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGaa                                                                >  1:128125/1‑214 (MQ=255)
                                                                                                                                                                                                                          tACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATCGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCATAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTTAGAAATCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTcggcgggtgtgc                       >  1:322110/1‑245 (MQ=255)
                                                                                                                                                                                                                          tACCGATCCCCAAAGGACGCTGTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGGTgcgc                       >  1:569366/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aGGACGCTGTTAATAAAGGAGAAAAATTCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGcgtcg          <  1:49004/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gTTAATAAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCGGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATatca  <  1:330703/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
GTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTGCGTAGCATACATGAGGTTTTGTATAAAAATGGCGGGCGATATCAACGCAGTGTCAGAAATCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCA  >  NC_000913/7695‑8183

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