Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 875,395 G→T 9.6% A21S (GCA→TCA)  yliF → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913875,3950GT 9.6% 213.9 / 4.8 84A21S (GCA→TCA) yliFputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (47/27);  new base T (2/8);  total (49/35)
Fisher's exact test for biased strand distribution p-value = 1.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACG  >  NC_000913/875246‑875544
                                                                                                                                                     |                                                                                                                                                     
tgtTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTcag                                                                                                                                                       >  2:780851/1‑150 (MQ=255)
        cTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGca                                                                                                                                               >  2:1651085/1‑149 (MQ=255)
        cTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCGGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGcc                                                                                                                                               >  2:1686187/1‑150 (MQ=255)
          cAGGGGTATCTGTACCAGAAAACCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTTAAAGAGTAGTTTACTTATGTCCAGAATCAAAAATTTAGTACTAAAAGTCATTATGCTAAATTTTATAATAATTTCAGCACTTGACtt                                                                                                                                             >  1:360959/1‑149 (MQ=255)
           aGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGTGATAAAAGTGTAAAATTAGTAAAAGAGTAGTTTACGTATGTACAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGc                                                                                                                                            <  2:970578/150‑1 (MQ=255)
           aGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGc                                                                                                                                            >  2:413440/1‑150 (MQ=255)
             gggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGTGATAAAAGGGGAAATTTAGTAAAAGAGTAGTTTACGAATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCgg                                                                                                                                          <  2:1675779/150‑1 (MQ=255)
             gggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCgg                                                                                                                                          <  2:837667/150‑1 (MQ=255)
             gggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAAATAGGAAAAGAGTAGTTTACCTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTCCTGATTTTTATCATGATATCATCAGTTGCCTGCgg                                                                                                                                          >  1:1497303/1‑150 (MQ=255)
              ggTATCTGTACCAGAATCCCATGCCATTCTCCGCCTGGGATAAAAGTGGTAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGAACTTACAGTCAGTCTGCTGATTTTTTTCATGATTTCAGCAGTTGCCTGCggg                                                                                                                                         <  1:490894/150‑1 (MQ=255)
              ggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGctgt                                                                                                                                         >  2:1276473/1‑147 (MQ=255)
              ggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCggg                                                                                                                                         <  2:1127759/150‑1 (MQ=255)
              ggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCggg                                                                                                                                         <  1:1608387/150‑1 (MQ=255)
              ggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATAAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCggg                                                                                                                                         >  1:1211783/1‑150 (MQ=255)
                 aTCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGAta                                                                                                                                      >  2:1595831/1‑149 (MQ=255)
                 aTCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATc                                                                                                                                      <  1:1145221/150‑1 (MQ=255)
                  tCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTTGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCt                                                                                                                                     <  2:734261/150‑1 (MQ=255)
                  tCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCt                                                                                                                                     <  2:2441/150‑1 (MQ=255)
                  tCTGTACCAGAAAACCATGCCATGCTCCGCCTGTGATAAAAGTGGAAAATTAGTAAAAGAGGAGTTTCCGTATGTCCAGATTCAATAAGTTCGTACTTACAGTCAGGCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCt                                                                                                                                     <  1:308385/150‑1 (MQ=255)
                    tgaaCCAGAAACCCATGACAATCTCAGCCTTGGATAAAAGTGTAAAATTAGTAAAAAAGTAGTTTAATTATGTCCAGAATCAATAAGTTCCTACTTACAGTCACTCTGCAGATTTTTAACATGATATCAGCAGTTGCCTAAGGGATCTaa                                                                                                                                   >  2:642875/4‑149 (MQ=255)
                    tGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTac                                                                                                                                   <  2:345818/150‑1 (MQ=255)
                     gTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTaca                                                                                                                                  <  1:833950/150‑1 (MQ=255)
                     gTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGACAGACTGATGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTaca                                                                                                                                  >  1:731986/1‑150 (MQ=255)
                      tACCAGAAACCCATGACATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTacac                                                                                                                                 >  2:1906442/1‑150 (MQ=255)
                       aCCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTacaca                                                                                                                                >  2:859237/1‑149 (MQ=255)
                         cAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGAACTACACTCa                                                                                                                              >  2:1878354/1‑150 (MQ=255)
                          aGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCaa                                                                                                                             >  2:843537/1‑150 (MQ=255)
                           gAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCaaa                                                                                                                            >  2:427169/1‑150 (MQ=255)
                               cccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGt                                                                                                                        <  1:1659425/150‑1 (MQ=255)
                                ccATGCCCTTCTCCGCCTGGCATAACAGTCGAAAAAAAGAACAATAGTCCTTTACTTAAGGCCAGAATCAATACCTTCGTACTTTCAGTCAGTCTGCTGATTTTTTTTTTTATTGATTCTGTGTCCAGCGTGATCTACACTCAAATGGTa                                                                                                                       <  2:101407/150‑1 (MQ=255)
                                   tGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAg                                                                                                                    <  1:659736/150‑1 (MQ=255)
                                      cATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGaa                                                                                                                 >  2:72184/1‑150 (MQ=255)
                                            ccGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGgtg                                                                                                           <  1:866812/150‑1 (MQ=255)
                                             cGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAAACAATAAGTTCGTACTTACAGTCAGTCTGCTGATTATTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAATGAACGGctgt                                                                                                          >  2:101551/1‑150 (MQ=255)
                                               ccTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGtat                                                                                                        >  2:984475/1‑150 (MQ=255)
                                                cttggATAAAAGTGGAAAATTAGTAAAAGAGAAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGAATTTTATCATGATTTCATCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACTGCTGtata                                                                                                       >  2:1888574/4‑150 (MQ=255)
                                                cTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGtata                                                                                                       >  2:56103/1‑150 (MQ=255)
                                                  gggATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGc                                                                                                     >  2:918929/1‑150 (MQ=255)
                                                    gATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTAGTTACAGTCAGTCTGCTGATTTTTATGATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCt                                                                                                   <  1:784454/150‑1 (MQ=255)
                                                    gATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATTTTTATTTATTCAGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCt                                                                                                   <  2:740718/150‑1 (MQ=255)
                                                      tAAAATTGCAAAAATAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTTTTTTTATTTATTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGa                                                                                                 <  2:386377/150‑1 (MQ=255)
                                                      tAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTACAGCAGTTGCCTGCGGGTTCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGa                                                                                                 >  2:702249/1‑150 (MQ=255)
                                                        aaaGTGGAAAATTAGTAAAAGAGTAGATTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATATTTATCATGATTTCAGCAGTTGCCTGCCGGATCTACACTCAAAAGGTAAAGGAACGGGTGTATAGCCTGaaa                                                                                               >  2:1621572/1‑150 (MQ=255)
                                                         aaGTGGAAAATTAGTAAAAGAGTAGGTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAAc                                                                                              <  1:343375/150‑1 (MQ=255)
                                                           gTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAg                                                                                            >  2:1508279/1‑150 (MQ=255)
                                                             ggAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGAATAGCCTGAAACAGTc                                                                                          >  1:1297559/1‑150 (MQ=255)
                                                             ggAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTATGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTc                                                                                          >  2:17445/1‑150 (MQ=255)
                                                              gAAAATTGGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTcc                                                                                         >  1:1492853/1‑150 (MQ=255)
                                                                    tattAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTTTTTTTATTTATTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTAtt                                                                                   <  2:354789/147‑1 (MQ=255)
                                                                       tAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGAt                                                                                >  2:683621/1‑150 (MQ=255)
                                                                            gagTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTTCTTTTTTTTTTTTTTATTTATTCTGTTTCCTGCGTGATCTAAACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGc                                                                           <  2:1644696/150‑1 (MQ=255)
                                                                                      cGTATGTCCAGAATCAATAAGTTAGTACTTACAGTATGTCTGCTGATTTTTATCATTATTTAATCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTGTAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTg                                                                 <  2:302307/150‑1 (MQ=255)
                                                                                        tATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTTATTTTTTTTTTTATTTATTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGc                                                               <  2:1214703/150‑1 (MQ=255)
                                                                                          tGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCaa                                                             <  1:1026383/150‑1 (MQ=255)
                                                                                           gtacAGAATCAATAAGTTCGTACTTACAGTCAGTCTCCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGTAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTTGCaaa                                                            <  1:566311/147‑1 (MQ=255)
                                                                                              cAGAATCAATAAGTTCGTACTTACAGTCAGTCTGATGATTTTTATCATGATTTCAGCAGTTGCCTGCTGGATCTAAACTCAAATGGTAAAGGAAAGGGTGTATATCCTGAAACAGTCCGGTATTGCTACTGCTTTTGCGGTGGCAAatat                                                         >  2:1626514/1‑150 (MQ=255)
                                                                                                gAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTg                                                       >  1:1150457/1‑150 (MQ=255)
                                                                                                     aaTAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGaa                                                  >  2:1330673/1‑150 (MQ=255)
                                                                                                          gTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCTGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCg                                             >  2:1600551/1‑150 (MQ=255)
                                                                                                          gTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCg                                             >  1:603160/1‑150 (MQ=255)
                                                                                                          gTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCg                                             >  2:831994/1‑150 (MQ=255)
                                                                                                              gTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGt                                         >  2:1901168/1‑150 (MQ=255)
                                                                                                                 cTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGc                                      >  2:995266/1‑150 (MQ=255)
                                                                                                                 cTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGc                                      >  2:1521904/1‑150 (MQ=255)
                                                                                                                     cagttagtCTGCTGATTTTTATGATGATTTCAGCAGTTGCCTGCGGGATCTTCACTCAAATGGTAATGGAAGGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTgg                                  <  1:243164/150‑1 (MQ=255)
                                                                                                                     cagtcagtCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTgg                                  <  1:413440/150‑1 (MQ=255)
                                                                                                                       gtcagtcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTg                                                                 >  2:32798/1‑117 (MQ=255)
                                                                                                                       gtcagtcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTg                                                                 <  1:32798/117‑1 (MQ=255)
                                                                                                                       gtcagtcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCa                                >  2:384868/1‑150 (MQ=255)
                                                                                                                       gtcagtcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCa                                >  2:583126/1‑150 (MQ=255)
                                                                                                                         cagtcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAAt                              >  2:451850/1‑150 (MQ=255)
                                                                                                                           gtcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTg                            >  2:1586511/1‑150 (MQ=255)
                                                                                                                           gtcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTg                            >  2:1229791/1‑150 (MQ=255)
                                                                                                                            tcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGa                           <  1:683621/150‑1 (MQ=255)
                                                                                                                             ctgctgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGAt                          >  1:1347460/1‑150 (MQ=255)
                                                                                                                              tgctgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGTCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATc                         >  2:1355250/1‑150 (MQ=255)
                                                                                                                              tgctgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATc                         >  2:1355102/1‑150 (MQ=255)
                                                                                                                               gctgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCt                        <  1:72184/150‑1 (MQ=255)
                                                                                                                                 tgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTa                      >  2:517591/1‑150 (MQ=255)
                                                                                                                                  gATTTTTATCATGTTTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTTGCAATTGATCTTAt                     <  1:642875/150‑1 (MQ=255)
                                                                                                                                    tttttATTTTGATTTCTTCTGTTTCCTTCGTGATCTAGACTCATATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCa                   <  2:1828563/150‑1 (MQ=255)
                                                                                                                                           cATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGca            >  2:756651/1‑149 (MQ=255)
                                                                                                                                            ttGATTTCTGCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTa           <  2:1224660/149‑1 (MQ=255)
                                                                                                                                                 ttCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATcc      <  1:978438/150‑1 (MQ=255)
                                                                                                                                                  tCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATccc     >  2:1032971/1‑150 (MQ=255)
                                                                                                                                                   cagcagTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTAGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACACTAAATCCCa    >  1:1516883/1‑150 (MQ=255)
                                                                                                                                                     gcagTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACg  >  1:1706119/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACG  >  NC_000913/875246‑875544

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: