Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,977,360 G→T P162Q (CCG→CAG)  flhC ← DNA‑binding transcriptional dual regulator FlhC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,977,3600GT97.9% 157.4 / ‑6.4 47P162Q (CCG→CAG) flhCDNA‑binding transcriptional dual regulator FlhC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (24/22);  minor base A (0/1);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCACTGTTCAAGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTC  >  NC_000913/1977116‑1977593
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         
cAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCAGGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATTTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGAc                                                                                                                                                                                                                                     <  1:281884/251‑1 (MQ=255)
  ggTAACATAATAAGATAAGCCCGACAACCTCCTTCCACTGTTGACGATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGTGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGGCTTTTTACTGCCCGGGATGGCTGTTGACAt                                                                                                                                                                                                                                   <  2:13494/251‑1 (MQ=255)
     aaCCTAATAAGATAAGCCCGACATCATCCTTCCACTGTTGACCATGACAGGATGGTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCGGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATGCTGGGAAAGTTTACGTCGTTTTACTGCCCGGGATGGCTGTTGACATaa                                                                                                                                                                                                                                 <  2:213641/250‑1 (MQ=255)
                  aagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGc                                                                                                                                                                                                                   >  2:340775/1‑251 (MQ=255)
                      aCGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGcc                                                                                                                                                                                                               >  1:10135/1‑251 (MQ=255)
                       cGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCa                                                                                                                                                                                                              <  2:366081/251‑1 (MQ=255)
                        gACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCa                                                                                                                                                                                                              >  1:104745/1‑250 (MQ=255)
                                      aCTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAAGTCGTTACCGCTGATGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGGTTACGTCTTTTTACTGCCAGGGATGGATGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCg                                                                                                                                                                                               >  1:31597/1‑251 (MQ=255)
                                         gTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTgg                                                                                                                                                                                            <  1:257045/251‑1 (MQ=255)
                                            gACCATGACAGGATGTTCCGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGACTGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCTTCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAACGTTTACGTCTTTTTACTGCTCGGGATGGCTGTTGACATAAGATGCAGGCAAATCTTCCAACAGGCTGGTGAGCGTGGGt                                                                                                                                                                                          >  1:62314/1‑250 (MQ=255)
                                                            tcaggcgtcaggcgTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGccgccg                                                                                                                                                                          <  1:464434/250‑1 (MQ=255)
                                                                 cgtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCgcagc                                                                                                                                                                      >  1:82669/1‑249 (MQ=255)
                                                                   tcAGGCGTTACCGCGCGATTGGGGCAACAAAAAGTAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAACGCTGACAACAGGCTGGTGAGCGTGGGAAATAAAATTGCCGCCGCAGCCGtt                                                                                                                                                                  >  2:366565/1‑251 (MQ=255)
                                                                         gTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACTGCTGCTGGAATGTTGCGCCTCCCCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCACCAGTTGTGGGATAATATCGGCAGGCTTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACCTAAGCTGCAGGCAAAGCTGCCAACAGGCTGGCAACCGTGGGACATAAAATTGCCGCCGCAGCAGGTTCAGCt                                                                                                                                                            >  2:302826/1‑251 (MQ=255)
                                                                                    aTTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGtt                                                                                                                                                    <  2:278050/248‑1 (MQ=255)
                                                                                    aTTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGtt                                                                                                                                                    >  1:278050/1‑248 (MQ=255)
                                                                                      tGGGGCAAAAAAAAGCAGCGGTCCGTAGTTCCCGTTGATCGCATGTTGCGCATCAACGTATCAGTTAAACAGCCTGTACACTCTGTTCATCCAGCAGTTGTGGGATAATCTCGGCAGGATTCTGGGAAAGTTTCCGTCTTTTTACCGCCCGGGATGGCTGTTGACATAACCTGCAGGCAAAGCTGACAACTGGCTGGTGGGCGTGGGTAATAACATTGCCGCCGCAGCAGTTGCAGCTGGAACGTTGCAGt                                                                                                                                               >  2:232197/1‑251 (MQ=255)
                                                                                       ggggCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTa                                                                                                                                              <  1:44697/251‑1 (MQ=255)
                                                                                                  aaGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTc                                                                                                                                    >  1:229486/1‑250 (MQ=255)
                                                                                                      agcGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTcaaca                                                                                                                                >  1:382035/1‑250 (MQ=255)
                                                                                                       gcGGTACTACGTTACCGCTGTTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACaaa                                                                                                                              <  2:439471/251‑1 (MQ=255)
                                                                                                          gTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCg                                                                                                                           >  1:272432/1‑251 (MQ=255)
                                                                                                                    aCCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGt                                                                                                                  <  1:243219/250‑1 (MQ=255)
                                                                                                                                        ccTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGTATAATATCGGCAGGATTCTTGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCcag                                                                                             <  2:547170/251‑1 (MQ=255)
                                                                                                                                        ccTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTGACTGCCCGTGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCcag                                                                                             <  1:83832/251‑1 (MQ=255)
                                                                                                                                                   cAGTTAAACAGCCTTTACTCTCTGTTCATCCAGCAGTTGTGGTATAATATCGTCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCctt                                                                                  <  1:505969/251‑1 (MQ=255)
                                                                                                                                                   cAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCGGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGTATGGCTTTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCct                                                                                   <  1:313849/250‑1 (MQ=255)
                                                                                                                                                    aGTTAACCAGCCTGTCCTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGCAAGTTTCCGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGATGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTgg                                                                                                                                                                                            >  2:55987/1‑144 (MQ=255)
                                                                                                                                                    aGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGTTGGTGAGCGTgg                                                                                                                                                                                            <  1:55987/144‑1 (MQ=255)
                                                                                                                                                    aGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCcagc                                                                                            >  1:320052/1‑240 (MQ=255)
                                                                                                                                                    aGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCcagc                                                                                            <  2:320052/240‑1 (MQ=255)
                                                                                                                                                      ttAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCcttctt                                                                               >  1:167429/1‑251 (MQ=255)
                                                                                                                                                       tAAACAGCCTGTACTCTCTGATCATCCAGCAGGTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTATTTATACTGACCGGGCGGCCTGTTGACATAAGCTGCAGGCAAAACTGCCAACAGGATGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGAAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCAAGTCACGGGGTAATGCCAGCAGAGGTCCtacttc                                                                              >  1:482773/1‑251 (MQ=255)
                                                                                                                                                                tGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTgg                                                                     >  1:428403/1‑251 (MQ=255)
                                                                                                                                                                   aCTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCa                                                                  <  2:180621/251‑1 (MQ=255)
                                                                                                                                                                      ctctGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTTCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTGCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCa                                                                  <  1:529545/248‑1 (MQ=255)
                                                                                                                                                                                    cagTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCACTGTTCAAGGTATAAAc                                                 >  1:522523/1‑251 (MQ=255)
                                                                                                                                                                                    cagTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCACTGTTCAAGGTATAAAc                                                 >  1:536370/1‑251 (MQ=255)
                                                                                                                                                                                        ggtgGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCACTGTTCAAGGTATAAACGGTa                                             <  2:10135/250‑1 (MQ=255)
                                                                                                                                                                                             gATAATAGCGGCAGGATTCGGGGAAAGTTAACGTCTTTTTACTGCCCGGGATGGCTGTTGACAAAAGCTGCATGCAAAGCTGCCAACAGGCTGGTGATCGTGGGTAAGAAAATGGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCTCTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGGGGGCCTTCTTCGGCGTGTGGGCACTGTTCAAGGTATAAACGGTAGGCtt                                        <  2:87920/251‑1 (MQ=255)
                                                                                                                                                                                                 atatCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCCCTGTTCAAGGGATAAACGGTAGGCTTTGAt                                    >  2:64126/1‑251 (MQ=255)
                                                                                                                                                                                                                  ggAAAGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCACTGTTCAAGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCAt                   >  1:173594/1‑251 (MQ=255)
                                                                                                                                                                                                                     aaGTTTCCGTCTGTTGGCTGCCCGGGATGGCAGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCTCCAATGTCCAGGCACGGGTTAAGGCCAGCAGTGGTCCTTCTTCTGCTTGGGGGCACTGTTCAAGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTac                <  1:302826/251‑1 (MQ=255)
                                                                                                                                                                                                                     aaGTTTACGTCTTTTTACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCACTGTTCAAGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTac                <  1:278916/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ttttACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCACTGTTCACGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTCCACAAACCGGttt    >  2:303010/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tttACTGCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCACTGTTCAAGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGtttt   >  1:512651/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tacggCCCGGGATGGCTGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCACTGTTCAAGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTACTCAAACCGGTTTTc  <  1:421150/246‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         
CAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACATAAGCTGCAGGCAAAGCTGCCAACAGGCTGGTGAGCGTGGGTAATAAAATTGCCGCCGCAGCAGTTGCAGCTGGAAAGTTGCAGTAATCCACTTTCAACAAACCGCACCAATGTCCAGGCACGGGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCACTGTTCAAGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTC  >  NC_000913/1977116‑1977593

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