Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,629,921 C→T E228K (GAG→AAG)  yhiI ← putative membrane fusion protein YhiI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,629,9210CT97.9% 168.8 / NA 48E228K (GAG→AAG) yhiIputative membrane fusion protein YhiI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/25);  minor base G (0/1);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGAACATCAGTTTCAGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTCGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGC  >  NC_000913/3629675‑3630164
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
ggAACATCAGTTTCAGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGc                                                                                                                                                                                                                                                    >  1:357165/1‑250 (MQ=255)
       cagttccggCCGTTCATCGCTGTTTTCGACGGTTTTGGGCGTGAACGGGGCGACCCTGGCGAAAAAACGGATGGTGGCAGGAATCCGCAGATCGGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGAt                                                                                                                                                                                                                                            <  2:70674/243‑1 (MQ=255)
                ccGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATCTTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCAcgc                                                                                                                                                                                                                                   >  2:319505/1‑251 (MQ=255)
                     cATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGcgtccg                                                                                                                                                                                                                              >  1:15601/1‑251 (MQ=255)
                             gTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGATTACGCAGATCTGGCGCGGCATCGAGGATCAGCAGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCGGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCgg                                                                                                                                                                                                                      >  1:179091/1‑251 (MQ=255)
                                       tttttGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGt                                                                                                                                                                                                            >  1:141484/1‑251 (MQ=255)
                                              cGTGAACTGGGCGACACTGGCGACAACACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGCCTCCACCGCCCAGTTTCAGAGTGCCCGCCTGCTCGCTTGGCAGGATGAAAGTCATCTAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGt                                                                                                                                                                                                     >  1:142827/1‑251 (MQ=255)
                                                 gAACTGGGCGACACTGGCGACAAATCTGATGTTTGCAGGAATACGCAGATCTGGCGGGGCATCGAGGATCATCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATATACGTCGCTGAGATCGACCATATTCAGCTCCCGACCGCCTGCCGCCAGTACTTCGCCTTGCTTGGCAACCCGATACTGCACTCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCAt                                                                                                                                                                                                  <  1:47744/251‑1 (MQ=255)
                                                 gAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCAt                                                                                                                                                                                                  >  2:267460/1‑251 (MQ=255)
                                                  aaCTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATc                                                                                                                                                                                                 <  1:267413/251‑1 (MQ=255)
                                                     ggggCGACCCTGGCGACAAACCTGATGGTTGCAGAAATACGCAGATCGGGCGCGGCATCGGGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGat                                                                                                                                                                                              <  2:93848/250‑1 (MQ=255)
                                                       ggCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGatat                                                                                                                                                                                            >  2:2892/1‑251 (MQ=255)
                                                               tGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGc                                                                                                                                                                                     >  2:305776/1‑250 (MQ=255)
                                                               gggCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCa                                                                                                                                                                                    <  1:157650/250‑1 (MQ=255)
                                                                  cGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTGCCCCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATg                                                                                                                                                                                 <  1:352879/251‑1 (MQ=255)
                                                                           tGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTc                                                                                                                                                                         >  1:203162/1‑250 (MQ=255)
                                                                           gggtggtgggaGGAATACGCAGATTTGGGGCGGCAGCGGGGATCAGCGGGGCTTCACCGCCCAGTTCCAGCGGGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCa                                                                                                                                                                        <  2:84920/241‑1 (MQ=255)
                                                                            gATGGTTGCAGGCATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATCGACGTCGCTGAGATCGACCCTATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTc                                                                                                                                                                         >  2:221365/1‑249 (MQ=255)
                                                                                    aaGGAATACGCAGATCTGCCGCGGCATCGAGTATCAGCCGGGCTTAACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTtgt                                                                                                                                                               <  1:165865/250‑1 (MQ=255)
                                                                                     aGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCCCCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCaa                                                                                                                                                                                   <  1:94348/230‑1 (MQ=255)
                                                                                     aGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCaa                                                                                                                                                                                   >  2:94348/1‑230 (MQ=255)
                                                                                      ggAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTTTGc                                                                                                                                                             >  2:271881/1‑251 (MQ=255)
                                                                                       gAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGcc                                                                                                                                                            <  1:48018/251‑1 (MQ=255)
                                                                                         aTACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGc                                                                                                                                                          >  1:24453/1‑251 (MQ=255)
                                                                                           ccGCAGATCTGGGGCGGCACCGGGGATCAGCCGGGCTTCCCCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCtt                                                                                                                                                        <  2:12579/250‑1 (MQ=255)
                                                                                              cAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCgc                                                                                                                                                     >  1:208926/1‑250 (MQ=255)
                                                                                                     ggCGCGGCATCGAGGATCATCCGGGCTTCCCCGCCCAGTTTCGGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCTACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAAGGCGTCGTTCAGTGGCTTGTGCCGCTTCGACGCgg                                                                                                                                                <  1:339427/249‑1 (MQ=255)
                                                                                                           ggagggaggacaagcggggCTGCACCGCCTAGTTTAAGCGGGCCCGCGTGTTCGGTTGACAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGGCTGGCTTGGCAACCCGATGCTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCAGCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTgcgc                                                                                                                                        <  2:329035/235‑1 (MQ=255)
                                                                                                                gagtaccggcgggccttccccccccATTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGaa                                                                                                                                   <  2:278416/229‑1 (MQ=255)
                                                                                                                    gccAGCCGGGCTCCCCCGCCCGTGTTCAGCGGGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGat                                                                                                                              <  2:278168/249‑1 (MQ=255)
                                                                                                                     ccggCCGTGCTCCACCGCCGAGTTTGAGGTGGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGata                                                                                                                             <  2:139145/248‑1 (MQ=255)
                                                                                                                       aGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATAtt                                                                                                                           >  1:262498/1‑251 (MQ=255)
                                                                                                                          gggcctcccccGCCCAGTTCCAGGGGGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGgt                                                                                                                        <  2:79961/242‑1 (MQ=255)
                                                                                                                                              gcggcgcgCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGcc                                                                                                     <  2:278415/244‑1 (MQ=255)
                                                                                                                                                              gTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCg                                                                                       <  1:76818/249‑1 (MQ=255)
                                                                                                                                                                      ggCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGc                                                                             <  1:151565/251‑1 (MQ=255)
                                                                                                                                                                       gCAGGAAGAAAGTCCTATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCt                                                                            >  1:98228/1‑251 (MQ=255)
                                                                                                                                                                           gaagaaAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGCCCGTCACGCGGGGCTTTCAGTTCGTTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTCTTTCCCTTCGACCCGGGGTTGCGCCTGAATGGATATGGGGCGTGGCCCCTCTATAGCCGCCTTAGAATCCTGAACCTTAGCTTTCg                                                                        >  2:251364/1‑251 (MQ=255)
                                                                                                                                                                             agaaAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCCCCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTGGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGcc                                                                      <  1:318581/251‑1 (MQ=255)
                                                                                                                                                                                                gCTGAGATCGACCATATTCAGCACCCGCCCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCggg                                                   <  2:60234/251‑1 (MQ=255)
                                                                                                                                                                                                     gaTCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGTCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGGTTGCGCCTGAATGATATTTGTGCGTGCCGCTTTTATAGCCGCTTTAGAAGCCGATACCTGATCTTTCGCCGATTCCAGCGCAGCTCCGGCGctc                                              >  2:243638/1‑251 (MQ=255)
                                                                                                                                                                                                           ccATATTCACCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAAGGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATTATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAgcg                                        <  2:205353/251‑1 (MQ=255)
                                                                                                                                                                                                                 tCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAgcggcggc                                   <  1:161997/250‑1 (MQ=255)
                                                                                                                                                                                                                                 ggccgccAGCACTTCGCCTGGCTTGGCAACCCGATACTGCGCGCGTCCGTCACGCGGTGCTTTCAGTTCGCTGTCATCGATATCTTCGGCAATGCGCCTTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGGGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTGTCGCAGATTCCAGCGCAGTTCGGTCGGTCTCAGCGGCGGCGCGATCGTCATCCAGCt                  <  1:169916/250‑1 (MQ=255)
                                                                                                                                                                                                                                  gccgccAGCACTTCGCCTGGCTTGGCAAGCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGGTTTCGCCTGAATGATATTGGTGCGTGCCGCTTCCATAGCCCCTTTTGACCCCCATACCTGAGCTTTCGCCGAATTCAGCGCAGCTCCGGCGCCCTCCGCGGAGACGCGCTCCTCATcccactg                 >  2:241802/1‑246 (MQ=255)
                                                                                                                                                                                                                                     gccAGCACTTCGCCGGGCTGGGAAACCCGATACTGCAGGCGTCCGGCACGCGGGGCTCTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTGCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCACAGCGGCGGCGCGAGCGTCATCCAGCTGTTg              <  2:355968/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aCTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGGTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGGTCGTCATCCAGCTGTTgcgc           >  2:87599/1‑248 (MQ=255)
                                                                                                                                                                                                                                                 ccTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGc  >  1:315914/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
GGAACATCAGTTTCAGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTCGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGC  >  NC_000913/3629675‑3630164

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: