Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,501,186 (T)6→7 intergenic (‑187/‑440) yi41 ← / → b4279 IS4 hypothetical protein/putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,501,1801.T99.1% 345.4 / ‑3.3 107intergenic (‑181/‑446)yi41/b4279IS4 hypothetical protein/putative transport protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (51/55);  total (52/55)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 119 total aligned reads displayed.

ACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCT  >  NC_000913/4501105‑4501250
                                                                            |                                                                       
aCAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTg                                                                          <  1:4802928/76‑1 (MQ=255)
   aaCTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGttt                                                                      >  2:4946826/1‑76 (MQ=255)
    aCTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGtttt                                                                    >  1:3204292/1‑76 (MQ=255)
      tATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGtttttt                                                                  <  2:3641912/76‑1 (MQ=255)
       aTCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGtttttgt                                                                 <  1:194198/76‑3 (MQ=255)
        tCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGtttttttg                                                                <  2:663474/76‑3 (MQ=255)
        tCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGtttttttg                                                                <  2:3258748/76‑3 (MQ=255)
         cAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGtttttttga                                                               >  1:1530926/1‑73 (MQ=255)
         cAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGtttttttga                                                               <  1:4327371/76‑4 (MQ=255)
            cTACAAAAATATGCTCAATTTGTGGCATAAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCt                                                             >  1:1955657/1‑76 (MQ=255)
             tACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTg                                                            <  2:1946600/76‑1 (MQ=255)
              aCAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTgg                                                           >  1:3044571/1‑76 (MQ=255)
               cAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGt                                                          <  1:1917221/76‑1 (MQ=255)
               cAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGt                                                          >  2:4943372/1‑76 (MQ=255)
               cAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGt                                                          >  2:5210253/1‑76 (MQ=255)
                 aaaaTATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGtat                                                        <  1:2207327/76‑1 (MQ=255)
                  aaaTATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGtata                                                       >  2:3282211/1‑76 (MQ=255)
                  aaaTATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGtata                                                       <  1:4450846/76‑1 (MQ=255)
                    atatGCTCAATTTGTGACATCAGTAACAAAACGCGTTGTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAAc                                                     <  2:4757317/76‑1 (MQ=255)
                     tatGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACa                                                    <  1:5189635/76‑1 (MQ=255)
                     tatGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACa                                                    <  1:978817/76‑1 (MQ=255)
                     tatGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACa                                                    >  2:2076417/1‑76 (MQ=255)
                     tatGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACa                                                    <  1:1702519/76‑1 (MQ=255)
                     aatGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGGTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACa                                                    <  1:5219182/75‑1 (MQ=255)
                      atGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAg                                                   >  1:2739857/1‑76 (MQ=255)
                      atGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAg                                                   >  1:1571866/1‑76 (MQ=255)
                      atGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAg                                                   <  1:2834073/76‑1 (MQ=255)
                          tCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGtat                                               <  1:4589022/76‑1 (MQ=255)
                           cAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGtata                                              <  2:636464/76‑1 (MQ=255)
                            aaTTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATaa                                             <  2:2094295/76‑1 (MQ=255)
                             aTTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATaaa                                            <  1:2521763/76‑1 (MQ=255)
                             aTTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATaaa                                            <  1:3512836/76‑1 (MQ=255)
                              tttGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAg                                           <  2:4745765/76‑1 (MQ=255)
                               ttGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAgg                                          >  1:4615052/1‑76 (MQ=255)
                               ttGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAgg                                          >  1:93658/1‑76 (MQ=255)
                                tgtgACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGt                                         <  2:2663757/76‑1 (MQ=255)
                                 gtgACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGta                                        <  2:4810456/76‑1 (MQ=255)
                                  tgACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGtat                                       >  1:4042652/1‑76 (MQ=255)
                                  tgACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGtat                                       >  2:5583261/1‑76 (MQ=255)
                                  tgACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGtat                                       <  2:1818677/76‑1 (MQ=255)
                                   gACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACGGGTATAAAGGtata                                      <  1:1468498/76‑1 (MQ=37)
                                   gACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGtata                                      >  2:4570465/1‑76 (MQ=255)
                                   gACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGtata                                      <  2:1610990/76‑1 (MQ=255)
                                   gACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGtata                                      >  1:3144395/1‑76 (MQ=255)
                                     cATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACa                                    <  2:5746372/76‑1 (MQ=255)
                                      aTCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATCCAg                                   <  2:1548892/76‑1 (MQ=37)
                                      aTCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAg                                   >  1:247624/1‑76 (MQ=255)
                                      aTCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAg                                   >  1:2691001/1‑76 (MQ=255)
                                      aTCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAg                                   >  1:723237/1‑76 (MQ=255)
                                       tCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGa                                  >  2:1369492/1‑76 (MQ=255)
                                       tCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGa                                  <  1:5774047/76‑1 (MQ=255)
                                         aGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGaaa                                >  2:4494524/1‑76 (MQ=255)
                                         aGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGaaa                                >  1:2284653/1‑76 (MQ=255)
                                         aGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGaaa                                >  2:3514460/1‑76 (MQ=255)
                                          gTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGaaaa                               >  2:3140258/1‑76 (MQ=255)
                                          gTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGaaaa                               >  1:2324556/1‑76 (MQ=255)
                                          gTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGaaaa                               <  2:3441300/76‑1 (MQ=255)
                                          gTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGaaaa                               >  1:3430255/1‑76 (MQ=255)
                                          gTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGaaaa                               <  1:3515745/76‑1 (MQ=255)
                                             aCAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCa                            >  1:3837161/1‑76 (MQ=255)
                                             aCAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCa                            >  1:2507986/1‑76 (MQ=255)
                                              cAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCaa                           <  2:2739857/76‑1 (MQ=255)
                                              cAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCaa                           <  2:2066158/76‑1 (MQ=255)
                                               aaaaCGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCaaa                          >  2:5367600/1‑76 (MQ=255)
                                               aaaaCGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCaaa                          >  2:1397450/1‑76 (MQ=255)
                                                 aaCGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGa                        >  2:511162/1‑76 (MQ=255)
                                                 aaCGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGa                        <  1:744684/76‑1 (MQ=255)
                                                  aCGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGaa                       >  2:4066571/1‑76 (MQ=255)
                                                   cgcgTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGaaa                      <  1:3540661/76‑1 (MQ=255)
                                                   cgcgTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGaaa                      <  1:3511324/76‑1 (MQ=255)
                                                   cgcgTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGaaa                      <  1:5367601/76‑1 (MQ=255)
                                                    gcgTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAAt                     <  1:3787447/76‑1 (MQ=255)
                                                     cgTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATAGAGAAAAGCAAAGAAATa                    <  1:271241/76‑1 (MQ=37)
                                                     cgTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATa                    >  1:5610155/1‑76 (MQ=255)
                                                      gTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATAc                   <  1:3368648/76‑1 (MQ=255)
                                                       ttttGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACt                  <  1:1369492/76‑1 (MQ=255)
                                                           gTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGGAAAGCAAAGAAATACTGCaa              >  1:4351685/1‑76 (MQ=37)
                                                            ttATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGACATACTGCaac             >  2:4641313/1‑75 (MQ=37)
                                                            ttATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCaaa             >  2:4634164/1‑76 (MQ=255)
                                                            ttATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCaaa             <  1:3811392/76‑1 (MQ=255)
                                                            ttATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCaaa             <  2:4110577/76‑1 (MQ=255)
                                                            ttATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCaaa             <  1:739944/76‑1 (MQ=255)
                                                            ttATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCaaa             <  2:2562238/76‑1 (MQ=255)
                                                             tATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACGGAAAAGCAAAGAAATACTGc                <  2:2709518/72‑1 (MQ=37)
                                                             tATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACGGAAAAGCAAAGAAATACTGc                >  1:2709518/1‑72 (MQ=37)
                                                             tATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAg            >  1:5833542/1‑76 (MQ=255)
                                                             tATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAg            <  2:4828381/76‑1 (MQ=255)
                                                              aTGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAgg           >  2:1894329/1‑76 (MQ=255)
                                                                gtgGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGaa         <  1:570852/76‑1 (MQ=255)
                                                                 tgGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGaaa        >  1:5764744/1‑76 (MQ=255)
                                                                 tgGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGaaa        >  2:4730161/1‑76 (MQ=255)
                                                                  ggATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGaaaa       <  2:93205/76‑1 (MQ=255)
                                                                    aTTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACa     <  1:4494525/76‑1 (MQ=255)
                                                                    aTTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACa     >  1:2156661/1‑76 (MQ=255)
                                                                    aTTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACa     >  1:1372731/1‑76 (MQ=255)
                                                                     ttGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAg    <  1:4570466/76‑1 (MQ=255)
                                                                      tGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGc   <  1:297066/76‑1 (MQ=255)
                                                                      tGCTTGTTTGTTTGATCGGGAATAACAGGTATAAAGGGATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGc   >  1:4215093/1‑76 (MQ=255)
                                                                       gCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCt  >  2:1272761/1‑76 (MQ=255)
                                                                       gCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCt  <  1:4943373/76‑1 (MQ=255)
                                                                            |                                                                       
ACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCT  >  NC_000913/4501105‑4501250

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