Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,591,930 G→C E92Q (GAG→CAG)  tsr → methyl‑accepting chemotaxis protein Tsr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,591,9300GC97.9% 143.2 / NA 48E92Q (GAG→CAG) tsrmethyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base C (21/25);  minor base A (0/1);  total (21/27)
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

GCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAG  >  NC_000913/4591685‑4592172
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
gcctACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAAGGACAAAGAAAATTTCCCTGTTGTACAAACCATTCGCGAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACAACCTCAACCGCGCGGGTATCCGCTACATGATGGAGCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGt                                                                                                                                                                                                                                                <  1:439635/248‑1 (MQ=255)
        ctgGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGGAGCGGGTCAACCGTTGCTGAGCTGATGCA‑GAGGGCCAGTAt                                                                                                                                                                                                                                        >  2:376899/1‑251 (MQ=255)
          ggttttggccggttttggccTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATTCCCTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTGCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTAttt                                                                                                                                                                                                                                      <  2:242501/251‑1 (MQ=255)
          ggTTTTGGCCGTTTTTGGTCGTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTAtt                                                                                                                                                                                                                                       <  1:364333/250‑1 (MQ=255)
                      ttttGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAg                                                                                                                                                                                                                           <  1:166013/250‑1 (MQ=255)
                                      cTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGc                                                                                                                                                                                                         >  1:75260/1‑251 (MQ=255)
                                            ccAGGCGGTCTTGTCTATAATGCCGTAAAGAAGTACAAAGAAAATTTCGCGTTTTTACAAACCATTCGCCAGCAGCAATCCACGCTTAATTGCAGCTGGTTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTa                                                                                                                                                                                                   <  2:459231/250‑1 (MQ=255)
                                             cAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTa                                                                                                                                                                                                   <  1:102701/250‑1 (MQ=255)
                                                    tCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGt                                                                                                                                                                                           >  1:16919/1‑251 (MQ=255)
                                                       gTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGc                                                                                                                                                                                        <  1:348565/251‑1 (MQ=255)
                                                       gTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGGTTTACAAACCATTCGCCAGCAGCAATCCACGCGGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGtt                                                                                                                                                                                          <  2:16919/249‑1 (MQ=255)
                                                            ttAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGt                                                                                                                                                                                   >  1:128504/1‑251 (MQ=255)
                                                                  ccTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCcg                                                                                                                                                                             >  1:177420/1‑251 (MQ=255)
                                                                          aaTGACAAAGAAAATTTCACTGGTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCagag                                                                                                                                                                     <  2:383126/251‑1 (MQ=255)
                                                                           aTGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGc                                                                                                                                                                    <  2:216649/251‑1 (MQ=255)
                                                                               cAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGcc                                                                                                                                                                                       <  2:381087/228‑1 (MQ=255)
                                                                               cAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGcc                                                                                                                                                                                       >  1:381087/1‑228 (MQ=255)
                                                                               aaaaGTAAATTTCACTTGTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCAccg                                                                                                                                                                <  2:97640/250‑1 (MQ=255)
                                                                                   gAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACcgccgc                                                                                                                                                            >  1:188371/1‑251 (MQ=255)
                                                                                     aaaTTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCa                                                                                                                                                                   >  2:459899/1‑242 (MQ=255)
                                                                                     aaaTTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCa                                                                                                                                                                   <  1:459899/242‑1 (MQ=255)
                                                                                      aaTTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCcgcggg                                                                                                                                                                                  >  2:256015/1‑224 (MQ=255)
                                                                                      aaTTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTg                                                                                                                                                                                  <  1:256015/226‑1 (MQ=255)
                                                                                       aTTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAAACGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGGCAGAGACCCGCCGCAggg                                                                                                                                                        >  2:104006/1‑249 (MQ=255)
                                                                                         ttCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCGGAATGGCCGCGGGGTCGGGTTGTTGCAGACGCGTAACACCCTCACCCGCGCGGGTACCCGCTACATGAGGGATCAGAAAAATATTGGTAGCGGTTCAACCGTGGCTGAGCTGATGCA‑GAGTGCCAGTATTCCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGTGTGGCCGCGTGACCCGCGTAAGAGCACCGCCGCAGCgg                                                                                                                                                       <  2:404226/250‑1 (MQ=255)
                                                                                            aCTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCACGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑CAAACTGGGCAGATTACGACGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCagc                                                                                                                                                   >  1:47336/1‑250 (MQ=255)
                                                                                              tGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGACTGGCAGCTGGGTCGCGTTGTTGAAGACGCGGAACACCCTCAACCGCGCGGGTATCCGCTTCATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGAACCGCGTCAGAGCACCGCCGACGCGGCag                                                                                                                                                    >  1:307473/1‑248 (MQ=255)
                                                                                               gTTTTATAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCGGGGTCGCGTTGTTGCAGACGCGTACCCCCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTAAACCTTTTCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCCCCGCCGCAGCGGCagag                                                                                                                                                  <  2:128504/249‑1 (MQ=255)
                                                                                                  ttACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCa                                                                                                                                              >  1:11978/1‑250 (MQ=255)
                                                                                                    aCAACCCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGAAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGAAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAAc                                                                                                                                           <  1:52429/251‑1 (MQ=255)
                                                                                                        aCCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGACACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGAtacaacgca                                                                                                                                        >  2:24350/1‑242 (MQ=255)
                                                                                                         ccATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAAt                                                                                                                                      >  1:485162/1‑251 (MQ=255)
                                                                                                                 gcgcaggaATACACGCTGATTGGCAGCTGGGTCGCGTTGGTGCAGACGCGTAACACCCTCAACCGCGAGGGAATCCGCTACATGTTGGATCAGAATAACATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGATGAAACAGGCGGAAACACAACTGGGCGGATTAAGAAGCCTTGCCGCGTGAACCGCGTCAGAGAACCGCAACAGGGGCCGAGATCAAACGCAATTACGat                                                                                                                                >  2:40486/3‑250 (MQ=255)
                                                                                                                        aaTCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACGGGGCAGATTACGAAGCGTTGCCGCGTGACCTTCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATcac                                                                                                                       <  1:215499/251‑1 (MQ=255)
                                                                                                                         aTCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTCCGATATTTATcac                                                                                                                       >  2:142926/1‑250 (MQ=255)
                                                                                                                          cccgcgcttaagggCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCACCCGCGGGGGTTGCCGCTACATGATGGATCAGAAAAATATGGGGAGCGGTACAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTAAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACaa                                                                                                                     <  2:177420/239‑1 (MQ=255)
                                                                                                                           ccACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCACACGTAATTACGATATTTATCACAAt                                                                                                                    >  1:320014/1‑251 (MQ=255)
                                                                                                                                 ggAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTATCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATAGTTATCACAATGCGCTg                                                                                                              <  1:376899/250‑1 (MQ=255)
                                                                                                                                      ggCAGTTGGGTCGCGTTTTTGCAGACGCGTAACCCCCTCACCCGCCCGGGTACCCGCTACATGATGGATCAGAATATAATTGGTAGCGGTTCAACCGTGGCGGATCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAACCGGGGCAGATTACGAAGCTTTGCCGCGTGACCCGCTTCATAGCTCCGCCGCAGCGGCAGAGATCAAACGTATGTACGATATTTATCACAATGCGCTGGCGGa                                                                                                         <  1:159754/251‑1 (MQ=255)
                                                                                                                                       gCAGCTGTGTCGCGTTGTTGCAGACGCTTAACACCCTAAACCGCGCGGTTATCCGCTACAGGATGGACCATAATAATAGTGGTAGCGGTTCAACCGTTGCTGAGCTGATGAA‑GAGTGCCAGTATTTCGCTGAAACAGGCGAAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGCCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAg                                                                                                        <  1:274249/251‑1 (MQ=255)
                                                                                                                                               gTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACCTGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCCCCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATcc                                                                                                 >  1:370718/1‑250 (MQ=255)
                                                                                                                                                                                cgcgcgGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCtt                                                               >  1:166110/1‑251 (MQ=255)
                                                                                                                                                                                  ggcgggttccccctccatgaggaaTCAAAATATTATTGGTGGCGATTCAACCGTTGCGGATCTGATGCG‑GATTCCCAGTATTCCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCTTTGCCTCGTGTCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACGATTAGGTGCAGGCAAAATCAACGAGTTCTTTg                                                             <  2:249375/228‑1 (MQ=255)
                                                                                                                                                                                   gcgGGTTTCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGGGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGa                                                            <  1:493244/251‑1 (MQ=255)
                                                                                                                                                                                                                    atTGGTAGCGGTTCACCCGTTGCGGAGCTG‑TGGAAGGGGGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACGGGGCAGATTACGAAGCGTTGCGTCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCga                           <  2:135781/251‑1 (MQ=255)
                                                                                                                                                                                                                          aGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGc                      >  1:136955/1‑250 (MQ=255)
                                                                                                                                                                                                                                     cGTTGCTGGGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACGGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTTATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCtt          <  1:247056/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ggctgaTGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAg  <  1:411632/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
GCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAG  >  NC_000913/4591685‑4592172

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: ‑