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→GCAtsrmethyl‑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

TCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAAT  >  NC_000913/4591729‑4592214
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
ccAGGCGGTCTTGTCTATAATGCCGTAAAGAAGTACAAAGAAAATTTCGCGTTTTTACAAACCATTCGCCAGCAGCAATCCACGCTTAATTGCAGCTGGTTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTa                                                                                                                                                                                                                                             <  2:459231/250‑1 (MQ=255)
 cAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTa                                                                                                                                                                                                                                             <  1:102701/250‑1 (MQ=255)
        tCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGt                                                                                                                                                                                                                                     >  1:16919/1‑251 (MQ=255)
           gTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGc                                                                                                                                                                                                                                  <  1:348565/251‑1 (MQ=255)
           gTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGGTTTACAAACCATTCGCCAGCAGCAATCCACGCGGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGtt                                                                                                                                                                                                                                    <  2:16919/249‑1 (MQ=255)
                ttAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGt                                                                                                                                                                                                                             >  1:128504/1‑251 (MQ=255)
                      ccTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCcg                                                                                                                                                                                                                       >  1:177420/1‑251 (MQ=255)
                              aaTGACAAAGAAAATTTCACTGGTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCagag                                                                                                                                                                                                               <  2:383126/251‑1 (MQ=255)
                               aTGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGc                                                                                                                                                                                                              <  2:216649/251‑1 (MQ=255)
                                   cAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGcc                                                                                                                                                                                                                                 >  1:381087/1‑228 (MQ=255)
                                   cAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGcc                                                                                                                                                                                                                                 <  2:381087/228‑1 (MQ=255)
                                   aaaaGTAAATTTCACTTGTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCAccg                                                                                                                                                                                                          <  2:97640/250‑1 (MQ=255)
                                       gAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACcgccgc                                                                                                                                                                                                      >  1:188371/1‑251 (MQ=255)
                                         aaaTTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCa                                                                                                                                                                                                             >  2:459899/1‑242 (MQ=255)
                                         aaaTTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCa                                                                                                                                                                                                             <  1:459899/242‑1 (MQ=255)
                                          aaTTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCcgcggg                                                                                                                                                                                                                            >  2:256015/1‑224 (MQ=255)
                                          aaTTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTg                                                                                                                                                                                                                            <  1:256015/226‑1 (MQ=255)
                                           aTTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAAACGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGGCAGAGACCCGCCGCAggg                                                                                                                                                                                                  >  2:104006/1‑249 (MQ=255)
                                             ttCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCGGAATGGCCGCGGGGTCGGGTTGTTGCAGACGCGTAACACCCTCACCCGCGCGGGTACCCGCTACATGAGGGATCAGAAAAATATTGGTAGCGGTTCAACCGTGGCTGAGCTGATGCAGAGTGCCAGTATTCCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGTGTGGCCGCGTGACCCGCGTAAGAGCACCGCCGCAGCgg                                                                                                                                                                                                 <  2:404226/250‑1 (MQ=255)
                                                aCTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCACGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAACAAACTGGGCAGATTACGACGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCagc                                                                                                                                                                                             >  1:47336/1‑250 (MQ=255)
                                                  tGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGACTGGCAGCTGGGTCGCGTTGTTGAAGACGCGGAACACCCTCAACCGCGCGGGTATCCGCTTCATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGAACCGCGTCAGAGCACCGCCGACGCGGCag                                                                                                                                                                                              >  1:307473/1‑248 (MQ=255)
                                                   gTTTTATAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCGGGGTCGCGTTGTTGCAGACGCGTACCCCCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTAAACCTTTTCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCCCCGCCGCAGCGGCagag                                                                                                                                                                                            <  2:128504/249‑1 (MQ=255)
                                                      ttACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCa                                                                                                                                                                                        >  1:11978/1‑250 (MQ=255)
                                                        aCAACCCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGAAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGAAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAAc                                                                                                                                                                                     <  1:52429/251‑1 (MQ=255)
                                                            aCCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGACACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGAtacaacgca                                                                                                                                                                                  >  2:24350/1‑242 (MQ=255)
                                                             ccATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAAt                                                                                                                                                                                >  1:485162/1‑251 (MQ=255)
                                                                     gcgcaggaATACACGCTGATTGGCAGCTGGGTCGCGTTGGTGCAGACGCGTAACACCCTCAACCGCGAGGGAATCCGCTACATGTTGGATCAGAATAACATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGATGAAACAGGCGGAAACACAACTGGGCGGATTAAGAAGCCTTGCCGCGTGAACCGCGTCAGAGAACCGCAACAGGGGCCGAGATCAAACGCAATTACGat                                                                                                                                                                          >  2:40486/3‑250 (MQ=255)
                                                                            aaTCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACGGGGCAGATTACGAAGCGTTGCCGCGTGACCTTCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATcac                                                                                                                                                                 <  1:215499/251‑1 (MQ=255)
                                                                             aTCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTCCGATATTTATcac                                                                                                                                                                 >  2:142926/1‑250 (MQ=255)
                                                                              cccgcgcttaagggCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCACCCGCGGGGGTTGCCGCTACATGATGGATCAGAAAAATATGGGGAGCGGTACAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTAAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACaa                                                                                                                                                               <  2:177420/239‑1 (MQ=255)
                                                                               ccACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCACACGTAATTACGATATTTATCACAAt                                                                                                                                                              >  1:320014/1‑251 (MQ=255)
                                                                                     ggAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTATCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATAGTTATCACAATGCGCTg                                                                                                                                                        <  1:376899/250‑1 (MQ=255)
                                                                                          ggCAGTTGGGTCGCGTTTTTGCAGACGCGTAACCCCCTCACCCGCCCGGGTACCCGCTACATGATGGATCAGAATATAATTGGTAGCGGTTCAACCGTGGCGGATCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAACCGGGGCAGATTACGAAGCTTTGCCGCGTGACCCGCTTCATAGCTCCGCCGCAGCGGCAGAGATCAAACGTATGTACGATATTTATCACAATGCGCTGGCGGa                                                                                                                                                   <  1:159754/251‑1 (MQ=255)
                                                                                           gCAGCTGTGTCGCGTTGTTGCAGACGCTTAACACCCTAAACCGCGCGGTTATCCGCTACAGGATGGACCATAATAATAGTGGTAGCGGTTCAACCGTTGCTGAGCTGATGAAGAGTGCCAGTATTTCGCTGAAACAGGCGAAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGCCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAg                                                                                                                                                  <  1:274249/251‑1 (MQ=255)
                                                                                                   gTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACCTGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCCCCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATcc                                                                                                                                           >  1:370718/1‑250 (MQ=255)
                                                                                                                                    cgcgcgGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCtt                                                                                                         >  1:166110/1‑251 (MQ=255)
                                                                                                                                      ggcgggttccccctccatgaggaaTCAAAATATTATTGGTGGCGATTCAACCGTTGCGGATCTGATGCGGATTCCCAGTATTCCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCTTTGCCTCGTGTCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACGATTAGGTGCAGGCAAAATCAACGAGTTCTTTg                                                                                                       <  2:249375/228‑1 (MQ=255)
                                                                                                                                       gcgGGTTTCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGGGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGa                                                                                                      <  1:493244/251‑1 (MQ=255)
                                                                                                                                                                        atTGGTAGCGGTTCACCCGTTGCGGAGCTGTGGAAGGGGGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACGGGGCAGATTACGAAGCGTTGCGTCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCga                                                                     <  2:135781/251‑1 (MQ=255)
                                                                                                                                                                              aGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGc                                                                >  1:136955/1‑250 (MQ=255)
                                                                                                                                                                                         cGTTGCTGGGCTGATGCAGAGTGCCAGTATTTCGCTGAAACGGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTTATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCtt                                                    <  1:247056/251‑1 (MQ=255)
                                                                                                                                                                                                 ggctgaTGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAg                                            <  1:411632/250‑1 (MQ=255)
                                                                                                                                                                                                          agagTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGCCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGAt                                   <  1:131974/251‑1 (MQ=255)
                                                                                                                                                                                                               tGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGc                               >  1:36683/1‑251 (MQ=255)
                                                                                                                                                                                                                        ttCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGat                       >  1:176451/1‑250 (MQ=255)
                                                                                                                                                                                                                           gCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCg                   <  1:134409/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ggCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTcta                           <  1:309714/234‑3 (MQ=255)
                                                                                                                                                                                                                                    ggCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGGGGCTTACCTGGGGCAAAACGATCGGCTctc                           >  2:309714/1‑232 (MQ=255)
                                                                                                                                                                                                                                       ggAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGAt        >  1:450426/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   gggCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAAt  >  1:400388/1‑245 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
TCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAAT  >  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: