Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 875,408 G→T 13.9% G25V (GGG→GTG)  yliF → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913875,4080GT13.9% 118.4 / 2.9 60G25V (GGG→GTG) yliFputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (28/22);  new base T (2/8);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 7.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAG  >  NC_000913/875261‑875547
                                                                                                                                                   |                                                                                                                                           
gTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCggg                                                                                                                                            <  1:200283/149‑1 (MQ=255)
  aTCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGAt                                                                                                                                          <  1:388949/149‑1 (MQ=255)
   tCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATc                                                                                                                                         <  2:11726/149‑1 (MQ=255)
     tGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGGCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGTATCTa                                                                                                                                       <  2:908780/149‑1 (MQ=255)
     tGTACCAGAAACCCATGCCATTCTACGACTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATTTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCATGCGGGATCTa                                                                                                                                       <  2:890401/149‑1 (MQ=255)
      ttACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTac                                                                                                                                      <  1:575729/148‑1 (MQ=255)
      gTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTac                                                                                                                                      >  2:265130/1‑149 (MQ=255)
          aaGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTc                                                                                                                                  <  1:1145693/148‑1 (MQ=255)
                    tGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCTTACTTACATTCAGTCTGCTGATATTTATCATGATTTCAGCAGATTCCTGCGGGATCTACACTCAAATGGTaaa                                                                                                                        >  2:124550/1‑149 (MQ=255)
                    tGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCATGCGGGATCTACACTCAAATGGTaaa                                                                                                                        >  2:613944/1‑149 (MQ=255)
                     gCCATTCTCCGCATGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCTTGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAg                                                                                                                       <  1:124550/149‑1 (MQ=255)
                       cATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGCATCTACACTCAAATGGTAAAGGa                                                                                                                     >  1:259624/1‑149 (MQ=255)
                         ttCTCCGCCTGGGATAAAAGTGGAAAATTACTAAAACAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTTTTTTTATTTATTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAAc                                                                                                                   <  2:54590/149‑1 (MQ=255)
                                 cTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTATGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTTATGATTTTTATCATGATTTCAGCAGTTGCCTGCGGTATCTACACTCAAATGGTAAAGGAACGGGTGtat                                                                                                           >  1:707304/1‑149 (MQ=255)
                                       tAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTg                                                                                                     >  1:253939/1‑149 (MQ=255)
                                           aGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTCCAGTCAGTCTGATTATTTTTTTTTTTATTTCTTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAAc                                                                                                 <  2:1007165/149‑1 (MQ=255)
                                           aGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCAAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGCTTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTAAACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAAc                                                                                                 >  1:879997/1‑149 (MQ=255)
                                              ggAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAAAGGGTGTATAGCCTGAAACAGt                                                                                              >  1:80712/1‑149 (MQ=255)
                                               gAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTc                                                                                             >  2:678455/1‑149 (MQ=255)
                                                  aaTTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGGGTATAGCCTGAAACAGTCCGt                                                                                          <  1:1018583/149‑1 (MQ=255)
                                                      aGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATGATGATTTAAGCTGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTAtt                                                                                      <  2:1095348/149‑1 (MQ=255)
                                                       gTAAAAGAGTAGTTTACGTATGTCAAGAATCAATAAGTTAGTACTTACAGTCAGTCTACTGATTTTTATAATGATTTCAGCAGTTGACTGCGTTATCTACACTCAAATGGTAAAGGAACGGGTGTATACCCTGAAACATTCAGTTAtta                                                                                     >  2:1088137/1‑148 (MQ=255)
                                                        tAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGa                                                                                    >  1:440785/1‑149 (MQ=255)
                                                          aaaGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATa                                                                                  >  2:846948/1‑149 (MQ=255)
                                                          aaaGAGTAGTTTACGTATGTCCAGAAGCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATa                                                                                  >  1:643720/1‑149 (MQ=255)
                                                                 agtATACGGATCTCCAGAATCAATAAGTTCGTACTTACAGGCAGTCTGCTGATTTTTTTTTTTATTTATTCAGTTGCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCttt                                                                           <  2:91039/149‑1 (MQ=255)
                                                                   tgtACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTg                                                                         <  1:16122/147‑1 (MQ=255)
                                                                         tATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTgg                                                                   >  1:407407/1‑149 (MQ=255)
                                                                          aTGTCCTGAATCAATAATTGCGTACTTTCAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATATACACTCAAATGGTAAAGGAACGGGTGTATAGCCGGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGc                                                                  <  1:1112511/149‑1 (MQ=255)
                                                                            gTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCaa                                                                >  1:927416/1‑149 (MQ=255)
                                                                             tCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAAAACTCAGTTATAGCTACTGCTTTTGAGGTGTCaaa                                                               >  1:221163/1‑149 (MQ=255)
                                                                                  aaTCAATAAGTTCGTACTTACAGTCAGTGTGCTGTTTTTTATAATTATTTCATATGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTg                                                          <  2:1147096/149‑1 (MQ=255)
                                                                                     cAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGAAAATATTGCTg                                                       >  1:1017068/1‑149 (MQ=255)
                                                                                      aaTAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATAATGCTTTTGAGGTGGCAAATATTGCtta                                                      >  2:730769/1‑147 (MQ=255)
                                                                                      aaTAAGTTAGTACTTCCAGTCAGGCTGCTGATTTTTATCATGATTTCAGCTGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGa                                                      <  2:707304/149‑1 (MQ=255)
                                                                                        tAAGTTCGTACTTACAGTCAGTCTGCTGAGTTTTATCATGATTTCCGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCTTTATTTATACTGCTTTTGCGGTGGCATATATTGCTGAat                                                    >  1:1116644/1‑149 (MQ=255)
                                                                                          aGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTGGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAatat                                                  <  1:276818/149‑1 (MQ=255)
                                                                                           gTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATACCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATc                                                 >  2:425655/1‑149 (MQ=255)
                                                                                            ttCGTACTTACAGTCAGTCTGCTGATTTTTTTTTTTATTTAATCAGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGTGTGTATAGCCTGAAACAGTACGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCg                                                <  2:879522/149‑1 (MQ=255)
                                                                                            ttCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCg                                                >  1:322559/1‑149 (MQ=255)
                                                                                              cGTACTAACAGTCAGTCAGCTGATTTTTATCATGCTTTCAGCAGTTGCCGGCGGGATCTACACTCAAATGGGAAAGGAACGGGTGTATACCCTGAACCAGTCCGTTCTTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGc                                              <  1:319471/149‑1 (MQ=255)
                                                                                                 aCTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTa                                           >  1:623035/1‑149 (MQ=255)
                                                                                                     aCATTCAGTCTTCTTATTTTTTTTTTTATTTATTCAGTTGCGTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGt                                       <  2:127440/149‑1 (MQ=255)
                                                                                                     aCAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAATCAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATAATGCTGAATATCTGCGTAGCGt                                       >  1:793682/1‑149 (MQ=255)
                                                                                                      cagtcagtCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGTTCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTg                                      <  2:322559/149‑1 (MQ=255)
                                                                                                           attcttttgATTTTTATAATTATTTAATAAGTTTCCTTATTGATATACACTAAAATGGTAAAGGAAAGGGTTTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAAt                                 <  2:138511/142‑1 (MQ=255)
                                                                                                            gtcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAAtt                                <  1:447536/149‑1 (MQ=255)
                                                                                                             tcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTg                               <  1:846948/149‑1 (MQ=255)
                                                                                                               tgctgATTTTTATCATGATTACATCAGTTGCCTGCGGGCTCTACACTCAACTGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATAGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGGGGCAATTGAt                             <  1:872460/149‑1 (MQ=255)
                                                                                                                gctgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATc                            >  2:511789/1‑149 (MQ=255)
                                                                                                                 ctgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCt                           <  2:80712/149‑1 (MQ=255)
                                                                                                                   gATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACTGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTa                         <  1:857670/149‑1 (MQ=255)
                                                                                                                     tttttATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATc                       <  1:265130/149‑1 (MQ=255)
                                                                                                                          aTCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAacac                  <  2:1131418/149‑1 (MQ=255)
                                                                                                                            cATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGc                >  2:552931/1‑149 (MQ=255)
                                                                                                                            cATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGc                >  1:547299/1‑149 (MQ=255)
                                                                                                                             aTGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGATATTGATACTGCTTTTGCGGTGGCAAATATTGCTTAATATCGGCTTAGCGTGGCAATTTATCTTATCAACACgat               >  2:1112890/1‑147 (MQ=255)
                                                                                                                                       ctgTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACg     <  2:1029648/147‑1 (MQ=255)
                                                                                                                                          ttGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACggag  >  2:249697/1‑149 (MQ=255)
                                                                                                                                          ttGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGACTGAAACAGTCCTTTATTGATACTGCTTTTGCGGTTGCAATTATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACggag  >  2:900040/1‑149 (MQ=255)
                                                                                                                                          ttGCCTGCGGGATCAACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGATATTGCGCTGGCAAATATTGATGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACCCGCTAAATCCCACggag  >  2:557164/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
GTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAG  >  NC_000913/875261‑875547

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: ‑