Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,591,974 G→A A106A (GCG→GCA)  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,9740GA98.0% 150.7 / NA 50A106A (GCG→GCA) tsrmethyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (23/26);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAAT  >  NC_000913/4591729‑4592214
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
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                                                                                                                                                                                                                                 >  1:381087/1‑228 (MQ=255)
                                   cAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGcc                                                                                                                                                                                                                                 <  2:381087/228‑1 (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)
                                                                                                                                                                                                          a‑gagTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGCCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGAt                                   <  1:131974/251‑1 (MQ=255)
                                                                                                                                                                                                               tGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGc                               >  1:36683/1‑251 (MQ=255)
                                                                                                                                                                                                                        ttCGCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGat                       >  1:176451/1‑250 (MQ=255)
                                                                                                                                                                                                                           gCTGAAACAGGCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCg                   <  1:134409/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ggCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTcta                           <  1:309714/234‑3 (MQ=255)
                                                                                                                                                                                                                                    ggCGGAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGGGGCTTACCTGGGGCAAAACGATCGGCTctc                           >  2:309714/1‑232 (MQ=255)
                                                                                                                                                                                                                                       ggAAA‑AAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGAt        >  1:450426/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   gggCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAAt  >  1:400388/1‑245 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
TCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAA‑AAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAAT  >  NC_000913/4591729‑4592214

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