Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,407,859 A→G 18.7% 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,8590AG18.7% 174.7 / 34.7 90intergenic (‑64/‑191)ydaM/ydaNdiguanylate cyclase, csgD regulator/putative Zn(II) transporter
Reads supporting (aligned to +/- strand):  ref base A (38/35);  new base G (11/6);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

GAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTT  >  NC_000913/1407760‑1407956
                                                                                                   |                                                                                                 
gAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGa                                                                                                   >  2:696582/1‑100 (MQ=255)
 aGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGAt                                                                                                  >  1:2029021/1‑100 (MQ=255)
 aGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGAt                                                                                                  >  2:3862085/1‑100 (MQ=255)
     aGTCGAGGAGATTGAAGTTGTGCGTAATCATTGAGATCGCGTTATTTTTATCAATTTTTGTTGAAAAAAAAAAAAAAAAAAGCCACGTCATATGGTCAGa                                                                                              <  2:1897362/100‑1 (MQ=255)
     aGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGa                                                                                              >  1:888663/1‑100 (MQ=255)
     aGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGa                                                                                              >  1:76872/1‑100 (MQ=255)
      gTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGAt                                                                                             >  1:1153458/1‑100 (MQ=255)
       tCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATa                                                                                            >  2:3692460/1‑100 (MQ=255)
         cAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAAt                                                                                          >  2:1655096/1‑100 (MQ=255)
         cAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAAt                                                                                          >  1:407552/1‑100 (MQ=255)
           gggTATTGACGGTATGCGTAATCATTCAAATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAATCCACGTCATATGGTCAGATa                                                                                            >  2:1040672/1‑96 (MQ=255)
           gggTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTc                                                                                        >  1:1573594/1‑100 (MQ=255)
            ggTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCt                                                                                       <  1:3002993/100‑1 (MQ=255)
            ggTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCt                                                                                       <  1:2040961/100‑1 (MQ=255)
                 tGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATaa                                                                                  >  1:2258295/1‑100 (MQ=255)
                  gAAGTTGTGCGTCATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAAt                                                                                 >  1:2713281/1‑100 (MQ=255)
                  gAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAAt                                                                                 >  1:2248654/1‑100 (MQ=255)
                  gAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAAt                                                                                 <  2:42835/100‑1 (MQ=255)
                     gTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATg                                                                                >  2:2302610/1‑98 (MQ=255)
                     gTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATCCTGATAATGAt                                                                              >  2:2302240/1‑100 (MQ=255)
                     gTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGAt                                                                              >  1:1383969/1‑100 (MQ=255)
                     gTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGAt                                                                              <  1:3962607/100‑1 (MQ=255)
                      ttGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATa                                                                             <  1:624970/100‑1 (MQ=255)
                      ttGTGCGTAATCATTCAGATCCCGTTATTTGTATCAATTTGTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATg                                                                                <  2:3593350/97‑1 (MQ=255)
                       tgtgCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCTACGTCATATGATCAGATAATTCTGATAATGATa                                                                             >  2:3567620/1‑99 (MQ=255)
                       tgtgCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAg                                                                            <  2:2452465/100‑1 (MQ=255)
                        gtgCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGa                                                                           >  1:3454359/1‑100 (MQ=255)
                         tgCGTAATCAGTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACTTCATATGATCAGATAATTCTGATAATGATAGAc                                                                          <  2:2724255/100‑1 (MQ=255)
                           ngTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGAAAAATAGCCACGTGATATGATCTGATAATTCTGATAATGATAGACGc                                                                        >  2:2910792/2‑100 (MQ=255)
                           cGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGc                                                                        <  2:2455631/100‑1 (MQ=255)
                              aaTCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATTATCAGATAATTCTGATAATGATAGACGCTAt                                                                     >  2:612066/1‑100 (MQ=255)
                               aTCATTGAGATCGCGTTAGTATGATCAATTTTTGTTGCTATAGGATCGCAAAAAAGCCACGTCATATGAGCAGATAATTCTGATAATGATAGACGCt                                                                       <  2:3963185/97‑1 (MQ=255)
                                tCATTGAGATCCCGTGATTTTTATCAATTTTTTTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTAttt                                                                   >  2:1220611/1‑100 (MQ=255)
                                 cATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTa                                                                  <  1:3489362/100‑1 (MQ=255)
                                 cATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTa                                                                  >  1:90419/1‑100 (MQ=255)
                                 cATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAAGGATAGACGCTATTTa                                                                  >  1:4079355/1‑100 (MQ=255)
                                 cATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTa                                                                  <  1:4105302/100‑1 (MQ=255)
                                     gagaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAacac                                                              >  1:1980859/1‑100 (MQ=255)
                                     gagaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAacac                                                              >  2:3509963/1‑100 (MQ=255)
                                     gagaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAacac                                                              <  1:1815424/100‑1 (MQ=255)
                                     gagaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAacac                                                              <  1:2997298/100‑1 (MQ=255)
                                     gagaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAacac                                                              >  2:2980919/1‑100 (MQ=255)
                                         tCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTca                                                          >  2:869014/1‑100 (MQ=255)
                                         tCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTca                                                          >  1:288863/1‑100 (MQ=255)
                                          cccGTTATTTTTATCAATTTTTTTTTTTTTTCTATCGCAAAAAAGCCACGTCATATGATCATATAATTCTGATAATGATAGACGCTATTTAACACTTcac                                                         >  2:3457385/1‑100 (MQ=255)
                                           ccGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTcaca                                                        >  2:2330644/1‑100 (MQ=255)
                                             gTTATTTTTTTTTTTTTTTGTTGCTATACGATCGCAAAAAAGTCACTTTATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAATACtt                                                            >  2:1611060/1‑94 (MQ=255)
                                                 tttttATCAATTTTTTTTTTTTTTTTATCGCAAAAAAGCCGCGTCATATGATCAGATTATTCTGATAATGATAGACGTTATTTAACACTTcttatttttt                                                  >  2:271121/1‑91 (MQ=255)
                                                 tttttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGttt                                                  >  1:1007606/1‑100 (MQ=255)
                                                  ttttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTg                                                 <  1:3991626/100‑1 (MQ=255)
                                                  ttttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTg                                                 <  2:732199/100‑1 (MQ=255)
                                                  ttttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTg                                                 <  1:2330657/100‑1 (MQ=255)
                                                  ttttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTg                                                 <  1:2847570/100‑1 (MQ=255)
                                                     tATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGtat                                              <  1:216353/100‑1 (MQ=255)
                                                      aTCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGtata                                             <  2:407549/100‑1 (MQ=255)
                                                      aTCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGtaaa                                             <  1:2803791/100‑3 (MQ=255)
                                                        cAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACg                                           <  2:1002148/100‑1 (MQ=255)
                                                                gttgCTATCCGATCGCAAAAAAGCCACGTCATATGATCTGATAATTCTGATAATGATATACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCa                                   <  2:2248641/100‑1 (MQ=255)
                                                                  tgCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtt                                 >  1:2465758/1‑100 (MQ=255)
                                                                  tgCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtt                                 <  2:3219813/100‑1 (MQ=255)
                                                                  tgCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtt                                 >  2:1888437/1‑100 (MQ=255)
                                                                   gCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAttt                                >  1:1413883/1‑100 (MQ=255)
                                                                    cTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtttt                               >  1:2033849/1‑100 (MQ=255)
                                                                    cTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtttt                               >  1:3847723/1‑100 (MQ=255)
                                                                    cTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtttt                               <  2:2391886/100‑1 (MQ=255)
                                                                    cTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtttt                               >  1:3344976/1‑100 (MQ=255)
                                                                    cTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAttt                                <  2:149620/99‑1 (MQ=255)
                                                                     tATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTg                              >  2:212099/1‑100 (MQ=255)
                                                                     tATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAAGGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTg                              >  2:1374978/1‑100 (MQ=255)
                                                                      aTCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTAACACGGTTTGTATACGGAAAAGCATTTTGc                             >  2:2822060/1‑100 (MQ=255)
                                                                            tcacAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTg                       <  2:827914/97‑1 (MQ=255)
                                                                            tCGCAAAAAAGCCCCGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAACAGCATTTTGCttttt                        >  2:124301/1‑99 (MQ=255)
                                                                            tCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTg                       <  1:337195/100‑1 (MQ=255)
                                                                            tCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTg                       <  1:2980937/100‑1 (MQ=255)
                                                                             cGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGt                      <  2:1607830/100‑1 (MQ=255)
                                                                              gCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTa                     >  2:1505759/1‑100 (MQ=255)
                                                                              gCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTa                     <  2:3847698/100‑1 (MQ=255)
                                                                               cAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTAt                    >  1:1308566/1‑100 (MQ=255)
                                                                                aaaaaaGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTAtt                   <  1:3567643/100‑1 (MQ=255)
                                                                                  aaaaGCCACGTCATATGATCAGATAATTCTGATTATGATAGACGCTAGTTATCACTTCACACGGTTTGTATACGGAAAAGCAAAATGCttttccgtatac                 <  2:3966651/100‑9 (MQ=255)
                                                                                      gCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAAttt             >  2:2256498/1‑100 (MQ=255)
                                                                                        cACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAg           <  1:502243/100‑1 (MQ=255)
                                                                                        cACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAg           <  2:2430262/100‑1 (MQ=255)
                                                                                         aCGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGa          <  2:1222412/100‑1 (MQ=255)
                                                                                           gTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGCAAAGCATGTTGCTTTTTGTATTCAATTTAGACa        >  2:1364027/1‑100 (MQ=255)
                                                                                           gTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACa        <  2:1894359/100‑1 (MQ=255)
                                                                                           gTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACa        <  2:314152/100‑1 (MQ=255)
                                                                                           gTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACa        >  2:1058646/1‑100 (MQ=255)
                                                                                             cATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGa      <  1:997388/100‑1 (MQ=255)
                                                                                                 tGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACaaaagt   >  2:2014874/1‑97 (MQ=255)
                                                                                                 tGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAAttt  <  1:3733914/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
GAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTT  >  NC_000913/1407760‑1407956

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: