Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,899,169 C→T 6.2% A248V (GCG→GTG)  yoaD → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase, regulator of cellulose production

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,899,1690CT 6.2% 206.7 / 3.1 64A248V (GCG→GTG) yoaDputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase, regulator of cellulose production
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base T (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

TGGGCTGGTAGATTCCCTTCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCTTTGCTTAATGCGCGCAGCCAG  >  NC_000913/1899026‑1899311
                                                                                                                                               |                                                                                                                                              
tGGGCTGGTAGATTCCCTTCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTTCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGtg                                                                                                                                            <  3:69238/148‑1 (MQ=255)
      ggTAGATTCCCTTCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGTGGTgctgct                                                                                                                                       >  6:260240/1‑147 (MQ=255)
      ggTAGATTCCCTTCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCAGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATAACCGGTGGTGGTgatgct                                                                                                                                       >  2:155395/1‑147 (MQ=255)
       gTAGATTCCCTTCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGAACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTa                                                                                                                                      <  1:102514/147‑1 (MQ=255)
           aTTCCCTTCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGttt                                                                                                                                 <  1:190107/148‑1 (MQ=255)
           aTTCCCTTCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGttt                                                                                                                                 <  1:380663/148‑1 (MQ=255)
           aTTCCCTTCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGtt                                                                                                                                  >  2:302910/1‑147 (MQ=255)
           aTTCCCTTCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCGGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGttt                                                                                                                                 <  2:350597/148‑1 (MQ=255)
            ttCCCTTCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGttt                                                                                                                                 >  8:234132/1‑147 (MQ=255)
            ttCCCTTCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTa                                                                                                                                >  6:175389/1‑148 (MQ=255)
            ttCCCTTCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCGGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGttt                                                                                                                                 <  2:356151/147‑1 (MQ=255)
               ccTTCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCAGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACtg                                                                                                                             <  4:203060/148‑1 (MQ=255)
                  tCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCg                                                                                                                                               <  4:236774/127‑1 (MQ=255)
                  tCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCg                                                                                                                                               >  3:236774/1‑127 (MQ=255)
                  tCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACtggtg                                                                                                                          >  2:68637/1‑148 (MQ=255)
                  tCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACtgg                                                                                                                            >  6:335459/1‑146 (MQ=255)
                   ccGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACtggt                                                                                                                           <  5:335459/146‑1 (MQ=255)
                      cAACCTGAAGACAATGAAAACTATCAGGTTTCATCGCAACGCTTTACTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCCCTATCTTCAAACACAATTACCGCTGGCGGAGCTGCTAAGTTTACTGGTGGGc                                                                                                                       >  2:164207/1‑147 (MQ=255)
                              aGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGGTGGCGGTGCTGCTTAGTTTACTGGTGGGCTACATCGcc                                                                                                              <  1:164207/148‑1 (MQ=255)
                              aGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGcc                                                                                                              <  5:175389/148‑1 (MQ=255)
                              aGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGcc                                                                                                              <  2:216825/148‑1 (MQ=255)
                                   aTGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCctggct                                                                                                         <  7:234132/148‑1 (MQ=255)
                                        aaCTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGa                                                                                                    <  2:91715/148‑1 (MQ=255)
                                              cAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCAAAAACAATTACCGATGCCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCt                                                                                               >  5:170999/1‑147 (MQ=255)
                                                    tcttcGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTATGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGa                                                                                         <  8:52316/147‑1 (MQ=255)
                                                                        tACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAaat                                                                     >  4:64574/1‑147 (MQ=255)
                                                                            aTTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAaatcaatc                                                                >  8:311647/1‑148 (MQ=255)
                                                                             ttAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCt                                                               <  7:311647/148‑1 (MQ=255)
                                                                             ttAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCt                                                               >  5:135687/1‑148 (MQ=255)
                                                                              tAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTcc                                                                             <  1:182484/133‑1 (MQ=255)
                                                                              tAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTcc                                                                             >  2:182484/1‑133 (MQ=255)
                                                                                           ccGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAAc                                                 <  1:302910/148‑1 (MQ=255)
                                                                                            cGGGGGCTACGGCGCTGGAATGGCACTATCTTGCAACACAATTACCGCTGGCGGTGCTGCTATGTTTACTGGTGGGCTACATCGCCTGTCTGGCGACCGATTACCGAATGAGCTTTTACATCGAAATCAATATGGGCCTGGCGCAACa                                                >  4:39708/1‑148 (MQ=255)
                                                                                               gggCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGttt                                                                                                                                 >  8:263654/1‑64 (MQ=255)
                                                                                               gggCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGttt                                                                                                                                 <  7:263654/64‑1 (MQ=255)
                                                                                                  cTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGtt                                          >  2:355158/1‑148 (MQ=255)
                                                                                                    aCGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTgg                                                             <  5:267301/127‑1 (MQ=255)
                                                                                                    aCGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTgg                                                             >  6:267301/1‑127 (MQ=255)
                                                                                                     cGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGa                                       <  1:68637/148‑1 (MQ=255)
                                                                                                        cgcTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAAtt                                    <  1:366291/148‑1 (MQ=255)
                                                                                                              catggcaCTATCTTCCATCACAATTACCGCTGGTGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGGAACATGAGTTCGAATTGTTCTg                              <  1:155395/148‑1 (MQ=255)
                                                                                                              catggcaCTATCTTCCAACACAATTACCGCTGGTGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTg                              <  5:260240/148‑1 (MQ=255)
                                                                                                              catggcaCTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTg                              <  3:64574/148‑1 (MQ=255)
                                                                                                                 ggcaCTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGGGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGt                             >  1:112290/1‑146 (MQ=255)
                                                                                                                 ggcaCTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGt                             >  8:211194/1‑146 (MQ=255)
                                                                                                                   caCTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGc                         >  6:92025/1‑148 (MQ=255)
                                                                                                                     cTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCt                       >  7:18226/1‑148 (MQ=255)
                                                                                                                      tATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCtt                      >  7:110184/1‑148 (MQ=255)
                                                                                                                       aTCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCttt                     >  8:99219/1‑148 (MQ=255)
                                                                                                                       aTCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCttt                     >  1:278088/1‑148 (MQ=255)
                                                                                                                       aTCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCttt                     >  1:95543/1‑148 (MQ=255)
                                                                                                                       aTCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCttt                     >  7:125740/1‑148 (MQ=255)
                                                                                                                       aTCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCttt                     >  7:118082/1‑148 (MQ=255)
                                                                                                                        tCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCtt                      <  2:95543/146‑1 (MQ=255)
                                                                                                                            ccAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCTTTGCt                  <  6:135687/146‑1 (MQ=255)
                                                                                                                             cAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCTTTGCTTaa               >  7:38465/1‑148 (MQ=255)
                                                                                                                             cAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCTTTGCTTa                <  6:170999/147‑1 (MQ=255)
                                                                                                                               acacAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCTTTGCTTAATg             <  5:92025/148‑1 (MQ=255)
                                                                                                                                cacaATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCTCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACa                                                <  7:295390/112‑1 (MQ=255)
                                                                                                                                cacaATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCTCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACa                                                >  8:295390/1‑112 (MQ=255)
                                                                                                                                   aaTTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCTTTGCTTAATgcgcg         <  2:278088/148‑1 (MQ=255)
                                                                                                                                      tACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAAc                                                 >  7:324726/1‑105 (MQ=255)
                                                                                                                                      tACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAAc                                                 <  8:324726/105‑1 (MQ=255)
                                                                                                                                          gCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCTTTGCTTAATGCGCGcagccag  >  4:252604/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                              
TGGGCTGGTAGATTCCCTTCCGCAACCTGAAGACAATGAAAACTATCAGGTTTCTTCGCAACGCTTTCCTTTTACCATTAACGTTAATGGTCCGGGGGCTACGGCGCTGGCATGGCACTATCTTCCAACACAATTACCGCTGGCGGTGCTGCTAAGTTTACTGGTGGGCTACATCGCCTGGCTGGCGACCGCTTACCGAATGAGCTTTTCCCGCGAAATCAATCTGGGCCTGGCGCAACATGAGTTCGAATTGTTCTGTCAGCCTTTGCTTAATGCGCGCAGCCAG  >  NC_000913/1899026‑1899311

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