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

GCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAG  >  NC_000913/4591685‑4592172
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
gcctACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAAGGACAAAGAAAATTTCCCTGTTGTACAAACCATTCGCGAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACAACCTCAACCGCGCGGGTATCCGCTACATGATGGAGCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGt                                                                                                                                                                                                                                                <  1:439635/248‑1 (MQ=255)
        ctgGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGGAGCGGGTCAACCGTTGCTGAGCTGATGCAGAGGGCCAGTAt                                                                                                                                                                                                                                        >  2:376899/1‑251 (MQ=255)
          ggttttggccggttttggccTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATTCCCTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTGCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTAttt                                                                                                                                                                                                                                      <  2:242501/251‑1 (MQ=255)
          ggTTTTGGCCGTTTTTGGTCGTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTAtt                                                                                                                                                                                                                                       <  1:364333/250‑1 (MQ=255)
                      ttttGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAg                                                                                                                                                                                                                           <  1:166013/250‑1 (MQ=255)
                                      cTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGc                                                                                                                                                                                                         >  1:75260/1‑251 (MQ=255)
                                            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                                                                                                                                                                                       <  2:381087/228‑1 (MQ=255)
                                                                               cAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGcc                                                                                                                                                                                       >  1:381087/1‑228 (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)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
GCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAG  >  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: