Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,407,859 A→G 63.1% intergenic (‑64/‑191) ydaM ← / → ydaN diguanylate cyclase, csgD regulator/putative Zn(II) transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,407,8590AG63.1% 79.7 / 90.4 95intergenic (‑64/‑191)ydaM/ydaNdiguanylate cyclase, csgD regulator/putative Zn(II) transporter
Reads supporting (aligned to +/- strand):  ref base A (17/18);  new base G (28/32);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

AGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAAT  >  NC_000913/1407761‑1407954
                                                                                                  |                                                                                               
aGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGAt                                                                                                <  1:3095729/100‑1 (MQ=255)
 gTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATc                                                                                               <  2:1038185/100‑1 (MQ=255)
  tAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCa                                                                                              <  2:620427/100‑1 (MQ=255)
  tAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCa                                                                                              >  2:2199623/1‑100 (MQ=255)
  tAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCa                                                                                              <  1:161741/100‑1 (MQ=255)
     gTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGAt                                                                                           >  1:395579/1‑100 (MQ=255)
     gTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGAt                                                                                           >  2:1633333/1‑100 (MQ=255)
       ccAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTTTTTTTTTTTTTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATaa                                                                                         >  1:2113668/1‑100 (MQ=255)
       ccAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATaa                                                                                         >  1:2772506/1‑100 (MQ=255)
       ccAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATaa                                                                                         >  1:3308403/1‑100 (MQ=255)
        cAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAAt                                                                                        <  2:369492/100‑1 (MQ=255)
          gggTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTTTTTTTTTTTTATCGCAAAAAAGCCACGTCATATGGTCAGATAATTc                                                                                      >  1:1092576/1‑100 (MQ=255)
          gggTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTc                                                                                      >  2:986374/1‑100 (MQ=255)
          gggTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTc                                                                                      >  1:3313931/1‑100 (MQ=255)
          gggTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTc                                                                                      <  1:2412516/100‑1 (MQ=255)
          gggTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTc                                                                                      >  2:2144373/1‑100 (MQ=255)
            gTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTg                                                                                    <  1:714298/100‑1 (MQ=255)
            gTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTg                                                                                    <  2:3308431/100‑1 (MQ=255)
            gTATTGAAGTTGTGAGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTAAAAAAAAAAAAAAAAGCCACGTCATATGATCAGATAATTCt                                                                                     <  1:1052239/99‑1 (MQ=255)
              aTTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGAt                                                                                  <  1:720203/100‑1 (MQ=255)
               ttGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATa                                                                                 >  1:2445490/1‑100 (MQ=255)
                tGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATaa                                                                                >  1:667982/1‑100 (MQ=255)
                tGAAGTTGTGCGTAATCATTGAGATCCCGTTATTGTTATCAATTTTTGTTGCTATCCAATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATaa                                                                                <  1:1238037/100‑1 (MQ=255)
                  aaGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATg                                                                              >  1:1696927/1‑100 (MQ=255)
                     ttGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATa                                                                           <  2:520845/100‑1 (MQ=255)
                      tgtgCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAg                                                                          >  1:1909465/1‑100 (MQ=255)
                       gtgCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGa                                                                         <  2:68734/100‑1 (MQ=255)
                         gCGTAATCATTCAGATCCCGTTATTTTGATCATTTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACg                                                                       <  1:3263818/100‑1 (MQ=255)
                          cGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGc                                                                      <  2:1638059/100‑1 (MQ=255)
                          cGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGc                                                                      >  1:2624891/1‑100 (MQ=255)
                              aTCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTAtt                                                                  <  1:1652165/100‑1 (MQ=255)
                               tCATTTAGATCCCGTTATTTTTATCAATATTTGTTGCTATCCGATCAAAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTAttt                                                                 <  2:1034759/100‑1 (MQ=255)
                               tCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTAtt                                                                  <  1:626242/99‑1 (MQ=255)
                                cATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTa                                                                >  1:1743953/1‑100 (MQ=255)
                                cATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTa                                                                >  2:2113841/1‑100 (MQ=255)
                                  ttGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAac                                                              >  1:976834/1‑100 (MQ=255)
                                  ttGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAac                                                              <  1:1039017/100‑1 (MQ=255)
                                  ttGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAac                                                              >  2:41658/1‑100 (MQ=255)
                                    gagaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAacac                                                            >  1:3086224/1‑100 (MQ=255)
                                    gagaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAacac                                                            <  2:1809550/100‑1 (MQ=255)
                                    gagaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAacac                                                            <  1:165740/100‑1 (MQ=255)
                                     agaTCCCGTTATTTTTTTTTTTTTTTTTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACt                                                           >  1:37047/1‑100 (MQ=255)
                                     agaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACt                                                           >  1:551461/1‑100 (MQ=255)
                                      gaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACtt                                                          >  1:1888381/1‑100 (MQ=255)
                                      gaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACtt                                                          >  1:132753/1‑100 (MQ=255)
                                       aTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTc                                                         <  2:2841882/100‑1 (MQ=255)
                                       aTCCCGTTATTTTTATCAATTTTTGTTACTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTc                                                         <  1:2258811/100‑1 (MQ=255)
                                        tCCCGTTATTTTTTTCAATTTTTGTTGCTATCCGATCTCAAAAAAGCCACGTCATATGATCCGATAATTCTGATAATGATAGACGCTATTTAACACt                                                           <  1:2634690/97‑1 (MQ=255)
                                        tCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTca                                                        <  1:2285790/100‑1 (MQ=255)
                                        tCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTca                                                        >  2:2890365/1‑100 (MQ=255)
                                         cccGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTcac                                                       <  1:929595/100‑1 (MQ=255)
                                          ccGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTcaca                                                      <  1:728973/100‑1 (MQ=255)
                                          ccGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTcaca                                                      >  2:2237832/1‑100 (MQ=255)
                                            gTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACg                                                    <  2:1909478/100‑1 (MQ=255)
                                             ttATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACgg                                                   <  2:2495428/100‑1 (MQ=255)
                                                tttttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGttt                                                >  1:2331480/1‑100 (MQ=255)
                                                 ttttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTCATAATGATAGACGCTATTTAACACTTCACACGGTTTg                                               <  1:1809663/100‑1 (MQ=255)
                                                  tttATCAATTTGTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGttt                                                <  1:1674627/98‑1 (MQ=255)
                                                   ttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGta                                             >  2:1514615/1‑100 (MQ=255)
                                                       cAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCATATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACg                                         <  1:1983737/100‑1 (MQ=255)
                                                       cAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACg                                         <  2:304690/100‑1 (MQ=255)
                                                       cAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACg                                         >  2:938921/1‑100 (MQ=255)
                                                           ttttGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGaaa                                     <  2:1972361/100‑1 (MQ=255)
                                                            tttGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGaaaa                                    <  2:2444409/100‑1 (MQ=255)
                                                              tgttgCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGc                                  >  2:2004090/1‑100 (MQ=255)
                                                                 tgCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtt                               >  1:957215/1‑100 (MQ=255)
                                                                  gCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAttt                              <  2:976837/100‑1 (MQ=255)
                                                                   cTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtttt                             >  2:85876/1‑100 (MQ=255)
                                                                     aaacgaaaaaaaaaaaGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGc                           <  1:957821/90‑1 (MQ=255)
                                                                     aTCCGATCGCAAAAGAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGc                           >  2:797901/1‑100 (MQ=255)
                                                                        cGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCttt                        >  1:662895/1‑100 (MQ=255)
                                                                        cGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCttt                        >  2:1670940/1‑100 (MQ=255)
                                                                          aTCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCttttt                      >  1:806719/1‑100 (MQ=255)
                                                                           tCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTg                     <  2:2829804/100‑1 (MQ=255)
                                                                            cGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGt                    <  1:760591/100‑1 (MQ=255)
                                                                            cGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGt                    >  2:1207594/1‑100 (MQ=255)
                                                                             gCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTa                   <  1:43450/100‑1 (MQ=255)
                                                                             gCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtttt                             >  1:1282577/1‑90 (MQ=255)
                                                                              cAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTAt                  >  2:1317518/1‑100 (MQ=255)
                                                                              cAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTAt                  >  1:1645514/1‑100 (MQ=255)
                                                                                aaaaaGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTc                >  1:2373050/1‑100 (MQ=255)
                                                                                 aaaaGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCa               <  2:2100018/100‑1 (MQ=255)
                                                                                 aaaaGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCa               <  1:2699337/100‑1 (MQ=255)
                                                                                  aaaGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTATCACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCaa              <  1:3139269/100‑1 (MQ=255)
                                                                                  aaaGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCaa              <  1:2113827/100‑1 (MQ=255)
                                                                                   aaGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAAt             >  1:863313/1‑100 (MQ=255)
                                                                                    aGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAAtt            >  1:1010662/1‑100 (MQ=255)
                                                                                      ccACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTa          <  1:2972192/100‑1 (MQ=255)
                                                                                       cACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAg         <  1:293498/100‑1 (MQ=255)
                                                                                        aCGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCTTTTTGCTTTTTGTATTCAATTTa          >  1:1843984/1‑98 (MQ=255)
                                                                                        aCGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAg         <  1:3380526/99‑1 (MQ=255)
                                                                                         cGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTATCACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGAc       <  1:1955308/100‑1 (MQ=255)
                                                                                          gTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGAcn      <  2:403262/100‑2 (MQ=255)
                                                                                          gTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACa      <  2:648723/100‑1 (MQ=255)
                                                                                              tatGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAAt  >  1:3069164/1‑100 (MQ=255)
                                                                                                  |                                                                                               
AGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAAT  >  NC_000913/1407761‑1407954

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