Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913177,1050TG42.0% 80.7 / 54.8 112T174P (ACA→CCA) yadShypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (10/55);  new base G (47/0);  total (57/55)
Fisher's exact test for biased strand distribution p-value = 4.80e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 115 total aligned reads displayed.

CGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGAT  >  NC_000913/177031‑177180
                                                                          |                                                                            
cGGTAGCTTAAGAGCCAAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGtt                                                                             <  1:2934288/76‑1 (MQ=255)
  gTAGCTTAAGATCCCAACAAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTc                                                                           <  1:5628087/76‑1 (MQ=255)
   tAGCTTAAGATGCCACCGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTcc                                                                          <  1:487835/76‑1 (MQ=255)
   tAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCGATAGCGTCACGACCATGCCCATCATACTGGCTGTTTcc                                                                          <  2:1861531/76‑1 (MQ=255)
    aGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTTCa                                                                         >  1:3105187/1‑76 (MQ=255)
     gCTTAAGATGCAAGCGAAGACCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAg                                                                       <  2:3692147/76‑1 (MQ=255)
      cTTAGGAGGCCAACGAATGGCGGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGt                                                                      <  2:523065/76‑1 (MQ=255)
       ttAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTTCCGTg                                                                     >  1:148531/1‑76 (MQ=255)
        tAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCCGTgg                                                                    >  1:4444754/1‑76 (MQ=255)
         aaGATGCCACCGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGt                                                                   <  2:2751809/76‑1 (MQ=255)
          aGATCCCAACAAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTa                                                                  <  1:5778192/76‑1 (MQ=255)
            acgCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACg                                                                <  1:3591608/74‑1 (MQ=255)
            aTGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTTCAGTGGTACg                                                                >  1:3975556/1‑76 (MQ=255)
            aTGCAAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTAca                                                                <  2:3600128/76‑2 (MQ=255)
              gcaaaCGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGa                                                              <  2:1788651/73‑1 (MQ=255)
              gCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGGGTTCAGTGGGACGGa                                                              >  1:1133345/1‑76 (MQ=255)
                cAACGAATAGCCGCCAGCCGAATCAATAGCGTCGCGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGaaa                                                            <  2:4780046/76‑1 (MQ=255)
                  acaaaTAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAt                                                          <  1:2155014/73‑1 (MQ=255)
                  aCGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAt                                                          <  1:792196/76‑1 (MQ=255)
                  aCGAATACCCCCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAt                                                          <  1:2196043/76‑1 (MQ=255)
                    gAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTTCAGGGGTACGGGAAAtgg                                                        >  2:62461/1‑75 (MQ=255)
                     aaTAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgtg                                                       <  1:3885579/76‑1 (MQ=255)
                     aaTAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTATTTCCAGTGGTACGGAAAAtgtg                                                       <  1:2363383/76‑1 (MQ=255)
                      aTAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTTCAGTGGGACGGGAAATgtgt                                                      >  1:3504986/1‑76 (MQ=255)
                       tAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCACTGGTACGGAAAATGTGta                                                     >  2:2133698/1‑76 (MQ=255)
                       tAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGGGTCCAGTGATACGGGAAAAGTGta                                                     >  1:2352189/1‑76 (MQ=255)
                        aGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGtaa                                                    <  1:4962260/76‑1 (MQ=255)
                        aGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTTCCGTGGGACGGAAAAAGTGtaa                                                    >  1:2643122/1‑76 (MQ=255)
                        aGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTCCCGTGGTACGGGAAATGTGgaa                                                    >  2:4507642/1‑76 (MQ=255)
                         gccgccAGCCGAGTCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGtaat                                                   <  2:3995590/76‑1 (MQ=255)
                          ccgccAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTCCAGTGGGACGGAAAAAGTGtaata                                                  >  2:4381907/1‑76 (MQ=255)
                           cgccGGCCGAATCAATAGCGTCCCGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGtaataa                                                 <  2:106900/76‑1 (MQ=255)
                           cgccAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGtaataa                                                 <  2:1580674/76‑1 (MQ=255)
                           caccAGCCGAATCAATAGAGGCACGACCATGCCCATCACACTGGCTGTTTCCAGTGGTACGGAAACTGTGtaataa                                                 <  2:2232068/74‑1 (MQ=255)
                            gccAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTTCCGTGGGACGGGAAATGTGgtataaa                                                >  1:657297/1‑75 (MQ=255)
                            gccAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTACAGTGGGACGGAAAATATGtaataaa                                                >  1:4880110/1‑75 (MQ=255)
                            gccAGCCGAATCAATAGCGTCACGACCATGACCATCATACTGGCTGTTTCCAGGGGGACGGGAAATGTGGAATAAg                                                >  1:260278/1‑76 (MQ=255)
                             ccAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGc                                               <  1:3542611/76‑1 (MQ=255)
                             ccAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGc                                               <  1:5341342/76‑1 (MQ=255)
                             ccAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGc                                               <  2:1027144/76‑1 (MQ=255)
                             acAGCCGAATCAAGAGCGGCACGACAAGGCCCATCAGACTGGCTGTTGCCAGGGGTACGGAAAAGGTGTAATAAGc                                               <  2:1063261/75‑1 (MQ=255)
                              cAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTCCAGTGGTACGGGAAATGTGtaataaaac                                              >  1:5209281/1‑73 (MQ=255)
                               aGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTTCAGTGGGACGGGAAATGTGTAATAAGACg                                             >  1:2042706/1‑76 (MQ=255)
                                gCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGt                                            <  2:3687608/76‑1 (MQ=255)
                                gCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGt                                            <  2:896758/76‑1 (MQ=255)
                                gCCGAATCAATAGCGTCACGACCATGCCCATCATACGGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGt                                            <  2:3714604/76‑1 (MQ=255)
                                 ccGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTa                                           <  1:2295062/76‑1 (MQ=255)
                                  cGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAg                                          <  2:5056428/76‑1 (MQ=255)
                                   gAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGc                                         <  2:4645545/76‑1 (MQ=255)
                                   gAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGc                                         <  2:761114/76‑1 (MQ=255)
                                   gAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTTCAGTGGGACGGGAAATGTGTAATAAAACGTAGc                                         >  1:1080321/1‑76 (MQ=255)
                                    aaTCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTCCGGACAATGTGTAATAAGCCGTAGCg                                        <  1:1472854/76‑1 (MQ=255)
                                    aaTCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTTCAGTGGGAAGGAAAAAGGGGAATAAGACGGAGCg                                        >  1:5745845/1‑76 (MQ=255)
                                    aaTCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTCCCCGGGGGCGGGAAAAGTGTAATAAGCCGCAGCg                                        >  2:4937459/1‑76 (MQ=255)
                                    aaTCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTCCAGTGGTACGGGAAATGTGGAATAAGCCGTAGCg                                        >  1:4135780/1‑76 (MQ=255)
                                    aaTCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTCCAGTGGGACGGGAAATGTGTAATAAGCCGCAGCg                                        >  2:4767598/1‑76 (MQ=255)
                                     aTCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGt                                       <  1:5812872/76‑1 (MQ=255)
                                     aTCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTCCCGTGGGGCGGGAAATGTGTAATAAGCCGGAgggg                                       >  2:2439270/1‑73 (MQ=255)
                                     aTCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTCCCCTGGGACGGAAAATGTGTAATAAGCCCTAGCgg                                       >  2:4635126/1‑75 (MQ=255)
                                      tCAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTTCAGTGGTACGGGAAATGTGTAATAAACCGTAGCGTg                                      >  1:3809528/1‑76 (MQ=255)
                                       cAATAGCGTCACGACCATGCCCATCATACTGGCTGGTTCCCGTGGTACGGGAAATGTGTAATAAACCGGAGCgggg                                     >  2:3857663/1‑76 (MQ=255)
                                        aaTAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGa                                    <  2:5290382/76‑1 (MQ=255)
                                        aaTAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGa                                    <  1:394099/76‑1 (MQ=255)
                                        aaTAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGGACGGAAAATGTGTAATAAACCGTAGCGTggg                                    >  2:2349802/1‑75 (MQ=255)
                                        aaTAGCGTCACGACCATGCCCATCATACTGGCTGGTTCCAGTGGTACGGAAAATGTGTAATAAAACGTAGCGTGGa                                    >  1:3252928/1‑76 (MQ=255)
                                        aaTAGCGTCACGACCATGCCCATCATACTGGCTGGTTCCAGTGGTACGGAAAATGTGTAATAAAACGTAGCGTGGa                                    >  1:3195964/1‑76 (MQ=255)
                                         aTAGCGTCACGACCATGCCCATCATACTGGCTGGTTCCAGTGGTACGGAAAAAGTGTAATAAGCCGCAGCGCGGAc                                   >  2:701172/1‑76 (MQ=255)
                                           aGCGTCACGACCATGCCCATCATACTGGCTGGTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGGGGGCaa                                 >  1:3557048/1‑76 (MQ=255)
                                           aGCGGCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACaa                                 <  2:4888422/76‑1 (MQ=255)
                                              gTCACGACCATGCCCATCATACTGGCTGGTTCCAGTGGGAAGGAAAATGTGTAATAAGGCGGAGCGTGGGCAAAAc                              >  1:4065227/1‑76 (MQ=255)
                                              gTCACGACCATGCCCATCATACTGGCTGGTTCCAGGGGTACGGGAAATGTGTAATAAACCGTAGCGGGGGAAATAc                              >  1:1190888/1‑76 (MQ=255)
                                               tCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACAATAcc                             <  2:3134472/76‑1 (MQ=255)
                                               tCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACAATAcc                             <  1:4007385/76‑1 (MQ=255)
                                               tCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACAATAcc                             <  1:2449605/76‑1 (MQ=255)
                                                cACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGGAAATGTGTAATAAGACGTAGGGTGGACCATAccc                            >  2:526718/1‑75 (MQ=255)
                                                  cGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACAATAccgcc                          <  1:3089568/76‑1 (MQ=255)
                                                  cGACCATGCCCATCATACTGGCTGGTTCCAGTGGGACGGGAAATGTGTAATAAGCCGTAGCGGGGGCAATAccgcc                          >  1:1293247/1‑76 (MQ=255)
                                                  cGACCATGCCCATCATACTGGCTGGTTCCAGTGGGACGGAAAATGTGTAATAAACCGTAGCGTGGACAATAacgcc                          >  1:2545219/1‑76 (MQ=255)
                                                  cGACCAGCCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACAATAccgcc                          <  2:2792473/76‑1 (MQ=255)
                                                   gACCATGCCCATCATCCGGGCTTTTTCCAGGGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACAATAccgccg                         <  2:3654937/76‑1 (MQ=255)
                                                   gACCATGCCCATCATACTGGCTGGTTTCAGTGGGACGGAAAATGTGTAATAAGCCGTAGCGGGGGGAATAccgccg                         >  1:1065240/1‑76 (MQ=255)
                                                   gACCATGCCCATCATACTGGCTGGTTCCCGTGGTACGGAAAATGTGTAATAAGCCGTCGCGTGGACAAAAacgccg                         >  2:1813775/1‑76 (MQ=255)
                                                   gACCATCCCCTCCATGCGGGCTGTTTCCAGTGGTACAGAAAAGTTGTAAAAAGCCATAGCGTGGACAATAccgccg                         <  2:5604530/76‑1 (MQ=255)
                                                    cccATGCCCATCATGCGGGTTTTTTCCAGTGGTACGAAAATTGTGAAATAAGCCGTAGCGTGGACAACACCGCCGa                        <  1:2471552/75‑1 (MQ=255)
                                                    aCCATGCCCATCATACTGGCTGGTTTCCGTGGTACGGGAAAAGTGTAAAAAACCCTAGCGGGGGCAAAAccgccca                        >  1:266936/1‑74 (MQ=255)
                                                      cATGCCCATCATACTGGCTGGTTCCCGTGGGACGGGAAATGTGTAATAAACCGTAGCGGGGGCAAAAACGCCGaca                      >  1:2135129/1‑74 (MQ=255)
                                                      cATGCCCATCATACTGGCTGGTTCCAGTGGGACGGGAAATGTGTAATAAACCGTAGCGTGGACAATACCGCCGata                      >  1:3531572/1‑76 (MQ=255)
                                                       aTGCCCATCATACTGGCTGGTTCCAGTGGTACGGGAAATGTGTAATAAGCCGTAGCGTGGGCAATACCGCCGataa                     >  1:4652616/1‑76 (MQ=255)
                                                            catcatACTGGCTGGTTTCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACAATAACCCCGctaataacg                >  2:2839963/1‑73 (MQ=255)
                                                            catcatACTGGCTGGTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACCATACCGCCGATAATAccg                >  2:3081385/1‑74 (MQ=255)
                                                              tcatACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGGGGACAATCACCCCGCAAATACACGc              >  2:603773/1‑76 (MQ=255)
                                                               cgtACGGGCTGTTTCGAGGGGTAAGAAAAATGTGAAATAAGCGGGAGGGTGGAAAATCCCGCCGATAATACAGGCa             <  2:4452529/74‑1 (MQ=255)
                                                                atACTGGCTGGTTCCCGTGGGACGGGAAATGTGTAATAAGCCGTAGGGTGGGCAATACCCCCGATAATACAGGGAg            >  1:5531091/1‑76 (MQ=255)
                                                                atACTGGTGTTTTCAGGGGGTAGGAAAAATGGAAAAAAAGCCGGAGGGGGGACAATACCGCCAATAATACAGGCAg            <  2:1097352/76‑1 (MQ=255)
                                                                  aCTGGCTGGTTCCCGTGGTACGGGAAATGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAggt          >  1:4130580/1‑74 (MQ=255)
                                                                       cTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGGTGCgca     >  1:343107/1‑74 (MQ=255)
                                                                        tGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACAATAccg                            >  1:2999112/1‑52 (MQ=255)
                                                                        tGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACAATAccg                            <  2:2999112/52‑1 (MQ=255)
                                                                         gTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGGTGCGCAGa   >  2:4250134/1‑76 (MQ=255)
                                                                          tttCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGAt  >  1:5078770/1‑76 (MQ=255)
                                                                          |                                                                            
CGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGTAATAAGCCGTAGCGTGGACAATACCGCCGATAATACAGGCAGTTGCGTAGAT  >  NC_000913/177031‑177180

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