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.T95.1% 241.3 / 6.8 81intergenic (‑181/‑446)yi41/b4279IS4 hypothetical protein/putative transport protein
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base T (39/38);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCT  >  NC_000913/4501106‑4501250
                                                                            |                                                                         
cAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGt                                                                          <  2:2187668/76‑1 (MQ=255)
 aaaCTATCAGCGACAAAAATATGCTGAAGTTGTGACGTCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGtt                                                                        <  2:3887159/76‑1 (MQ=255)
   aCTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGtttt                                                                     >  2:2530548/1‑76 (MQ=255)
     tATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGtttttt                                                                  >  1:3418749/1‑76 (MQ=255)
     tATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGtttttt                                                                  >  1:1821987/1‑76 (MQ=255)
      aTCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAAAGCGTTTTTTTATGTGGATTGCTTGttttttt                                                                 >  1:2013131/1‑75 (MQ=255)
           cTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCt                                                             <  2:3385948/76‑1 (MQ=255)
           cTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCt                                                             >  1:3236616/1‑76 (MQ=255)
            tgcAAAAATATGCTCAAGTGGTGACATCAGGAACAAAACGCGTGTTGTTATGTGGATGGCTTGTTTTTTTGATCTg                                                            <  2:1362706/74‑1 (MQ=255)
            tACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTg                                                            <  1:2701289/76‑1 (MQ=255)
            tACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTg                                                            <  2:1348485/76‑1 (MQ=255)
             aCAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTgg                                                           <  2:4846865/76‑1 (MQ=255)
             aCAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTgg                                                           >  1:4835734/1‑76 (MQ=255)
             aCAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTgg                                                           <  2:4135319/76‑1 (MQ=255)
              cAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGt                                                          >  2:4059032/1‑76 (MQ=255)
               aaaaaTATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGta                                                         >  2:4167053/1‑76 (MQ=255)
                aaaaTATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGtat                                                        >  2:4008139/1‑76 (MQ=255)
                aaaaTATGCTCAATCTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGtat                                                        <  1:899293/76‑1 (MQ=255)
                 aaaTATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGtata                                                       <  1:1966193/76‑1 (MQ=255)
                  aaTATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATaa                                                      >  2:2797256/1‑76 (MQ=255)
                   atatGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAAc                                                     <  2:2374455/76‑1 (MQ=255)
                   atatGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAAc                                                     >  1:2939703/1‑76 (MQ=255)
                    tatGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACa                                                    >  2:3757604/1‑76 (MQ=255)
                    tatGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACa                                                    <  2:1737101/76‑1 (MQ=255)
                      tGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAgg                                                  >  1:2937921/1‑76 (MQ=255)
                      tGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAgg                                                  >  1:4267427/1‑76 (MQ=255)
                         tCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGtat                                               >  1:2245394/1‑76 (MQ=255)
                          cAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGtata                                              <  1:176719/76‑1 (MQ=255)
                          cAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGtata                                              <  1:1484502/76‑1 (MQ=255)
                           aaTTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATaa                                             >  1:1757304/1‑76 (MQ=255)
                           aaTTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATaa                                             <  1:1484852/76‑1 (MQ=255)
                           aaTTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATaa                                             >  1:654433/1‑76 (MQ=255)
                            aTTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATaaa                                            >  2:2777568/1‑76 (MQ=255)
                             tttGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAg                                           >  1:3471326/1‑76 (MQ=255)
                             tttGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAg                                           <  1:4863676/76‑1 (MQ=255)
                              ttGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAgg                                          >  1:4118328/1‑76 (MQ=255)
                               tgtgACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGt                                         <  2:2507259/76‑1 (MQ=255)
                               tgtgACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGt                                         <  2:4057040/76‑1 (MQ=255)
                                gtgACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGta                                        <  2:3892112/76‑1 (MQ=255)
                                 tgACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGtat                                       >  1:2487033/1‑76 (MQ=255)
                                 tgACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGtat                                       >  2:4133768/1‑76 (MQ=255)
                                 tgACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGtat                                       >  2:4127099/1‑76 (MQ=255)
                                      tCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGa                                  <  1:1333958/76‑1 (MQ=255)
                                       cAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGaa                                 <  1:4116848/76‑1 (MQ=255)
                                       cAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTCTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGaa                                 <  2:273347/76‑1 (MQ=255)
                                           aaCAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGc                             <  1:1412135/76‑1 (MQ=255)
                                            aCAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCa                            <  1:3617030/76‑1 (MQ=255)
                                             cAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCaa                           >  1:1357403/1‑76 (MQ=255)
                                             cAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCaa                           >  1:1362364/1‑76 (MQ=255)
                                              aaaaCGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCaaa                          >  1:2039313/1‑76 (MQ=255)
                                              aaaaCGCGTTTTGTGATGAGGTTTGCTTGGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCaaa                          >  2:2589819/1‑76 (MQ=25)
                                               aaaCGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAg                         <  1:576713/76‑1 (MQ=255)
                                               aaaCGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAg                         <  2:4145569/76‑1 (MQ=255)
                                                aaCGCGTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGa                        <  1:4659865/76‑1 (MQ=255)
                                                aaCGCGTTTTGTTATGTGGATTGCTTGTTTTTCTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGa                        <  1:1909582/76‑1 (MQ=255)
                                                  cgcgTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGaaa                      <  1:2738776/76‑1 (MQ=255)
                                                   gcgTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAAt                     <  2:444069/76‑1 (MQ=255)
                                                   gcgTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAAt                     >  1:1413980/1‑76 (MQ=255)
                                                   gcgTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAAt                     <  2:2097361/76‑1 (MQ=255)
                                                   gcgTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAAt                     >  1:3840096/1‑76 (MQ=255)
                                                   gcgTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAAt                     >  2:3509829/1‑76 (MQ=255)
                                                    cgTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATa                    <  1:569565/76‑1 (MQ=255)
                                                    agTTTTGTTCTGTGGATTGCTTGTTTTTTTGATCTGGTATAACTGGTATAAAGGTATACATAACCGCAACGAAATa                    >  1:828094/2‑76 (MQ=18)
                                                     gTTTTGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATAc                   <  2:2599504/76‑1 (MQ=255)
                                                      ttttGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACt                  >  1:1243225/1‑76 (MQ=255)
                                                       tttGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTg                 >  2:4001761/1‑76 (MQ=255)
                                                       tttGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTg                 <  2:3841613/76‑1 (MQ=255)
                                                        ttGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGc                <  1:3790739/76‑1 (MQ=255)
                                                         tGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCa               >  2:86241/1‑76 (MQ=255)
                                                         tGTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCa               >  2:725516/1‑76 (MQ=255)
                                                          gTTATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCaa              <  2:3908723/76‑1 (MQ=255)
                                                            tATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAg            >  1:1121164/1‑76 (MQ=255)
                                                            tATGTGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAg            <  1:2454040/76‑1 (MQ=255)
                                                            tATGTGGATTGCTTGTTCTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAg            >  1:2010157/1‑76 (MQ=255)
                                                            tATGAGGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAg            >  1:1098937/1‑76 (MQ=37)
                                                              tgtgGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGa          <  2:52650/76‑1 (MQ=255)
                                                                tgGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGaaa        >  2:3432534/1‑76 (MQ=255)
                                                                tgGATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGaaa        <  2:3615198/76‑1 (MQ=255)
                                                                 ggATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGaaca       >  1:1458382/1‑74 (MQ=37)
                                                                 ggATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGaaaa       <  2:1465796/76‑1 (MQ=255)
                                                                 ggATTGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGaaaa       >  1:1417431/1‑76 (MQ=255)
                                                                    ttGCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAg    <  2:2540481/76‑1 (MQ=255)
                                                                      gCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCt  <  2:152405/76‑1 (MQ=255)
                                                                      gCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCt  >  1:3652279/1‑76 (MQ=255)
                                                                      gCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCt  >  2:2784695/1‑76 (MQ=255)
                                                                      gCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCt  >  2:2957012/1‑76 (MQ=255)
                                                                      gCTTGTTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCt  >  2:4599688/1‑76 (MQ=255)
                                                                            |                                                                         
CAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCT  >  NC_000913/4501106‑4501250

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