Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,592,082 G→A L142L (CTG→CTA tsr → methyl‑accepting chemotaxis protein Tsr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,592,0820GA93.7% 170.7 / 2.6 64L142L (CTG→CTAtsrmethyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (32/27);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

GTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCG  >  NC_000913/4591836‑4592327
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
gttgCAGACGCTTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATGGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAGTTCGCTGGCGGAGCTGATCCAACTATTAg                                                                                                                                                                                                                                                   <  1:530119/251‑1 (MQ=255)
            tAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGCAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCCAACGTAATTACGACATTTATCACAATGCGCTGGCGGAGCTGATCCAAATATTAGGTGCAGGCaaa                                                                                                                                                                                                                                        >  2:599287/1‑250 (MQ=255)
                          gcgcgGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCttt                                                                                                                                                                                                                         <  1:292396/251‑1 (MQ=255)
                               ggTACCCGCTAAATGATGGATCAAAATAATTGTGGGAGCGGTTCAACCTTTGCTGTGCTGATGCAGAGTGCAATTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCTTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCa                                                                                                                                                                                                                    <  2:670027/251‑1 (MQ=255)
                                        tACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGAccc                                                                                                                                                                                                            >  1:265409/1‑250 (MQ=255)
                                          cATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAgg                                                                                                                                                                                                         >  1:359873/1‑251 (MQ=255)
                                              atgGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGCAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGAGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGat                                                                                                                                                                                                      >  1:101278/1‑250 (MQ=255)
                                                  aacAGAATAATATTGGAAGAGGTTCAACCGTTGCTGCGCTGATGCAGAATGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGCGATCAAACGAAATTACGATATTTATCACACTGCGCTGGCGGGGCTGATCCAACTATTAGGGGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGatata                                                                                                                                                                                                   >  2:70946/3‑248 (MQ=255)
                                                      gAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACg                                                                                                                                                                                             >  1:565967/1‑251 (MQ=255)
                                                             atTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCgg                                                                                                                                                                                      >  1:10806/1‑250 (MQ=255)
                                                                ggtagggtttCAACCGTGTCGGAGCTGAGGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACGGGGCAGATTACGAAGCGTGGCCGCGTGCCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCgaga                                                                                                                                                                                    <  2:12759/242‑1 (MQ=255)
                                                                ggTAGCGGTTCAACCGTGGCTGAGCTGAGGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTCCGAAGCGTTGCCGCGTGACCCGCTTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGaa                                                                                                                                                                                   <  1:667071/251‑1 (MQ=255)
                                                                  tAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATCTCAGGACGGTTTCGAGAAGc                                                                                                                                                                                 >  1:242068/1‑251 (MQ=255)
                                                                      ggTTCAACCGTTGCTGAGCTGATGCAGAGAGCCAGTCTTTCGCTGAAACCGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGCGCACCGCCACAGCGGCAGAGCTCAAACGTAATTACGGTATTTATCACAATGCGCTGGCGGAGCTGATCCCACTATTAGGTGCAGGCAAACTCAACGAGTTCTTTGATCAACCGACCCCGGGATATCTGGACGGGTTCGAGAAgatgtc                                                                                                                                                                             >  2:130292/1‑246 (MQ=255)
                                                                       gttaaaCCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGTGGAAACAGGCGGAAAAAAACGGGTAAGATTACGAAGCGTTGCCTGGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTAt                                                                                                                                                                            <  1:222358/247‑1 (MQ=255)
                                                                          cAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCACCGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTAtgtg                                                                                                                                                                         >  1:240741/1‑251 (MQ=255)
                                                                            aCCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGc                                                                                                                                                                       >  1:474489/1‑251 (MQ=255)
                                                                                ttgcggagcggatgcagatggCCAGTATTTCCGTGAAACAGGGGGAAAAAAGCTGGGCGGATTACGAAGCGTTGCCTCGTGACCCGCGTAAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCATCTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCtt                                                                                                                                                                     <  1:478844/229‑1 (MQ=255)
                                                                                 tGCTGAGCTGCTGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGGTATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCCGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCCGGGATATCAGGACGTATTCGCGATGCAGTATGTGGCTTAc                                                                                                                                                                   >  2:327132/1‑250 (MQ=255)
                                                                                     gagctaaTGCAGAGTGCAAGTATTGCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCTTGACCCGCGTAAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATgg                                                                                                                                                               <  1:327132/250‑1 (MQ=255)
                                                                                          gaTGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGtt                                                                                                                                                                                                                                                     <  1:621232/159‑1 (MQ=255)
                                                                                          gaTGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGtt                                                                                                                                                                                                                                                     >  2:621232/1‑159 (MQ=255)
                                                                                          gaTGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCaa                                                                                                                                                          >  1:536509/1‑250 (MQ=255)
                                                                                            tgcagagTGCCAGTATTTAGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGGTCTTTGATCAGCCGACCCAGGGAGATCAGGACGGCTTCGAGAAGCAGGAGGAGGCTTACATGGAGCAAAAc                                                                                                                                                       >  2:74241/4‑251 (MQ=255)
                                                                                               agagTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGAt                                                                                                                                                    >  1:638101/1‑251 (MQ=255)
                                                                                                gagTCCCAGTATTTCGCTGAAACAGGCGGAAAAAAACGGGGCAGATTACGAAGCGTTGCCGCGGGCCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCGTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATc                                                                                                                                                   <  2:265409/251‑1 (MQ=255)
                                                                                                 agTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCg                                                                                                                                                  <  1:168594/251‑1 (MQ=255)
                                                                                                        gTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTctat                                                                                                                                           >  1:563995/1‑248 (MQ=255)
                                                                                                          aTTTCGCTGAAACAGGCGGAAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGGGCAGGCAAAATCAACGAGTTCTTTGATCAGCAGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGGGCAAAACGATCGGCTCTATg                                                                                                                                          >  2:186427/1‑251 (MQ=255)
                                                                                                           ttccGCTGAAACAGGCGGAAAAAAACGGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGc                                                                                                                                                                                 <  1:96898/207‑1 (MQ=255)
                                                                                                            ttcttataaaCAGGCGGAAAAAATCTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCGCCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCTGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGGGGCTTACATGGAGCAAAACGATCGGCTCTATGat                                                                                                                                        <  2:617465/243‑1 (MQ=255)
                                                                                                                   caaCAGGCGGAAAAAAACTGGGCAGATTACGAAACGATGCCGCGTGACCCGCGTCAGAGAACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGTGCTGATCCAACTCATAGGTGCAGGAAAAATCAGCGAGTTCTTTGATCAGCCCACACAGGGATAGCAGGACGGTTTCGAGAAGCAGAAGGAGGCTTACATGGAGCAAAACGaaggggtatatgttatacccgt                                                                                                                                >  2:410139/2‑230 (MQ=255)
                                                                                                                       aGGCGGAAAAAAACTGGGCGGATTACGAATCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGc                                                                                                                            >  1:346894/1‑251 (MQ=255)
                                                                                                                        ggCGGAAAAAACCTGGGCAGATTACGAAGCGTTGCCGCGTAACCCGCGTCAGAGCACCGCCGCAGCGGCAGATATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCg                                                                                                                           <  2:175786/251‑1 (MQ=255)
                                                                                                                        ggCGGAAAAAAACGGGGCAGATTAAGAACGTTTGCCGTGCGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTAGCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCg                                                                                                                           <  1:21888/251‑1 (MQ=255)
                                                                                                                                 aaaCTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATg                                                                                                                  >  2:597824/1‑251 (MQ=255)
                                                                                                                                        ggcagATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTc                                                                                                                                                                                        <  2:99734/171‑1 (MQ=255)
                                                                                                                                        ggcagATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTc                                                                                                                                                                                        >  1:99734/5‑175 (MQ=255)
                                                                                                                                                  gAAGCGTGGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATAGCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGc                                                                                                  <  1:604977/251‑1 (MQ=255)
                                                                                                                                                         tGCCGCGGGCCCCGCGTCAGGGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGAtgtg                                                                                            <  2:565967/250‑1 (MQ=255)
                                                                                                                                                                gTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGGGGATTCtg                                                                                     >  2:200548/1‑251 (MQ=255)
                                                                                                                                                                 tGACCCGAGTCAGAGCACCGCCGCAGCGGCAGAGATCAACCGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACACTGCCTCCTACAGCCAGGCGCGGTGGATTCtgg                                                                                    >  2:627561/1‑251 (MQ=255)
                                                                                                                                                                   ccccgtgtcagggacccgccgCAGCGGCAGAGATCAAACGTAATGACGATATTTATCACAATGCGCTGGCGGAGATGATCCAACTATTAGGTGCAGGCAAAAACAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCGATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCtggtg                                                                                  <  2:623581/236‑1 (MQ=255)
                                                                                                                                                                   ccccGCGTCAGAGTACCGACGCAGCGGCAGAGATCAAACGTAATTACGGTATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTGTGATCAGCCGACTCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATACCAATGCCTCCTACAGCCAGGCGATGTGGATTCtggtg                                                                                  <  2:474489/250‑1 (MQ=255)
                                                                                                                                                                          tCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCCCAATGCGCTGGCGGAGCTGATCCAACGATTAGGGGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCAAGGGATAGCAGGCCGGTTTCGAGAAGCAGTATGTGGCGTACAGGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCAATAACATTGCCACCTACAGCCAGGCGGTGGGGATTCTGGGGGGCggg                                                                            >  2:214044/1‑248 (MQ=255)
                                                                                                                                                                           cAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGCCCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGtgat                                                                          >  1:623099/1‑251 (MQ=255)
                                                                                                                                                                              agCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTgatgat                                                                       <  1:158947/251‑1 (MQ=255)
                                                                                                                                                                                ccccgccgCAGCGGCAGAGATCAAACGTAATTACGAGATTTATCAAAATGCGCGGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGAGTCTGGTGGGCGTGATGAtcg                                                                     <  2:113084/249‑1 (MQ=255)
                                                                                                                                                                                     ccgcAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTa                                                                <  1:597824/251‑1 (MQ=255)
                                                                                                                                                                                      cgcAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCCACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGCCCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCCAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTCCAGCCAGGc                                                                                                  >  1:121033/1‑216 (MQ=255)
                                                                                                                                                                                      cgcAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCTATGTGGATTCTGGTGGGCGTGATGATCGTCGTAc                                                               <  1:214044/251‑1 (MQ=255)
                                                                                                                                                                                              cAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTc                                                       >  1:20592/1‑251 (MQ=255)
                                                                                                                                                                                                     cAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTc                                                 <  1:375028/250‑1 (MQ=255)
                                                                                                                                                                                                     cAAAAGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATAAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTCCTGGCGGTCCTCTTcc                                                >  1:470428/1‑250 (MQ=255)
                                                                                                                                                                                                            aaTTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCttttgccgctg                                         >  1:670529/1‑242 (MQ=255)
                                                                                                                                                                                                                 cGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCAACTTCGCCGTCTGGTTCg                                    >  1:164928/1‑251 (MQ=255)
                                                                                                                                                                                                                           cacaATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGGATTTAAGc                          >  1:183639/1‑251 (MQ=255)
                                                                                                                                                                                                                             caATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATTATAGCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCt                        <  1:625662/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTTGGGATTAAAGCCTccactgg                 >  1:637739/1‑245 (MQ=255)
                                                                                                                                                                                                                                        cGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAgcg             <  1:563079/251‑1 (MQ=255)
                                                                                                                                                                                                                                                aTCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAACACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCCGGCGATGAGGATTCTGGTGGGCGTGATGCTCGTCGTACTGGCGGTGATCTTATCCGTCGGGTTCTGTTTTAAAGCCTCCCGGGCAGCGcccgtgc      >  1:608472/1‑245 (MQ=255)
                                                                                                                                                                                                                                                  ccAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATc   <  1:200548/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccAACTATTAGGGGCAGGTCAAATCAACGAGGTGTTGGAGCAGCCGACCCAGGGATATAAGGACGGTTTCGAGAAGCAGGATGTGGCTTACATGGAGCAAAACGATCGGCTCTAGGATATCGCCGTCAGCGCTAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGGGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAAt    <  2:179052/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   cAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTGCCGCCTGGTTCGGTTTTAAAGCCTCGCTGGGAGCGCCAATGaaacc  >  1:241371/1‑248 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
GTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCG  >  NC_000913/4591836‑4592327

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: