Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,516,689 T→C 100% T382A (ACG→GCG)  yfeA ← putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,516,6890TC98.6% 449.8 / ‑3.4 137T382A (ACG→GCG) yfeAputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (69/65);  total (71/66)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

Only 100 of 138 total aligned reads displayed.

AGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGC  >  NC_000913/2516592‑2516787
                                                                                                 |                                                                                                  
  tCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACgg                                                                                                <  1:1047129/100‑1 (MQ=255)
   cGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATAAAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGc                                                                                               >  1:1071340/1‑100 (MQ=255)
    gCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCa                                                                                              >  1:1944328/1‑100 (MQ=255)
     cccTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCaa                                                                                             >  2:1550701/1‑100 (MQ=255)
      ccTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCaaa                                                                                            >  2:3257685/1‑100 (MQ=255)
      ccTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTAGTTTTCCTGCATCAACGGCAGCAGCCCACGGCaca                                                                                            >  2:2922069/1‑98 (MQ=255)
           cGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAg                                                                                       >  1:1149254/1‑100 (MQ=255)
           cGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAg                                                                                       >  1:1808351/1‑100 (MQ=255)
            gTTTCCGGCCCTCTCAGCACTAACAGCAGTTCATTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGc                                                                                      <  1:902843/100‑1 (MQ=255)
            gTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTTATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGc                                                                                      <  2:1777680/100‑1 (MQ=255)
            gTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGc                                                                                      <  1:3029899/100‑1 (MQ=255)
            gTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCATGCATCAACGGCAGCAGCGCACGGCAAATTGAGc                                                                                      >  1:809103/1‑100 (MQ=255)
             tttCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCg                                                                                     <  2:2483631/100‑1 (MQ=255)
              ttCCGGCCCGCTCAGCACTAACAGCAGTTTATTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCgg                                                                                    <  1:3323640/100‑1 (MQ=255)
               tCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGCCAGCAGCGCACGGCAAATTGAg                                                                                       >  2:2972982/1‑96 (MQ=255)
                ccGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGAt                                                                                  >  1:1313465/1‑100 (MQ=255)
                ccGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGAt                                                                                  <  2:979809/100‑1 (MQ=255)
                ccGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGAt                                                                                  <  1:396369/100‑1 (MQ=255)
                 cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                                 >  2:3821409/1‑100 (MQ=255)
                 cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                                 <  2:3865157/100‑1 (MQ=255)
                 cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                                 >  2:1162858/1‑100 (MQ=255)
                 cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                                 >  1:1480660/1‑100 (MQ=255)
                 cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                                 <  1:1691800/100‑1 (MQ=255)
                  ggCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATAc                                                                                >  2:2858603/1‑100 (MQ=255)
                  ggCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATAc                                                                                >  2:106435/1‑100 (MQ=255)
                   gCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACa                                                                               >  2:160686/1‑100 (MQ=255)
                    cccGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACa                                                                               <  2:3836375/99‑1 (MQ=255)
                      cGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCATCGCACGGCAAATTGAGCGGATACAGTg                                                                            <  2:1830524/100‑1 (MQ=255)
                      cGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTg                                                                            <  1:2183179/100‑1 (MQ=255)
                       gCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGa                                                                           <  1:911677/100‑1 (MQ=255)
                        cTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGaa                                                                          <  1:206825/100‑1 (MQ=255)
                        cTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTATCCTGCATCAACGGCAGCAGCGCACGGCAAATAGAGCGGATACAGTaaa                                                                          >  2:2447469/1‑100 (MQ=255)
                          cAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcg                                                                        >  1:2394577/1‑100 (MQ=255)
                          cAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcg                                                                        <  2:3593146/100‑1 (MQ=255)
                           aGCACTAACAGCAGTTCACTACCAGGCAATTGATACAACTTTTCGTTTTCCTGCATTAACGGCAGCAGCGCACGGCAAATTGATCGGATACAGTGAAcg                                                                        <  2:670337/99‑1 (MQ=255)
                            gCACTAACAGCAGTTCACTACCTGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                      <  1:1940622/100‑1 (MQ=255)
                            gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                      >  2:2566549/1‑100 (MQ=255)
                            gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                      <  2:349233/100‑1 (MQ=255)
                            gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                      >  2:1159600/1‑100 (MQ=255)
                            gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                      >  1:125567/1‑100 (MQ=255)
                             cACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACgcgc                                                                     >  1:1249416/1‑100 (MQ=255)
                              tcTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGca                                                                    <  2:2600224/99‑1 (MQ=255)
                               cTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcat                                                                   >  1:3162417/1‑100 (MQ=255)
                               cTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcaa                                                                   >  1:2600248/1‑99 (MQ=255)
                                tAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatc                                                                  >  2:2686877/1‑100 (MQ=255)
                                tAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatc                                                                  >  1:3635973/1‑100 (MQ=255)
                                 aaCAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatca                                                                 >  1:966413/1‑100 (MQ=255)
                                      cagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTGCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAg                                                            <  1:1018155/100‑1 (MQ=255)
                                       agTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACTACAACAACACAAAAAAAATTGAACGGATACAGTGAACGCGCATCATTAAGc                                                           >  1:353469/1‑100 (MQ=255)
                                       agTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc                                                           >  1:2045003/1‑100 (MQ=255)
                                       agTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc                                                           <  2:3674814/100‑1 (MQ=255)
                                        gTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGcc                                                          >  2:3215098/1‑100 (MQ=255)
                                        gTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGcc                                                          >  1:591373/1‑100 (MQ=255)
                                         ttCACTACCCGGCAATTGATACAACTTTTCGTTTTTCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCg                                                         <  1:675795/100‑1 (MQ=255)
                                          tCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGt                                                        <  1:2320893/100‑1 (MQ=255)
                                           cACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTa                                                       >  2:2895813/1‑100 (MQ=255)
                                           cACTACCCGGCAAGTGATACAACTTTTCGTTGTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTa                                                       >  1:2635611/1‑100 (MQ=255)
                                             cTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCATCTGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAAt                                                     <  1:167744/100‑1 (MQ=255)
                                             cTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAAt                                                     <  1:119145/100‑1 (MQ=255)
                                             cTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAAt                                                     >  2:222782/1‑100 (MQ=255)
                                                cccGACAATGGATACAACTATTGGTATGCCGGTATCGACGGCAGAAGCGCACGGCAAATTGAGCGGATACAGAGAACGGGAATCAATAAGAGGTAATga                                                   >  2:2850917/1‑99 (MQ=255)
                                                 ccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacg                                                 <  2:836607/100‑1 (MQ=255)
                                                 ccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTAAAAACGCATCATTAAGCCGTAATgacg                                                 >  1:2694948/1‑100 (MQ=255)
                                                  cGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacga                                                >  1:892708/1‑100 (MQ=255)
                                                     cAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTTAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACt                                              >  1:442063/1‑99 (MQ=255)
                                                     cAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTc                                             <  2:1808336/100‑1 (MQ=255)
                                                      aaTTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCa                                            <  1:1038548/100‑1 (MQ=255)
                                                        ttGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATa                                          <  1:2025237/100‑1 (MQ=255)
                                                        ttGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGAGCATCATTAAGCCGTAATGACGACTCATa                                          >  1:130698/1‑100 (MQ=255)
                                                           aTACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAAc                                       <  2:148146/100‑1 (MQ=255)
                                                           aTACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAAc                                       <  1:2895839/100‑1 (MQ=255)
                                                                aCTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAg                                  <  1:3854459/100‑1 (MQ=255)
                                                                  ttttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGAt                                >  1:3223453/1‑100 (MQ=255)
                                                                     tCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTAt                             <  1:2873153/100‑1 (MQ=255)
                                                                       gTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCa                           >  1:3775719/1‑100 (MQ=255)
                                                                        ttttCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCaa                          >  1:832333/1‑100 (MQ=255)
                                                                         tttCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAAt                         <  2:1956287/100‑1 (MQ=255)
                                                                          ttCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGTCTCATAAACTCAAGATTATCAATg                        <  2:964914/100‑1 (MQ=255)
                                                                            ccTGAATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATgcg                      <  2:3329039/100‑1 (MQ=255)
                                                                              tGCATCCACCGCAGCAGCCCTCGGCAAATTGAGAGGATAAAGTGAACGCGCATCATTAAGCCGTAATGACGAATCATAAACTCAAGATTATCAAttaac                     >  2:1451790/1‑95 (MQ=255)
                                                                                cATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAGGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcagg                  <  2:535464/100‑1 (MQ=255)
                                                                                cATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcagg                  <  1:492448/100‑1 (MQ=255)
                                                                                 aTCAACGGCAGCAGCGTACGGCAAATTGAGCGGATAAAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggc                 >  2:369139/1‑100 (MQ=255)
                                                                                 aTCAACGGCAGCAGCGCACGGTAAATTGATCGGATATAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggc                 <  1:372172/100‑1 (MQ=255)
                                                                                 aTCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggc                 <  2:2912804/100‑1 (MQ=255)
                                                                                 aTCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggc                 <  2:272012/100‑1 (MQ=255)
                                                                                 aTCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggc                 <  1:3821448/100‑1 (MQ=255)
                                                                                 aTCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATACgcaggc                 >  2:161361/1‑100 (MQ=255)
                                                                                  tCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATAATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggca                >  2:778601/1‑100 (MQ=255)
                                                                                   cAACGGCAGCGGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggcag               >  2:2757677/1‑100 (MQ=255)
                                                                                   cAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggcag               >  2:3848825/1‑100 (MQ=255)
                                                                                     aCGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGgcagca             <  1:2175870/100‑1 (MQ=255)
                                                                                        gcagcaGCTCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCaaaa          <  1:2757701/100‑1 (MQ=255)
                                                                                          agcagcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACTACTCATAAACTCAAGATTATCACTGAGCAGGCAGCAAAAAc        <  2:2127105/100‑1 (MQ=255)
                                                                                          agcagcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAAc        >  1:2804460/1‑100 (MQ=255)
                                                                                           gcagcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAAct       <  2:2232293/100‑1 (MQ=255)
                                                                                            cagcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAActc      <  1:2427545/100‑1 (MQ=255)
                                                                                             agcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAActct     <  2:2476380/100‑1 (MQ=255)
                                                                                             agcGAACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGAAAGAAACAAAAActct     <  2:1161395/100‑1 (MQ=255)
                                                                                                 tACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGc  >  2:2266139/1‑99 (MQ=255)
                                                                                                 |                                                                                                  
AGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGC  >  NC_000913/2516592‑2516787

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