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,9300GC92.2% 138.7 / 3.8 53E92Q (GAG→CAG) tsrmethyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base C (31/16);  total (34/19)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

CTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAA  >  NC_000913/4591686‑4592175
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
cTTACTGCTGGTTTAGGCCGTTTTTGGCCTTTTACAACTGCAATAAGTCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCCCGCTGAATGGCAGCTGGGTCGCGATGTTGCAAACGAGTAACACACTCAACCGCGCGGGTATCCGCTACGTGATGGATCAGAATAATCTTAGTAGCGGTTCAACCGTTGCTGAGCTGTTGAAGAGTg                                                                                                                                                                                                                                                 >  2:478844/1‑251 (MQ=255)
           ttttGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCg                                                                                                                                                                                                                                      >  1:180391/1‑251 (MQ=255)
            tttGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCCGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATACACTTGGTAGCGGTTCAACCGCTGCTGAGCTGATGCAGAGTGCAAGTAttttcc                                                                                                                                                                                                                                     >  2:667071/1‑248 (MQ=255)
                  cGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGCGTGCCAGTATTTCGCTGaacc                                                                                                                                                                                                                               >  2:241446/1‑249 (MQ=255)
                  cGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAAc                                                                                                                                                                                                                               <  1:65990/251‑1 (MQ=255)
                       ttGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATACGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCg                                                                                                                                                                                                                          >  1:50273/1‑251 (MQ=255)
                               ttACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAACTTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCACCCGCGCGGTGATCCGCTACATGATGGATCCGAATAATATTGGTCGCGGTTCAACCGTTGCTGAGCTGATGCAGCGTGCCAGTATTTCGCTGCACCAGGCGGcaaaaa                                                                                                                                                                                                                   >  2:305730/1‑250 (MQ=255)
                                    aCTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGACAGTATTTCGCTGAAACAGGCGGAAAAAAAATggg                                                                                                                                                                                                             >  2:563079/1‑251 (MQ=255)
                                       gACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGca                                                                                                                                                                                                           <  1:241446/250‑2 (MQ=255)
                                             aGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACg                                                                                                                                                                                                    <  1:305730/251‑1 (MQ=255)
                                                cgggctgatCTTTAATACATTAAAGAATGACAAAGAAAATTTCCCTGTTTTACAAAACCTTCGCCAGCAACAATCCACTCTGAATAGCAGCTGGGTCGCGTTGATGCAGACGCGTCACACCCTCAACAGCGAGGGTATACGCTACATGCAACAACAGAATAATCTTGGTAGAGGTTCAACCGTTGCTGAGCTGATGCCGAGTGCCAGTATTTCGCTGAAACAGGCGGAACACAACTggtagaagtcaga                                                                                                                                                                                                   >  2:625662/9‑238 (MQ=255)
                                                      gTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTg                                                                                                                                                                                            <  1:226505/250‑1 (MQ=255)
                                                                 ccTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCcg                                                                                                                                                                                >  1:86895/1‑251 (MQ=255)
                                                                   ttAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAAATGGGCAGATTACGAAGCGTTGCCGCGTGAccctcg                                                                                                                                                                              >  1:9895/1‑248 (MQ=255)
                                                                       aGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGaa                                                                                                                                                                                                                                 >  2:475879/1‑196 (MQ=255)
                                                                       aGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGaa                                                                                                                                                                                                                                 <  1:475879/196‑1 (MQ=255)
                                                                                aaGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCAccgcc                                                                                                                                                                 >  1:179052/1‑251 (MQ=255)
                                                                                      aTTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCg                                                                                                                                                           >  1:64087/1‑251 (MQ=255)
                                                                                                  tACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGATGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGCGCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCaa                                                                                                                                                >  2:168594/1‑250 (MQ=255)
                                                                                                    cAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGAGGGAGAGTGCAAGTATTTCGCTGAAACAGGCGGAAAAAAACGGGGCGGATTACGAAGCGTTGCCGCGTGACCGGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACg                                                                                                                                             <  2:61408/251‑1 (MQ=255)
                                                                                                          aTTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAAtt                                                                                                                                        >  1:640318/1‑250 (MQ=255)
                                                                                                             cGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACg                                                                                                                                     >  1:430614/1‑250 (MQ=255)
                                                                                                                      caATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATc                                                                                                                            >  1:598610/1‑250 (MQ=255)
                                                                                                                            ccGCTGAATGGCAGCTGGGTCGCGTTGTGGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTAAATGATGGATAAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATgc                                                                                                                     <  1:484928/250‑1 (MQ=255)
                                                                                                                               cTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGAAGCGGCAGAGATCAAACGTAATTACGATTTTTATCACAATGCGCt                                                                                                                  >  2:292396/1‑251 (MQ=255)
                                                                                                                                 gAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTg                                                                                                                 <  1:339553/250‑1 (MQ=255)
                                                                                                                                      gccgCTGGGTCGCGTGTGTGCAGACGCTTAACACCCTCAACCGCGCGGGTAGCCCCTACAGGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGAAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACGGGGAAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAg                                                                                                           <  2:572196/248‑1 (MQ=255)
                                                                                                                                          cTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGCGCTGATGCAGAGTGCCAGTATTTCGATGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGTCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGCTCAACCGTAATTTCGATATTATTCACAATGAGCCGGCAGAGCTgg                                                                                                       >  2:158947/1‑250 (MQ=255)
                                                                                                                                            gggTGGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGGTGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTGCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGGCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATc                                                                                                     <  1:612537/251‑1 (MQ=255)
                                                                                                                                              gTCGCTTTGTTGCAGACGCTTAACCCCCTCACCCGCGCGGGTATCCCCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGATGATGCAGAGTGCCAGTATTTCGCTGAACCAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGCCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCa                                                                                                   <  1:558259/251‑1 (MQ=255)
                                                                                                                                                   gttgttGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACtat                                                                                              >  1:12759/1‑249 (MQ=255)
                                                                                                                                                      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)
                                                                                                                                                                                                                                        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)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
CTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAA  >  NC_000913/4591686‑4592175

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: