Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,987,499 C→A intergenic (+71/+8) ftnB → / ← yecJ putative ferritin‑like protein/DUF2766 domain‑containing protein YecJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,987,4990CA98.1% 162.3 / NA 53intergenic (+71/+8)ftnB/yecJputative ferritin‑like protein/DUF2766 domain‑containing protein YecJ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (26/25);  total (28/25)
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.00e+00

GAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGA‑ACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAA  >  NC_000913/1987252‑1987728
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
gAACTGTAGGATGATTCAACCGTCAATTGCCTGCGCGATCTGGAAAAAGAACTGCAGCATGATTGTCTGGTGCTGCAAACCATTCTTGATGATGTGCGCAGTGCGAAACTTGCGGGTAGGTGCTCTGTGCAGACCGACCAACA‑GTTGCTGAATGTCGTGTCACACCAGCTGCATTTATCATCATCGGCGCTAATGCATTGCGCCGATGAAGTGTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGt                                                                                                                                                                                                                                     <  2:412926/251‑1 (MQ=255)
        tgatgaTTCACCCGTCAATTTCCTGCGCGATCTGGCAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTa                                                                                                                                                                                                                              >  2:102864/1‑250 (MQ=255)
                   ccGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCAcc                                                                                                                                                                                                                  >  1:32712/1‑251 (MQ=255)
                    cGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCAccc                                                                                                                                                                                                                 <  1:520064/251‑1 (MQ=255)
                    cGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCAccc                                                                                                                                                                                                                 >  1:458336/1‑251 (MQ=255)
                                    gATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTg                                                                                                                                                                                                  >  2:294266/1‑250 (MQ=255)
                                    gATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTg                                                                                                                                                                                                  >  1:301103/1‑250 (MQ=255)
                                          gAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCaa                                                                                                                                                                                           <  2:32712/251‑1 (MQ=255)
                                          gAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCa                                                                                                                                                                                            >  1:194764/1‑250 (MQ=255)
                                            aaaaGAACAGCTGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATa                                                                                                                                                                                         <  2:137179/251‑1 (MQ=255)
                                             cacgAACAGCAGCATGATGGTCTGTTGCTGCACACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGGCCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCCTTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGGAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTagatcggaagggggtcg                                                                                                                                                                                        >  2:74842/4‑235 (MQ=255)
                                             aaaGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACCCCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGt                                                                                                                                                                                                         <  1:74842/234‑1 (MQ=255)
                                                 aaCAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCAcc                                                                                                                                                                                                                  <  1:457007/221‑1 (MQ=255)
                                                 aaCAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCAcc                                                                                                                                                                                                                  >  2:457007/1‑221 (MQ=255)
                                                      cagcaTGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATgg                                                                                                                                                                               >  1:351732/1‑251 (MQ=255)
                                                         caTGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGCAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTaa                                                                                                                                                                                                       >  2:684014/1‑224 (MQ=255)
                                                         caTGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTaa                                                                                                                                                                                                       <  1:684014/224‑1 (MQ=255)
                                                               ggTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTACTGGGATTCACCCTGTGGGGTAAACTTGCGTTCAGTAATTGGGATGGCTTTTTggt                                                                                                                                                                      >  2:19283/1‑250 (MQ=255)
                                                                    gTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTGCAATAAGCGCGATGGCTTTTTGGATTGcc                                                                                                                                                                 <  1:550606/251‑1 (MQ=255)
                                                                      tgctgcAAACCATTCTTGATGAAGTGCGCAGTGCGACACTTGCGGGTATGTGCCCTGTGCAGACCGACCACCATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGCTTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCg                                                                                                                                                               >  2:281610/1‑251 (MQ=255)
                                                                               ccATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGAc                                                                                                                                                      >  1:162820/1‑251 (MQ=255)
                                                                               ccATTCTTGATGAAGGGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGAc                                                                                                                                                      <  2:162820/251‑1 (MQ=255)
                                                                                     ttgttggagagcgcAGTGCGAACCTTACAGGGATGTACCATGTACAGACCGACAAAACTGTT‑CTGAATGTCGTGTCACACCAGCAGCATTGATCATCATCGACGCTAATGCATTGCGCCGATGAAGGATTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGGGGATTTTTTACTGTGATTAACACTGAGGGGTCACCTTGCGGTCAATAAGAGCGATGGCTTTTTGGATTGCCCGCTTGGTGATCggatc                                                                                                                                                >  2:67497/11‑248 (MQ=255)
                                                                                                 gcAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTGTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCggg                                                                                                                                    <  1:281610/251‑1 (MQ=255)
                                                                                                        gAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAgcg                                                                                                                                                                                     >  2:47424/1‑195 (MQ=255)
                                                                                                        gAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAgcg                                                                                                                                                                                     <  1:47424/195‑1 (MQ=255)
                                                                                                        gAAACTTGCGGGTATGTGACCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCAGCGGCGCTAATGCAGTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTa                                                                                                                             <  2:194764/251‑1 (MQ=255)
                                                                                                               gCGGGTATTGGCCCTGGGCAGATCGACCAACAGGTG‑CTGAATTGCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTGTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTGTTGGATTGCCCGCATGGTGACCGGTTCTGCGGCGGCGGGATGGTTAGTAAAg                                                                                                                       <  2:444316/250‑1 (MQ=255)
                                                                                                                 gggTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGCTTGCCCGCATGGTGACCGGGTc                                                                                                                                                >  1:118758/1‑223 (MQ=255)
                                                                                                                 gggTATGTGCCCTGTGCAGACCGACCAACATGGT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAAGAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTc                                                                                                                                                <  2:118758/223‑1 (MQ=255)
                                                                                                                              gtgCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCt                                                                                                       <  1:581237/251‑1 (MQ=255)
                                                                                                                                                gTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCTGCGCTAAGGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCa                                                                                      <  1:520579/250‑1 (MQ=255)
                                                                                                                                                      gAATGGCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACt                                                                                <  2:351732/251‑1 (MQ=255)
                                                                                                                                                                 cacaCCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGgcg                                                                      <  1:5595/250‑1 (MQ=255)
                                                                                                                                                                    acCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGAt                                                                  <  1:294266/251‑1 (MQ=255)
                                                                                                                                                                     ccAGCTGCTTTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGAGGGATGGTTAGTAAAGTCGATATTCTTCCGCTGACTGGGCATTTTTTCACGAAATTCAACGGGCGGGAtt                                                                 >  2:632707/1‑251 (MQ=255)
                                                                                                                                                                     ccAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGAtt                                                                 >  1:343377/1‑251 (MQ=255)
                                                                                                                                                                             aTTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCATCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGATTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGTCGGCGGGATGGTTAGTAAAGTCGATATTCTTAAGCTGACTGGAAATTTTTTCACGAAATTCAAAGGGCGCGATTACATCg                                                           >  1:45359/1‑249 (MQ=255)
                                                                                                                                                                                gATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGa‑a                                                      >  2:668918/1‑250 (MQ=255)
                                                                                                                                                                                 atcatcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGCTTACATCGAGA‑Ac                                                     >  2:274594/1‑250 (MQ=255)
                                                                                                                                                                                 atcagcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCGCCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTGCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGA‑Ac                                                     <  1:663791/250‑1 (MQ=255)
                                                                                                                                                                                    atcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACAGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAACGTCGATATTCTTCAGATGACTGGAAATTTTTTCACGCACTTCAACGGGCGAGATTAAATCGAGA‑ACATcc                                                 >  2:416467/1‑251 (MQ=255)
                                                                                                                                                                                                aaTGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGAAATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGTGA‑ACCTCCAGAAtttgtttg                                     >  2:202322/1‑251 (MQ=255)
                                                                                                                                                                                                 aTGCATTGCGCCGCTGAAGGTTTTGGGAACCCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTAACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGAACGCATGGAGACCAGGTCTGAGGCGGCGGGATGGTTATGAAAGACTATATTCTTCTGATGACTGGACATTTTTTCAAGACATGTAACGGGCGAGCTTACATCG‑GATACATCCCGAAtttgtttg                                     >  2:336352/1‑250 (MQ=255)
                                                                                                                                                                                                 aTGCATTGCGCCGATGAAGGTTTTGAGAACCCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGA‑ACATCCAGAAtttgtttg                                     >  1:519798/1‑250 (MQ=255)
                                                                                                                                                                                                    cATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGA‑ACATCCAGAATTTGTTTGATaa                                 <  1:102864/251‑1 (MQ=255)
                                                                                                                                                                                                       agagcCGATGAAGGGTTTGAGAAACCGCTGCCTGATCTGTTTGAAGCAGAGGTTTTTGTAATGGGATGCACACGGTGGGGTAAACTGGAGTTCAATAAGCGCGATGGCTTTGTGGATTGCCCGGATGTTGACCGGGTCTGCGGCGGCGGGAGGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTGTTGCACGATCTTCAACGGGCTCGATTACATCGAGA‑ACATCCAGAATTTGTTTGATAACCa                              <  1:300367/248‑1 (MQ=255)
                                                                                                                                                                                                        gcgtcGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGA‑ACATCCAGAATTTGTTTGATAACCa                              <  1:231560/246‑1 (MQ=255)
                                                                                                                                                                                                          gcCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGA‑ACATCCAGAATTTGTTTGATAACCAACt                           >  1:328258/1‑251 (MQ=255)
                                                                                                                                                                                                                         gagaAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAACGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGCTTACATCGAGA‑ACATCCAGAATTTATTTGATAACCACCTGGCAAGCCACcaca             >  2:163854/1‑250 (MQ=255)
                                                                                                                                                                                                                                    gCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCGCCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGTACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGA‑ACAGCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCaa  <  2:286798/250‑1 (MQ=255)
                                                                                                                                                                                                                                            tGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTa                                                                                                                             >  2:363967/1‑120 (MQ=255)
                                                                                                                                                                                                                                            tGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTa                                                                                                                             <  1:363967/120‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
GAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTT‑CTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGA‑ACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAA  >  NC_000913/1987252‑1987728

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