Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913301,4100CT57.6% 309.1 / inf 513E130K (GAG→AAG) yagTputative xanthine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (109/107);  new base T (165/129);  total (274/239)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-02
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 538 total aligned reads displayed.

CATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCC  >  NC_000913/301335‑301485
                                                                           |                                                                           
catcatGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGAtctc                                                                             <  1:3473876/76‑1 (MQ=255)
catcatGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGAtctc                                                                             <  1:2387376/76‑1 (MQ=255)
 atcatGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCg                                                                            <  1:5677509/76‑1 (MQ=255)
 atcatGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCg                                                                            <  1:2786103/76‑1 (MQ=255)
   catGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGc                                                                          >  2:3096551/1‑76 (MQ=255)
     tGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGccc                                                                        >  1:2000709/1‑76 (MQ=255)
      gCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGcccc                                                                       >  1:1503747/1‑76 (MQ=255)
      gCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGcccc                                                                       >  2:2968595/1‑76 (MQ=255)
      gCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGCACTTCAATGGTGGTGATCTTGGcccc                                                                       >  1:2200287/1‑76 (MQ=255)
        ttGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCCATGGTCGCGATCTCTGcccccc                                                                     >  2:5699579/1‑75 (MQ=255)
        ttGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCt                                                                     <  2:4847563/76‑1 (MQ=255)
        ttGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCt                                                                     >  1:3353160/1‑76 (MQ=255)
        ttGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCt                                                                     <  1:4256612/76‑1 (MQ=255)
          gATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGa                                                                   <  2:3646810/76‑1 (MQ=255)
          gATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGa                                                                   >  1:4644003/1‑76 (MQ=255)
            tAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATg                                                                 >  2:3808077/1‑76 (MQ=255)
             aaaGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCACCTGATGc                                                                >  1:4493079/1‑76 (MQ=255)
             aaaGGACGCCTGCATGGGGTGAAGATGAGCTGGCGAACCCACACCTGCAAGGGGGGTGATCTCGGGCCCCCGATga                                                                >  1:669541/1‑75 (MQ=255)
               aGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCAt                                                              <  2:2808193/76‑1 (MQ=255)
                ggCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATg                                                             <  1:283368/76‑1 (MQ=255)
                  ccgccTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGAc                                                           <  1:5343731/76‑1 (MQ=255)
                  ccgccTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGAc                                                           <  1:3317354/76‑1 (MQ=255)
                   cgccTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACt                                                          >  1:817790/1‑76 (MQ=255)
                   cgccTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACt                                                          <  1:4791286/76‑1 (MQ=255)
                   cgccTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACt                                                          >  2:4418275/1‑76 (MQ=255)
                    gccTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTg                                                         <  2:345323/76‑1 (MQ=255)
                     ccTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGc                                                        >  2:2206609/1‑76 (MQ=255)
                        gCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAg                                                     <  2:5434248/76‑1 (MQ=255)
                         cATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGc                                                    <  1:4890878/76‑1 (MQ=255)
                         cATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGc                                                    <  1:4087314/76‑1 (MQ=255)
                         cATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGc                                                    >  1:5845519/1‑76 (MQ=255)
                          aTGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCCTGACTGCAAGCg                                                   >  2:5574935/1‑76 (MQ=255)
                             gggTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTca                                                <  2:2477058/76‑1 (MQ=255)
                             gggTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTca                                                <  2:5382457/76‑1 (MQ=255)
                             gggTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTca                                                <  1:5027156/76‑1 (MQ=255)
                             gggTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCAGCCCCCTGATGCATGACTGCAAGCGTca                                                >  1:1910498/1‑76 (MQ=255)
                              ggTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTcag                                               <  1:293602/76‑1 (MQ=255)
                               gTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTcagg                                              <  1:1361704/76‑1 (MQ=255)
                                tGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTcaggc                                             <  1:5749513/76‑1 (MQ=255)
                                 gAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGCCTGCAAGCGTcaggca                                            <  2:3337290/76‑1 (MQ=255)
                                 gAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTcaggca                                            >  2:1203331/1‑76 (MQ=255)
                                 gAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTcaggca                                            <  1:5262780/76‑1 (MQ=255)
                                    gATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCaggcaggc                                         >  1:3949270/1‑76 (MQ=255)
                                    gATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCaggcaggc                                         <  2:835931/76‑1 (MQ=255)
                                     aTTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAggcaggca                                        <  2:328273/76‑1 (MQ=255)
                                      ttATCTGGCTAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCAt                                       >  1:1690919/1‑76 (MQ=255)
                                      ttATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCAt                                       >  1:1582871/1‑76 (MQ=255)
                                       tATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCAtt                                      >  1:1755976/1‑76 (MQ=255)
                                       tATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGCCCGGCAtt                                      >  2:393444/1‑76 (MQ=255)
                                       tATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCAtt                                      <  1:3676137/76‑1 (MQ=255)
                                        aTCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTa                                     >  1:1136024/1‑76 (MQ=255)
                                         tCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCAttca                                    <  1:37900/76‑3 (MQ=255)
                                           tGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGc                                  >  1:5373855/1‑76 (MQ=255)
                                           tGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGc                                  <  2:2239893/76‑1 (MQ=255)
                                             gCGAACCCTGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCt                                <  2:5001407/76‑1 (MQ=255)
                                             gCGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGcca                                <  2:1523323/76‑2 (MQ=255)
                                             gCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCt                                >  2:3544875/1‑76 (MQ=255)
                                              cGAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTg                               >  2:2681487/1‑76 (MQ=255)
                                              cGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTg                               >  1:3021093/1‑76 (MQ=255)
                                              cGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTg                               >  2:5264598/1‑76 (MQ=255)
                                               gAACCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCCACCCCCAGCCCGGCATTAAGCCtcc                              >  2:107569/1‑74 (MQ=255)
                                               gAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGTCTGc                              >  1:4359489/1‑76 (MQ=255)
                                               gAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCtcc                              >  2:5217055/1‑74 (MQ=255)
                                                aaCCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCg                             >  1:2569721/1‑76 (MQ=255)
                                                aaCCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCg                             >  2:4242368/1‑76 (MQ=255)
                                                aaCCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCg                             >  1:3157114/1‑76 (MQ=255)
                                                aaCCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCg                             >  2:706721/1‑76 (MQ=255)
                                                 aCCCAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGa                            >  1:74415/1‑76 (MQ=255)
                                                  cccAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGAc                           >  1:938999/1‑76 (MQ=255)
                                                   ccAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGaac                          >  2:2552854/1‑74 (MQ=255)
                                                   ccAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGAcc                          >  1:519083/1‑76 (MQ=255)
                                                    cAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCa                         >  1:1231105/1‑76 (MQ=255)
                                                    cAGGCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCa                         >  2:1920757/1‑76 (MQ=255)
                                                       gCCTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTg                      >  2:4059629/1‑76 (MQ=255)
                                                       gCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTg                      >  2:1464059/1‑76 (MQ=255)
                                                       gCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCACTg                      >  2:2399775/1‑76 (MQ=255)
                                                        ccTTCAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGa                     <  2:4744346/76‑1 (MQ=255)
                                                        ccTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGa                     <  1:2657537/76‑1 (MQ=255)
                                                           tcgaTGGTGGTGATGTTGGCCGCCTGGTGCATGAGTTCAATCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCa                  <  2:1283947/73‑1 (MQ=255)
                                                            cAATGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTca                                                <  1:5136734/45‑1 (MQ=255)
                                                            cAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAACCGTCAGCCAGGCATTAAGCCTCCCACCATCGAcccc                 >  2:139172/1‑74 (MQ=255)
                                                             aaTGGTGGTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCCGCGACCCTTGACCCGc                >  2:1730078/1‑76 (MQ=255)
                                                             aaTGGTGGTGATCTTCGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCCGCGACCATTGACCAGc                >  2:3731346/1‑76 (MQ=255)
                                                              aTGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCa               >  1:2474464/1‑76 (MQ=255)
                                                               tggtggTGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCAc              <  1:3600203/76‑1 (MQ=255)
                                                                ggtggtGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGCCCATTGACCAGCACg             <  2:2063858/76‑1 (MQ=255)
                                                                ggtggtGATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACg             >  1:2090708/1‑76 (MQ=255)
                                                                  tggtgATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGt           <  1:2274146/76‑1 (MQ=255)
                                                                  tggtgATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGt           <  1:4550671/76‑1 (MQ=255)
                                                                   ggtgATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTa          >  2:1784823/1‑76 (MQ=255)
                                                                   ggtgATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCACCACCGtc          >  2:5092402/1‑75 (MQ=255)
                                                                   ggtgATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGCTa          >  2:461635/1‑76 (MQ=255)
                                                                   gaggATCTTGGCCCGCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTa          <  2:1609202/73‑1 (MQ=255)
                                                                    gtgATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTAc         <  2:1986440/76‑1 (MQ=255)
                                                                     tgATCTTGGCCCCCTGATGCATGACTGCAAGCGGCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACa        <  2:2934111/76‑1 (MQ=255)
                                                                      gATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAg       >  2:2325445/1‑76 (MQ=255)
                                                                      gATCTTGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAg       >  1:1020505/1‑76 (MQ=255)
                                                                      gATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGAACGGTACAg       >  2:5819157/1‑76 (MQ=255)
                                                                          tcGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCAGTGACCAGCACGGTACAGGCTc   >  1:5494417/1‑76 (MQ=255)
                                                                           cGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGCACACGCTcc  >  2:3208723/1‑76 (MQ=255)
                                                                           |                                                                           
CATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCC  >  NC_000913/301335‑301485

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