Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 875,393 C→A 8.6% S20* (TCA→TAA)  yliF → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913875,3930CA 8.6% 277.7 / 4.6 89S20* (TCA→TAA) yliFputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (46/33);  new base A (2/7);  total (49/40)
Fisher's exact test for biased strand distribution p-value = 7.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCC  >  NC_000913/875248‑875540
                                                                                                                                                 |                                                                                                                                                   
ttCAGGCGCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTCGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATTAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTcagc                                                                                                                                                  <  2:267779/149‑1 (MQ=255)
  cAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTcagca                                                                                                                                                 <  4:311225/148‑1 (MQ=255)
    ggCTCAGGTGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAGTGGAAAATTAGTAAAAGGGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCAGATTTTTATCATGATTTCAGCAGTTg                                                                                                                                             >  2:39890/1‑149 (MQ=255)
    ggCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGtt                                                                                                                                              <  5:119654/149‑1 (MQ=255)
           gggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCg                                                                                                                                       <  6:47251/149‑1 (MQ=255)
           gggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCg                                                                                                                                       >  2:192211/1‑149 (MQ=255)
           gggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCg                                                                                                                                       >  2:197800/1‑149 (MQ=255)
                   gTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTAAGCAGTTGCCTGAGGGAtata                                                                                                                                >  6:331873/1‑148 (MQ=255)
                    tACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATc                                                                                                                                  <  1:192211/145‑1 (MQ=255)
                    tACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATc                                                                                                                                  <  1:197800/145‑1 (MQ=255)
                         gAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCa                                                                                                                          >  4:352425/1‑148 (MQ=255)
                           aaCCCATGCCATTCTCCGCCTGTGATAAAAGTGGAAAATTAGTAAAAGAGTACTTTACTTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTTTTTTTAATTATTCTGTTTCCTGCGTGATCTACACTCAAAt                                                                                                                       <  2:59105/149‑1 (MQ=255)
                           aaCCCACGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTTTTTTTATTGATTCTGTTTCCTGCGTGATCTACACTCAAAt                                                                                                                       <  2:59104/149‑1 (MQ=255)
                            aCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATg                                                                                                                      >  5:221999/1‑149 (MQ=255)
                            aCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCATGCGGGATCTACACTCAAATg                                                                                                                      >  3:68678/1‑149 (MQ=255)
                              ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGt                                                                                                                    <  5:331873/149‑1 (MQ=255)
                                    cATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTca                                                                                                                                                    <  3:218615/111‑1 (MQ=255)
                                    cATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTca                                                                                                                                                    >  4:218615/1‑111 (MQ=255)
                                     aTTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTTCAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGaa                                                                                                             >  1:217915/1‑149 (MQ=255)
                                      ttCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCTTACTTACAGTCAGTCTGCAGATTTTTATCATGATTTGAGAAGTTGCCTTCGGGATCTACACTCAAATGGTAAAGGAAc                                                                                                            >  5:76506/1‑149 (MQ=255)
                                         tcCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCACAATCAATAAGTTCTTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTAAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACggc                                                                                                         >  7:294569/1‑148 (MQ=255)
                                           cGCCTGGGATAAAAGTGGAAATTTAGTAAAAGCTTAGTTTACGTATGTCCAGAATCGATAAGTTCGTACTTACAGTCAGTCTCCTGATTTTTTTTTTTATTGATTCTGTGTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGgt                                                                                                        <  6:55157/148‑1 (MQ=255)
                                           cGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACgggag                                                                                                       >  2:141768/1‑147 (MQ=255)
                                           cGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGgt                                                                                                        <  5:180715/148‑1 (MQ=255)
                                                gggATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAg                                                                                                  >  2:172496/1‑149 (MQ=255)
                                                 ggATAAAAGTTGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCAGGATTTCAGCCGTTGCCTGCGGGATCTACACTCAAATGGTAATGGAACGGGTGTATAGc                                                                                                 >  8:376703/1‑149 (MQ=255)
                                                  gATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCAAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGcc                                                                                                >  8:324140/1‑149 (MQ=255)
                                                     aaaaGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTg                                                                                              >  6:321821/1‑148 (MQ=255)
                                                       aaGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGaaa                                                                                           <  6:221999/149‑1 (MQ=255)
                                                                  tAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTAt                                                                                >  5:379846/1‑149 (MQ=255)
                                                                    gTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTAtt                                                                               <  7:414332/148‑1 (MQ=255)
                                                                    gTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTAtt                                                                               <  7:352949/148‑1 (MQ=255)
                                                                      aaaaGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGa                                                                             >  6:93441/1‑148 (MQ=255)
                                                                       aaaGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTc                                                                                                                                                     >  2:249672/1‑75 (MQ=255)
                                                                       aaaGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTc                                                                                                                                                     <  1:249672/75‑1 (MQ=255)
                                                                       aaaGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAg                                                                                                                >  4:8639/1‑112 (MQ=255)
                                                                       aaaGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAg                                                                                                                <  3:8639/112‑1 (MQ=255)
                                                                       aaaGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCAGCTTATTTTTATCATGAATTCAGCAGTTGCCTGCGGGATCTACAATCAAATGGTAATGGAACGGTTGTATAGCCTGAAACAGTCGGATATAGATa                                                                           >  3:351626/1‑149 (MQ=255)
                                                                               gtTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCtttt                                                                   >  3:54683/1‑149 (MQ=255)
                                                                                tttACGTATGTCCAGACTCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACgg                                                                                                          <  5:136439/109‑1 (MQ=255)
                                                                                tttACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACgg                                                                                                          >  6:136439/1‑109 (MQ=255)
                                                                                tttACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCACTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTg                                                                  <  7:382232/149‑1 (MQ=255)
                                                                                tttACGTATGTCCAGAATCAATAAGTTCGTACTTAAAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTg                                                                  >  8:8989/1‑149 (MQ=255)
                                                                                  tACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGc                                                                 <  3:71445/148‑1 (MQ=255)
                                                                                  tACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCg                                                                >  5:57062/1‑149 (MQ=255)
                                                                                       aTGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCTAATGGTAAAGGAACGCGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTgg                                                            >  1:335066/1‑148 (MQ=255)
                                                                                        tGTCCAGAATCAATAAGTTCGCACTTACAGTCAGTCTGCTGATTTTTTTCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGc                                                           <  3:1991/148‑1 (MQ=255)
                                                                                           ccAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCt                                                                                               >  4:299176/1‑109 (MQ=255)
                                                                                           ccAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCt                                                                                               <  3:299176/109‑1 (MQ=255)
                                                                                           ccAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCaaa                                                        >  2:153136/1‑148 (MQ=255)
                                                                                            cAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAata                                                      >  1:236554/1‑149 (MQ=255)
                                                                                            cAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAata                                                      >  6:254024/1‑149 (MQ=255)
                                                                                             aGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAatat                                                     >  8:31711/1‑149 (MQ=255)
                                                                                               aaTCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTg                                                   <  5:229718/149‑1 (MQ=255)
                                                                                               aaTCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGAAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGATTTTGCGGTGGCAAATATTg                                                   >  3:335523/1‑149 (MQ=255)
                                                                                                tctaaTAAGTTAGGACTTAAAGTCATTCTGCTGTTTTTTTTCTTTATTTCATATGTTTCCTGCGGGATCTACACTCAAATTGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATAATGCTTTTGCGGTGGCAAATAtt                                                    <  6:29892/144‑1 (MQ=255)
                                                                                                 tCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCt                                                 <  1:172496/149‑1 (MQ=255)
                                                                                                 tCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCt                                                 <  1:141768/149‑1 (MQ=255)
                                                                                                 tCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTAAACTAAAATGGTAAAGGAACGGGTGTATATCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCt                                                 >  1:207933/1‑149 (MQ=255)
                                                                                                   aaTAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCt                                                 >  8:344341/1‑147 (MQ=255)
                                                                                                   aaTAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCt                                                 <  7:344341/147‑1 (MQ=255)
                                                                                                    aTAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTa                                                                                                                                >  5:386119/1‑67 (MQ=255)
                                                                                                    aTAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTa                                                                                                                                <  6:386119/67‑1 (MQ=255)
                                                                                                    aTAAGTTCGTACTTAAAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGTTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATAGATACTGCTTTTGCGGTGGCAAATAtttaga                                                >  5:29892/1‑143 (MQ=255)
                                                                                                       aGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAatat                                           <  3:352425/149‑1 (MQ=255)
                                                                                                        gTTCGTACTTACAGTCAGTCTGCTGATTTTTTTTTTTAATGATTCTGTGTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATc                                          <  2:217915/149‑1 (MQ=255)
                                                                                                        gTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATc                                          <  5:321821/149‑1 (MQ=255)
                                                                                                         ttCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCg                                         >  7:27418/1‑149 (MQ=255)
                                                                                                                ttACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGc                                                                 >  4:81949/1‑118 (MQ=255)
                                                                                                                ttACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGc                                                                 <  3:81949/118‑1 (MQ=255)
                                                                                                                ttACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGc                                  >  8:108634/1‑149 (MQ=255)
                                                                                                                ttACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGc                                  >  1:275169/1‑149 (MQ=255)
                                                                                                                  aCAGTCAGTCTGCTGATTTTTATTATTATTTAATCAGTTTCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGc                                  <  4:54683/147‑1 (MQ=255)
                                                                                                                  aCAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGt                                >  8:301699/1‑149 (MQ=255)
                                                                                                                  aCAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGt                                >  8:400564/1‑149 (MQ=255)
                                                                                                                   cagtcagtCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTg                               <  7:400564/149‑1 (MQ=255)
                                                                                                                   cagtcagtCTGATTATTTTTTTTTTGATTTCTGCTGTTTCTTGCGTGATCTACACTGATATGGTAAAGGAACGGGTGTATAGCGTGAAACTGTCCGTTATTGATACTGATTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGt                                <  6:76506/148‑1 (MQ=255)
                                                                                                                      tcagtcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGGGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGc                             >  7:222186/1‑148 (MQ=255)
                                                                                                                        agtcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAAt                          <  6:57062/149‑1 (MQ=255)
                                                                                                                             gctgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCAACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGAt                      >  8:151780/1‑148 (MQ=255)
                                                                                                                              ctgATTTTTTTTTTTATTTAATCTGTTGCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCt                    <  8:27418/149‑1 (MQ=255)
                                                                                                                              ctgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATc                     <  5:254024/148‑1 (MQ=255)
                                                                                                                               tgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGa                       <  5:114112/145‑1 (MQ=255)
                                                                                                                               tgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGa                       >  6:114112/1‑145 (MQ=255)
                                                                                                                                gATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAATGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGGGTGGCAATTGATCTTa                  >  1:139488/1‑149 (MQ=255)
                                                                                                                                gATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTa                  >  1:131715/1‑149 (MQ=255)
                                                                                                                                 aTTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTa                  <  7:108634/148‑1 (MQ=255)
                                                                                                                                 aTTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTAt                 <  7:31711/149‑1 (MQ=255)
                                                                                                                                 aTTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTAt                 >  5:97125/1‑149 (MQ=255)
                                                                                                                                  tttttATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATc                <  5:93441/149‑1 (MQ=255)
                                                                                                                                     ttATCTTGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCaa              <  8:222186/148‑1 (MQ=255)
                                                                                                                                       attaTGATTTAAGCAGTTTCCTGCGGGATATACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAaca            <  2:335066/145‑1 (MQ=255)
                                                                                                                                          aTGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGc         >  7:261411/1‑148 (MQ=255)
                                                                                                                                          aTGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGc         >  3:276679/1‑148 (MQ=255)
                                                                                                                                           tGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCt        >  8:298441/1‑148 (MQ=255)
                                                                                                                                                 cagcagTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATcc  >  3:347262/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCC  >  NC_000913/875248‑875540

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