Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,857 A→T 5.9% intergenic (+24/+22) qseC → / ← ygiZ quorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,171,8570AT 5.9% 531.3 / 17.7 152intergenic (+24/+22)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (69/74);  new base T (6/3);  total (75/77)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.00e-01

Only 100 of 165 total aligned reads displayed.

CACAGGGCAGCGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGAAGGTATACTGCAAATATTTCCTCCCGACTCAAGCATCCAGTCATACTTATTGAA  >  NC_000913/3171589‑3172127
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                
cccagggcagggTGGTATACGTCACGCGTAATGCTGAGAATTTCAGCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑Ga                                                                                                                                                                                                                                                                           <  4:19401/273‑1 (MQ=255)
  cggggCAGGGTGGGCGCCGTAACGCTGNNTGCTGATGATTTCATCGTGGGGGATAACGGNCCCGGGGGGAAACCAGAGGCCCTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGNGGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑Ga                                                                                                                                                                                                                                                                           <  3:43234/271‑1 (MQ=255)
   aGGGCAGCGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTTCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCCATGTCTATGAGAGCCGTTaaa                                                                                                                                                                                                                                                                                >  2:252999/1‑268 (MQ=255)
          cGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTTGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTTAGCTGGTAACTTCTCTGTCTATTGGAGCCGTTAAAAC‑GACTCTCAAAGa                                                                                                                                                                                                                                                                 >  2:282236/1‑275 (MQ=255)
                               tGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaa                                                                                                                                                                                                                                                                                >  4:219231/1‑240 (MQ=255)
                                                               ccGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCCCCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGAAGGGCATGACAGTTTAATTTGGGAATGCGGAAAAAGGTGGATTTGAGGCGAAGGGAAGCTGGTAACTTCTCGGTCTATGAGAGCCTTTAAAAA‑GACCCCCATTGGTTTTGCCAATGACCAAAAATAACACGGCACCAAAAAAGTCAACACACCACaac                                                                                                                                                                                                            >  4:19621/1‑274 (MQ=255)
                                                                                    tggccagAATTGGCGCACGCTTCTATCGCCCACCCGGACAAATCGCTACCGGCATCGGGCGTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAAAAAGGTGGATTTGAGGCGAAGGTAATCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGaa                                                                                                                                                                                      <  2:46002/269‑1 (MQ=255)
                                                                                         cgAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑Ga                                                                                                                                                                                                                                                                           <  4:272046/186‑1 (MQ=255)
                                                                                          gAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GAct                                                                                                                                                                                                                                                                         >  2:354930/1‑187 (MQ=255)
                                                                                              tGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCa                                                                                                                                                                                                                                                                     <  4:354930/187‑1 (MQ=255)
                                                                                                                        gACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTa                                                                                                                                                                                                                                                                                  >  7:686900/1‑149 (MQ=255)
                                                                                                                          aagcccGCTGCCGGTAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaa                                                                                                                                                                                                                                                                                >  5:504284/5‑149 (MQ=255)
                                                                                                                           aaaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                                                                                                                                                               >  7:237736/1‑149 (MQ=255)
                                                                                                                           aaaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                                                                                                                                                               >  8:1124869/1‑149 (MQ=255)
                                                                                                                           aaaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGACCAGCGAATCGCCAAATTGCATGGAATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                                                                                                                                                               <  5:1233577/149‑1 (MQ=255)
                                                                                                                               cGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                                                                                                                                                               <  8:1146896/145‑1 (MQ=255)
                                                                                                                               cGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                                                                                                                                                               >  7:1146896/1‑145 (MQ=255)
                                                                                                                                gCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GAct                                                                                                                                                                                                                                                                         >  8:405145/1‑149 (MQ=255)
                                                                                                                                  tACCGGCAGCGGGCTTTGGGTATTGATGGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCtt                                                                                                                                        <  4:282236/275‑1 (MQ=255)
                                                                                                                                    ccGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCa                                                                                                                                                                                                                                                                     <  6:498170/149‑1 (MQ=255)
                                                                                                                                        cggCGGGCTTGGGCTATCGATTGTGCAGCGAATCGCCAAACTGCATGGCAGTAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagata                                                                                                                                  <  2:19621/273‑1 (MQ=255)
                                                                                                                                         aGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAAC‑GGCTCTCATAGac                                                                                                                                                                                                                                                                <  6:357834/149‑2 (MQ=255)
                                                                                                                                             ggCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTa                                                                                                                                                                                                                                                            <  8:1357197/149‑1 (MQ=255)
                                                                                                                                                 tGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTaa                                                                                                                                                                                                                                                        <  7:1438039/149‑1 (MQ=255)
                                                                                                                                                  gggCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAAt                                                                                                                                                                                                                                                       <  8:121235/149‑1 (MQ=255)
                                                                                                                                                   ggCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATa                                                                                                                                                                                                                                                      <  7:1039017/149‑1 (MQ=255)
                                                                                                                                                      tATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCa                                                                                                                                                                                                                                                   >  6:1189521/1‑149 (MQ=255)
                                                                                                                                                       aTCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCaa                                                                                                                                                                                                                                                  >  6:1467250/1‑149 (MQ=255)
                                                                                                                                                       aTCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCaa                                                                                                                                                                                                                                                  >  6:1355581/1‑149 (MQ=255)
                                                                                                                                                          gATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAat                                                                                                                                                                                                                                               >  8:1288/1‑149 (MQ=255)
                                                                                                                                                           tttGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCGGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAata                                                                                                                                                                                                                                              <  8:630730/148‑1 (MQ=255)
                                                                                                                                                           aTTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAata                                                                                                                                                                                                                                              <  7:854270/149‑1 (MQ=255)
                                                                                                                                                           aTTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAata                                                                                                                                                                                                                                              >  7:343630/1‑149 (MQ=255)
                                                                                                                                                           aTTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAata                                                                                                                                                                                                                                              >  7:1346725/1‑149 (MQ=255)
                                                                                                                                                             tGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAAtata                                                                                                                                                                                                                                            <  7:741688/149‑1 (MQ=255)
                                                                                                                                                               tCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                                                                                                                                                                                          <  7:529240/149‑1 (MQ=255)
                                                                                                                                                               tCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                                                                                                                                                                                          >  5:115969/1‑149 (MQ=255)
                                                                                                                                                                  aGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCg                                                                                                                                                                                                                                       >  7:1253875/1‑149 (MQ=255)
                                                                                                                                                                    cGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTc                                                                                                                                                                                                                                     <  7:839430/149‑1 (MQ=255)
                                                                                                                                                                     gAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAACAGCAAAATATAAACCGtac                                                                                                                                                                                                                                    <  6:1016192/149‑3 (MQ=255)
                                                                                                                                                                      aaTCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                                                                                                                                                                                                   <  7:439964/149‑1 (MQ=255)
                                                                                                                                                                          gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaa                                                                                                                                                                                                                               >  6:38694/1‑149 (MQ=255)
                                                                                                                                                                             aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTa                                                                                                                                                                                                                                                            <  6:1201778/117‑1 (MQ=255)
                                                                                                                                                                             aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTa                                                                                                                                                                                                                                                            >  5:1201778/1‑117 (MQ=255)
                                                                                                                                                                             aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                                                                                                                                                            <  7:342451/149‑1 (MQ=255)
                                                                                                                                                                             aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                                                                                                                                                            <  7:1006299/149‑1 (MQ=255)
                                                                                                                                                                             aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                                                                                                                                                            <  5:694300/149‑1 (MQ=255)
                                                                                                                                                                             aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                                                                                                                                                            >  6:105297/1‑149 (MQ=255)
                                                                                                                                                                               aTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                                                                                                                                                                                                          <  5:1467250/149‑1 (MQ=255)
                                                                                                                                                                               aTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                                                                                                                                                                                                          >  8:1010796/1‑149 (MQ=255)
                                                                                                                                                                               aTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                                                                                                                                                                                                          <  5:54697/149‑1 (MQ=255)
                                                                                                                                                                                 tGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                                                                                                                                                                                        <  5:1355581/149‑1 (MQ=255)
                                                                                                                                                                                      ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCacca‑ac                                                                                                                                                                                                                  <  6:371422/149‑1 (MQ=255)
                                                                                                                                                                                           gAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAAC‑GGCTCTCATAGacagag                                                                                                                                                                                                                                                            >  8:884795/1‑98 (MQ=37)
                                                                                                                                                                                           gAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACaa                                                                                                                                                                                                             >  8:1253960/1‑149 (MQ=255)
                                                                                                                                                                                            aaTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACaaa                                                                                                                                                                                                            <  8:194141/149‑1 (MQ=255)
                                                                                                                                                                                              tGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAAcc                                                                                                                                                                                                          >  5:593630/1‑149 (MQ=255)
                                                                                                                                                                                               gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCa                                                                                                                                                                                                         <  8:686900/149‑1 (MQ=255)
                                                                                                                                                                                                ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaa                                                                                                                                                                                                                             <  6:105949/129‑1 (MQ=255)
                                                                                                                                                                                                ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAg                                                                                                                                                                                                        >  5:1071899/1‑149 (MQ=255)
                                                                                                                                                                                                 tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGa                                                                                                                                                                                                       <  7:1224600/149‑1 (MQ=255)
                                                                                                                                                                                                 tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGa                                                                                                                                                                                                       <  7:459381/149‑1 (MQ=255)
                                                                                                                                                                                                     tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACaaa                                                                                                                                                                                                   >  7:823001/1‑149 (MQ=255)
                                                                                                                                                                                                            aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACcgag                                                                                                                                                                                            >  5:716844/1‑149 (MQ=255)
                                                                                                                                                                                                              gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACcgagcg                                                                                                                                                                                          >  7:1447863/1‑149 (MQ=255)
                                                                                                                                                                                                               cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACcgagcga                                                                                                                                                                                         <  7:412139/149‑1 (MQ=255)
                                                                                                                                                                                                                 gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGAAACTTCCCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAAGATAAACCGTCCCCAAAAAAGCCACCA‑ACCCCAAACCAGACAAATCCCAAGGgaga                                                                                                                                                                                       >  8:1032379/1‑149 (MQ=255)
                                                                                                                                                                                                                   aCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAATAAAGCCACCA‑ACCACGAACCAGACAAAGACCGAGCGAGAAt                                                                                                                                                                                     <  6:567344/149‑1 (MQ=255)
                                                                                                                                                                                                                    cAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATa                                                                                                                                                                                    >  7:800730/1‑149 (MQ=255)
                                                                                                                                                                                                                       ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAAt                                                                                                                                                                                                                                                       >  3:27135/1‑80 (MQ=255)
                                                                                                                                                                                                                        gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACata                                                                                                                                                                                >  5:275361/1‑149 (MQ=255)
                                                                                                                                                                                                                         tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACcgagc                                                                                                                                                                                           <  8:374408/137‑1 (MQ=255)
                                                                                                                                                                                                                           gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCa                                                                                                                                                                                                                                                   <  1:27135/80‑1 (MQ=255)
                                                                                                                                                                                                                           gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATAc                                                                                                                                                                             <  8:823001/149‑1 (MQ=255)
                                                                                                                                                                                                                             tttGAGGCGAAGGTAAGCCGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATAccc                                                                                                                                                                           >  6:164565/1‑149 (MQ=255)
                                                                                                                                                                                                                                 aGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccaca                                                                                                                                                                                                                      >  1:52923/1‑102 (MQ=255)
                                                                                                                                                                                                                                   gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACt                                                                                                                                                                     >  7:799443/1‑149 (MQ=255)
                                                                                                                                                                                                                                       aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGaa                                                                                                                                                                 <  7:1253960/149‑1 (MQ=255)
                                                                                                                                                                                                                                               tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTccc                                                                                                                                                         >  7:1091612/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      ttCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATaaa                                                                                                                                                  >  8:1072079/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      ttCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATaaa                                                                                                                                                  >  6:1132466/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      ttCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATaaa                                                                                                                                                  >  6:800954/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           tGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtt                                                                                                                                             <  7:1072079/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt                                                                                                                                           <  8:837425/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               tATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTc                                                                                                                                          >  7:608024/1‑148 (MQ=255)
                                                                                                                                                                                                                                                               tATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTc                                                                                                                                          <  8:608024/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                aTGAGAGTCGTTTTAAC‑GGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCtt                                                                                                                                        >  7:244688/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                aTGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCtt                                                                                                                                        >  5:1300220/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  gagaGTCGTTTTAAC‑GGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaa                                                                                                                                      <  8:244688/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gagaGCCGTTAAAAC‑GTCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaa                                                                                                                                      <  7:1141776/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gagCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaaga                                                                                                                                    >  6:3413/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     agCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTACCGc                                                                                                                                                       <  8:232685/129‑1 (MQ=255)
                                                                                                                                                                                                                                                                     agCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTACCGc                                                                                                                                                       >  7:232685/1‑129 (MQ=255)
                                                                                                                                                                                                                                                                       ccGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagataa                                                                                                                                 <  6:761543/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       ccGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGAAGGTATACTGCAAATATTTCCTCCCGACTCAAg                       >  2:125327/1‑255 (MQ=255)
                                                                                                                                                                                                                                                                        cGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagataag                                                                                                                                <  5:378138/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cttAAAACAGACACTCATAAAATATACTAATAGCAAAAAAAAAACCGTCCCCAAAAAAACCACCA‑ACCACAAACCAGNAAAAAACCGAACGAGAATACATATACCCAAAACTCGAACAATTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGAAGGTATACTGCAAATATTTCCTCCCGACTCAAGCATCCAGTCATACTTATTGaa  <  3:69507/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                          ttAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGaa                                                                                                                              >  8:49597/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          ttAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGaa                                                                                                                              <  5:3413/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aaaaC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGaata                                                                                                                            >  5:224758/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                
CACAGGGCAGCGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAC‑GACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCA‑ACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGAAGGTATACTGCAAATATTTCCTCCCGACTCAAGCATCCAGTCATACTTATTGAA  >  NC_000913/3171589‑3172127

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑